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Summary. — Elementary group theory is shown to provide easily the 
independent cartesian orthogonal components of tensor properties in 
material systems of finite symmetry. The method is applied to second, 
third and fourth order polar and axial tensors of no intrinsie symmetry 
in trigonal and hexagonal systems, and the final results are tabulated. 


1. - Physical properties of matter represent the coupling action of matter 
between influences exerted and effects obtained (VoIGT (1)); however not all 
quantities entering in effect-influence relations are physical properties of matter. 
In this respect one distinguishes between field quantities and matter quantities. 
While matter quantities are in themselves consequences of the intimate structure 
and cohesion of matter, field quantities simply represent fields induced in 
matter by external causes: typical examples are surface stress (matter quan-- 
tity (2)) and elastic stress (field quantity). 

Generally more than one number is required to specify a physical property 
at a point in a material system: whenever this occurs, the property is called 
anisotropic since it depends on direction except for the systems in which for 
this property all directions are equivalent by symmetry. When only one 
number is necessary in any material system, the property does not depend 
on direction and is called scalar. The number specifying a scalar property 


(*) This work was largely performed at the Department of Physics, University of 
Illinois, Urbana, IL, U.S.A., during 1950. 

(1) W. Voigt: Lehrbuch der Krystallphysik (Leipzig, 1928), p. 122. 

(2) R. SHUTTLEWORTH: Proc. Phys. Soc., 63, 444 (1950). 


49 - Il Nuovo Cimento. 
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at any given point transforms as a tensor of order zero of polar or axial nature. 
The numbers specifying anisotropic properties frequently transform as com- 
ponents of polar or axial tensors of order greater than zero: since effects and 
influences themselves can generally be represented by tensors, this occurs 
when the effect-influence relation, if it is not linear, can be expanded in power 
series of the influence. When this cannot be done the description of aniso- 
tropie properties with tensors fails, as it does for cleavage and slip phenomena 
(HERMANN (*)). Here we shall consider only tensor properties of matter (matter 
tensors). 

Time reversal symmetry is responsible for the intrinsie symmetry of matter 
tensors. Indeed the existence of intrinsic symmetry for static properties is 
proved by the method of thermodynamic potentials; and the applicability of 
this method is conditioned by macroscopic reversibility (Voir (4)). Similarly 
the intrinsic symmetry of dynamic properties is expressed by Onsager’s re- 
lations; and the validity of Onsager’s relations follows from microscopic rever- 
sibility (ONSAGER (°)). 

Macroscopic space symmetry is responsible for the additional reduction of 
independent matter tensor components in symmetrical material systems: 
indeed matter tensor components must be invariant under the existing sym- 
metry. Although crystals are usually referred to, this statement applies as 
well to any system. The number of independent components of a matter tensor 
in a system of given symmetry is invariant under coordinate transformations (°), 
when reference is made to a general tensor of given order, type (polar or axial) 


(3) C. HERMANN: Zeits. f. Kryst., 89, 32 (1934). 

(4) See reference (+), sections 101, 102, 103: for a specific example section 166. 

(5) L. ONSAGER: Phys. Rev., 37, 405 (1931); 38, 2265 (1931): see also H. B. G. 
Casimir: Rev. Mod. Phys., 17, 343 (1945); $. R. De Groor: Thermodynamics of Irre- 
versible Processes, Section 2 (Amsterdam, 1951); W. Boas and J. K. MACKENZIE: Progr. 
Metal Physics, 2, 92-94 (1950). 

(6) The writer has been unable to find this theorem in the standard literature on 
matter tensors. The proof is so simple that the theorem actually appears obvious. 
The invariance of matter tensor components implies in any given frame, for any given 
symmetry operation (or group of symmetry operations) the existence of linear rela- 
tionships among the components: the number q of independent relationships equals 
the rank of the secular matrix. Since the rank of this matrix is invariant under mul- 
tiplieation by the non-singular transformation matrix (7), the number of independent 
components r = p—q (p = total number of components) is invariant under coordi- 
nate transformations. 

Besides its theoretical interest, the theorem is useful since it enables the numbers 
of independent components of matter tensors to be found by using the reference 
frames which are most convenient to handle. 

(7) See, for instance, M. Böcner: Kinführung in die höhere Algebra (Berlin, 1910), 
Satz 7, p. 85. 
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and intrinsic symmetry. The » independent components themselves can be 
chosen with a certain amount of freedom in each reference frame: some of them 
can be substituted by others through the solutions of the p— r (p = total 
number of eomponents) independent invariance relationships. Within this 
freedom, they are the same in frames which are not distinguishable in their 
relations to the existing symmetry operations; or which can be obtained one 
from the other by symmetry operations; or finally which differ from one another 
by rotations of any angle around symmetry axes of order higher than that 
of the tensor (HERMANN (®)). 

Several methods have been used to study this effect of macroscopie space 
symmetry on tensor properties of matter: while equivalent theoretically, they 
differ in practical utility. NEUMANN (8) imposes identity to the defining effect- 
influence relations written in frames equivalent by symmetry. Vorer (9), 
Love (7°) and most of the early workers in the field (11) impose invariance 
under the symmetry operations to scalars built out of tensor components. 
POCKELS (©) applies Neumann’s and Voigt’s methods. HERMANN (°) was ap- 
parently the first to impose invariance directly on each tensor component: 
he simplifies the equations of invariance by using as reference frames the prin- 
cipal axes of the symmetry operations, generally complex. WOOSTER (3), 
Bonp (1) and MASON (15) also impose invariance directly on each component: 
WoostER and MAson derive the equations of invariance by tensor algebra, 
and BOND by matrix algebra. Fumı (%) eliminates completely the algebraic 
calculations needed to obtain and solve the equations of invariance with his 
direct inspection method which uses the correspondence between cartesian or- 
thogonal tensor components and coordinate products. In the usual frames 
the method can be applied te obtain directly the independent tensor com- 


(8) F. NEUMANN: Vorlesungen über die Theorie der Elastieität (Leipzig, 1885), 
chap. 12. 

(%) See reference (1), sections SI, 82: for specifie applications, sections 147, 169, 
286 and 414. 
10) A, BE. H. Love: A Treatise on the Mathematical Theory of Elasticity, (Cambridge 
1927), 4th edition: section 105. 

ne See, for instance, M. ABRAHAM: Hne. d. Math. Wiss., IV. 14. III (Leipzig, 
2-08). 
12) > Pockets: Lehrbuch der Krystalloptik (Leipzig, 1906). 
18) W. A. Woosrer: A Textbook on Crystal Physics (Cambridge 1938): chap. I. 
er L. Bonn: Bell System Techn. Journ., 22, 1 (1943). 
15) X ms P. Mason: Bell System Techn. Journ., 26, SO (1947); and Piezoelectric 
cae and their Application to Ultrasonics (New York, 1950). 
(6) F. G. Fun: a) Acta Oryst., 5, 44 (1952); b) Acta Cryst. (1952) to appear. The 
use of the correspondence between the transformation properties of material constants 
and of coordinate products has been advocated lately also by H. JURETSCHRE: Acta 
Oryst., 5, 148 (1952). 
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ponents only for triclinic, monoclinic, orthorhombic, tetragonal and cubic 
symmetries in which the conventional coordinates do not transform into 
linear combinations of themselves under the symmetry operations: symme- 
tries C,;, S, and C,, can be treated directly only by using somewhat peculiar 
frames. However the method enables also to derive the independent com- 
ponents in the usual frames for all trigonal and hexagonal groups from the 
corresponding components for symmetry C;. Direct inspection is by far the 
simplest of the methods mentioned and it eliminates the possibility of errors 
which have been relatively frequent with the other method (see table X (17-1418). 
A significant application of direct inspection has been made to obtain the 
third-order elastic coefficients for all the crystallographic groups (%). 

With the exception of one, the group theoretical treatments of tensor pro- 
perties of symmetrical systems published so far give only the number of inde- 
pendent components: this number is logically taken equal to the number of 
independent totally symmetrical linear combinations of tensor components in 
the symmetry group of interest. To obtain this number BHAGAVANTAM (2) 
sets up the corresponding polar tensor representation: for polar tensors he 
looks then for the number of totally symmetrical components, while for axial 
tensors, he looks for the number of one-dimensional components antisym- 
metric in improper rotations. Each symmetry group is treated separately. 
JAHN (#1) avoids the lengthy character calculations of BHAGAVANTAM by re- 
ducing the polar or axial tensor representations of the full rotation group 
through his formulae for the symmetrical products of representations (22), de- 
rived with Tisza’s expression for the characters (2): he comes down to the 
group of interest through correlation tables. 

Apparently group theory has been applied to list independent components 
of matter tensors only in a paper on elastic constants by JAHN (2). The con- 
dition that the elastic potential energy be totally symmetrical in the sym- 
metry group of the system enables JAHN to write it down in terms of the 
linear combinations of strain components subtending the completely reduced 
strain tensor representation: this gives him the number of elastic constants 
in each symmetry group. The relations between them are found by com- 


(7) F. POCKELS: Ann. d. Phys. Chem., 37, 144 (1889). 

(18) W. P. Mason: Bell System Techn. Journ., 29, 184-86 (1950). 

(12) F. G. Fumi: Phys. Rev., 83, 1274 (1951); 86, 561 (1952). 

(2°) S. BHAGAVANTAM: Proc. Ind. Acad. Sei., 16, 359 (1942); S. BHAGAVANTAM 
and D. SURYANARAYANA: Acta Cryst., 2, 21 (1949). 

(21) H. A. JAHN: Acta Oryst., 2, 30 (1949). 

(22) H. A. JAHN: Zeits. f. Kryst., 98, 191 (1937). The proof of the theorem on 
which the method is based has never been published by Dr. JAHN (private commu- 
nication). 

(73) L. Tisza: Zeits. f. Phys., 82, 48 (1933): Eqs. 13-16. 
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paring the most general expression of the elastic potential in terms of strain 
components with the one obtained. 

The present paper is an attempt to treat the tensor properties of material 
systems of finite symmetry by group theory in a general fashion (**). 


2. — The components in any reference frame of a matter tensor in a material 
system must be invariant under the existing symmetry operations: the same 
holds, of course, for their linear combinations. Then each non-zero component 
(and in particular each independent one) must be invariant if one considers 
its relations with the other non-zero components and the non invariant linear 
combinations of non-zero components must vanish. Similarly each non-inva- 
riant component must vanish and, as a consequence, the same holds for any 
linear combination of them. 

The theory of three-dimensional rotation groups is a very suitable instru- 
ment for identifying the non-zero components in a chosen frame and for then 
obtaining a complete set of independent relationships among them; so leading 
to the independent components. One must simply a) construet all the inde- 
pendent linear combinations of tensor components belonging to the totally 
symmetrical representation; and b) construct all the independent non-totally 
symmetrical linear combinations which can be formed with the components 
contained in the bases a). These latter components are all the non-zero ones: 
the vanishing of the linear combinations b) gives a complete set of independent 
homogeneous equations among them. To find a set of independent compo- 
nents it is generally sufficient to look over these linear equations, but a cri- 
terion is available in case of need: a set of non-zero components leaving after 
cancellation of their columns a non-vanishing minor of the secular matrix is 
a set of independent components. The number of independent components 
obviously equals the difference between the number p — v (v = number of va- 
nishing components) of non-zero components and the number q— v (g ='num- 
ber of independent relationships among all components) of independent rela- 
tionships among them; that is the number p — q of totally symmetrical bases a). 

In this way one can study a matter tensor in a material system of any 
given symmetry. However, the existence of definite correlations among the 
irreducible representations of rotation groups suggests the possibility of ob- 
taining almost immediately the results in all subgroups of a group treated 
with proper completeness. To achieve this end it is necessary «) to construct 
the independent linear combinations of tensor components belonging to all 


(24) The treatment has been extended to infinite groups by Dr. H. J. JURETSCHKE: 
H. J. JUREISCHKE: Application of Group Theory to the Description of Physical Pro- 
perties of Orystals, Notes for Physics 8493, Polytechnic Institute of Brooklyn. Refe- 
rence is given to the present paper of the writer as due to appear. 
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the one-dimensional representations of the group which are correlated to the 
totally symmetrical representation of some of the subgroups; and ß) to con- 
struct with the components contained in the bases «) all the independent 
linear combinations belonging to the other irreducible representations of the 
group which can be formed. For each subgroup, the non-zero components 
are those contained in the bases «) which belong to its totally symmetrical 
representation: the vanishing of all the other linear combinations which contain 
these components gives a complete set of independent relationships among 
them. 

Without loss of generality, we shall restrict ourselves to cartesian ortho- 
gonal reference frames. It is then possible and convenient to replace tensor 
components in the procedure with coordinate products of corresponding order 
and intrinsic symmetry. Indeed polar tensor components in cartesian orthogonal 
frames transform as do coordinate products. Axial tensor components in 
cartesian orthogonal frames, instead, transform as do products of coordinates 
only under proper rotations, while under improper rotations their transform- 
ation equations differ from those for products of coordinates by a factor — 1. 
So, when axial tensors are considered, one must look in a) for linear combi- 
nations of coordinate products belonging to the one-dimensional representation 
antisymmetric in improper rotations; and in «) one will be interested in all 
the one-dimensional representations of the group which are correlated to such 
a representation of some of the subgroups. We shall denote the cartesian 
orthogonal components of polar and axial tensors with the coordinate products 
that are their indices. 

Since we wish to consider only finite symmetry groups, the usual formula 


(1) Ds == 


(n, = number of ith irreducible components of the reducible representation; 
h; = number of elements of class j in the group of order N; y'(&;) = cha- 
racter of an operation of class j in the reducible representation; y,(R,;) = cha- 
racter of an operation of class j in the 7th irreducible representation) provides 
us the number of bases a), or «), which can be constructed with coordinate 
products of a given order and intrinsic symmetry: the characters y'(R;) of 
the polar tensor (coordinate products) representations are quickly obtained 
by general or symmetrical (TıszA (2?)) direct products. 

The numbers of bases b), or ß), are easily found. Inspection of the cha- 
racter tables and consideration of the frame chosen gives the types of bases; 
the corresponding numbers of linear combinations are obtained by subtracting, 
from the numbers of coordinate products of interest, the numbers of bases a), 
or a), already constructed with them. Examples will be given in section 3. 


er 
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TABLE I (*). — Correlations between irreducible representations 


of trigonal and hexagonal groups. 


1 Dan | Dg—62 [o,n-6/m|o,,-6mm @,=6'| Se-8 | Dyn-62m Cyn-6 Dis - 3m Ds-32|C,,- 3m! 05-3 

— I un | Le |— ae) | = | 
A, | ae a AY 1 ERA ae, Amon 410) 

Ui | 1 

A | 1 A u As | A A u A 1 | A 4 lu A 1 A, A 
Agg | A, | A, A LEP Os Pea AE À Ao, A, Ay A 
Ae Ags) AS ATT Am EA, Dalal ed Ace || Als A, A 

| Bo | Bi B, B, BAER AB SAS AA A, A 

| Bu B, 18%, | B; | B Ab A, za A À; 4, A 
Be | A En Re A Wl a ras A 

| Boy, B, | B, | B, B Ay A, Al Aou À, A, A 

| Ex | B- 2, ae Se 0 boa SE E37 LE 2 E E BR! 
Bm ENL FE NEE, Bun Ho PEL E' E, B E E 
Pee ied, Ol Babi eg By HSE E, BE E Di 

El Bis Bu) BU lee Bl, BY BY | Be E E E 

| | 


Where alternative choices of corresponding elements are possible, the follow- 
ing choices have been made 


| Don OF v Dan Dsù D; ( 13, 
| 30, 30, 30, 30, 
| 30, 
3i0, 30, 30, 30, 
310, 30, 30, 


(*) Each row of the table lists the irreducible representations of the trigonal and hexagonal 
groups obtained from the representation of Dg, heading the row by eliminating the symmetry 
elements of Dg, which are not contained in the various groups. The nomenclature is the one 
adopted by EYRING, WALTER and KIMBALL, Quantum Chemistry (New York, 1944). 


The construction of the bases can be done with Wigner’s theorem (?°) by 
summing the products of the characters of the representation of interest with 
the transforms of a coordinate product under the symmetry operations of the 
given group. To apply the theorem, it is not necessary to build complete 
transformation tables for the coordinate products since regularities exist. 
Frequently a simpler procedure can be used: indeed higher order bases are 
obtained by multiplication of bases formed with lower order coordinate pro- 
ducts, when all representations in the direct products, but at most one, are 
one-dimensional. At times Wigner’s theorem is preferred just to obtain simpler 
bases (and so linear equations of more immediate solution), by applying it 
to proper sums or differences of coordinate products: a case will be met with 
in table IX. When no intrinsic symmetry is present, other bases may be 


(25) E. WIGNER: Gruppentheorie (Braunschweig, 1931), p. 123 ff. 
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TABLE II. — Character table of Dg. 
= ——— 9 = | = = 
Ee ail C, 205 2) 20, ye #230, 30, 
+ = | 
| | 
A, 1 | 1 1 I 1 1 | 
A, 1 I 1 I == | =] 
B, | 1 — 1 1 =" 1 | =i 
Be 1 =i | l | =] == À | 1 | 
E, 2 il a =] | 1 ON 0 
B, 2 2 — 1 — 1 0 | 0 


obtained from one constructed with either of the procedures, simply by per- 
mutation. : 

There is an apparently more systematic way of treating tensor properties. 
of matter by group theory: it consists in studying fully the polar tensor re- 
presentation corresponding to the given matter tensor. One applies formula (1) 
to obtain types and numbers of the independent bases of the group which 
can be constructed with all the coordinate products of the given order and 
intrinsic symmetry; then one constructs them, by multiplication or with 
Wigner’s theorem. This is certainly inconvenient whenever a good number 
of coordinate products do not need to be considered, that is when many tensor’ 


TABLE 111. — Transformation table of x, y, 2 in D, 


(0), y//C,(1). 21/06). 


= dE 7 5 | F 
E oat 0: 05 | ot OF 
FE V3 1 Van V3 V3 
i 0? DEAR Enr 3 ne ment rase 
ea RER! AB a Rigi ay 
Ne RAR en Aer eat Zr Zi aaa oe 
| | | | 
| FA | 2 | 2 | 2 oe LA 
oy |  c'{2) 0,8) | 0,2) 0,2) 0,8) 
1 ALS u u | 1 1/3 1 v3 
A — 1 —y =) | = @ — 4 
eg u ea ar, 2 2 7 
V3 er V3 V3 1 V3 1 
N De | } ~ 1 1 = 
7 2 x | 2 TG =e al 7 3 x SE à x ay 
—2 —z =% ea Suet — à — 2 
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components vanish in the group, or in the group and subgroups, of interest. 
Otherwise this method of treatment is practically equivalent to the other and 
may be preferable on the ground of personal taste. Examples of its use will 
be given in section 3. 


TABLE IV. — Polar second order tensor. 
6 Representation | Corresponding Independent non-zero 
rou : | : 
ZEIT of interest | representations of D,, | components | 
| | ev ne een | 
ee 2 ee 
| | | 
| 
Se | A rt Pipes PIE = VU CV EVE 
| | 
Da | Ay, Ay,, By, LE — YY, 22 


One can predict a priori that the tensor does not distinguish (i = inversion) 


beiweenseD, (DUC V Da A), ande DE Oss. Dan; 


mn. 


» OCMC = C3, X 1) > Os, sn: 
» SEO) Dt, 
» DD, x u 0,,x%) » 1D. C85. 


When one knows the totally symmetrical bases for D,, and S,, one can state 
that the tensor does not distinguish also (o,= symmetry plane xy) 


between Dg,(D3qXo,) and Da 


and » O67(S6 X On) » Sty. 


It is worthwhile mentioning a point of importance in practice. Unneces- 
sary work is avoided if one remembers that even coordinate products are even 
under inversion and odd ones are odd: it is also useful to keep in mind the 
even or odd behaviour of coordinate products under some simple reflections 
and rotations. 


3. Applications. 


We choose to study by group theory the tensor properties (up to the forth 
order) of material systems having trigonal and hexagonal symmetry groups, 
all of which are subgroups of D,,. Since even coordinate products are even under 
inversion and odd ones are odd, one can always do the computations neces- 
sary to obtain the x) and ß) bases for D,, = D, x7 in the group D,. 
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In D, we adopt the cartesian orthogonal frame x//C;, y//C5, 2//C,. This 
choice fixes the reference frames in the other groups through the correspon- 
dence of the symmetry elements: the cases where alternative choices are pos- 
sible are settled in the lower part of table I. Attention should be called to 
the eriteria of equivalence of reference frames given in section 1, page 741. 


TABLE V. — Axial second order tensor (2°) 
| Group Rent on Corresponding | Independent non-zero 
| of interest representations of De, components 
ie —— ee AS) 2 

D, Ay | A AN, RN RE 

ee as | Ae Al ry = — yx 

OF | À A An Asem Aes TE YY, Les Ly —— Ye 

D, | A, | A A Bu» By, LL — YY, 2 

Oz A, | Ay, Ang» Bu Boy | LY = —yYy% 

05 A ll ls, | zu = yy, 22; xy = — ya 
| One can predict a priori that the tensor vanishes identically in the centro- 
| symmetrical classes D,,, Ogns Se, Daa- 

When one knows the bases «) of interest in ©, and D, one can state that 
the tensor will vanish also in O,, = C;,Xo», Der = Ds X 6. 


It is necessary to consider in «) all the one-dimensional representations 
of Dj, — 6/mmm since they are all correlated to representations of interest for 
some of the subgroups (table I). The numbers of bases «) which can be 
constructed are provided by formula (1). Since D, does not contain improper 
rotations, we need the expression of the characters y'(R,) only for proper ro- 
tations (#7): for second order tensors of no intrinsic symmetry 


(2) L'(Rig)) = (1 + 2 cos p}? = 1 + 4 cosy + 4 cos?g, 
for third order tensors of no intrinsie symmetry 


(3) L'(R;(g)) = (1 + 2 cos p}° = 1 + 6 cos y + 12 cos? p + 8 cos, 


(5) Despite the difference in transformation behaviour under improper rotations, 
we denote the cartesian orthogonal components of axial tensors with the coordinate 
products that are their indices. 

(27) For improper rotations, the characters are in order 


(2') (Rp) = (— 1 + 2 cos p)? = 1 — 4 cos y +4cos?p, 
{3’) 1 (B,(v)) = (— 1 + 2 cos p} = — 1 + 6 cosy — 12 co p + 8cos?p, 
(4) x (Ri(p)) = (— 1 + 2 cos y)* = 1 — 8 cos p + 24 cos? p — 32 cosy + 16 cost. 
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and for fourth order tensors of no intrinsie symmetry 
(4) x(R;(p)) = (1 + 2 cosp)i= 1 + 8 cos p + 24 cos? + 32 cos?p + 16 cost p 


The results for D,, are in order 24,, + Ag; 34,, + 442, + Bi, + Bau; and 
1041, + 94%, + 4B,, + 4B,,. 

Concerning the bases B), inspection of the character table of D, (table IT) 
and of the transformation table of x, y and z in it (table III) shows that 
the coordinate products contained in A bases enter also Z, bases and only 


constructed 


TABLE VI. — Third order polar tensor representation of Dg, (*)- 
Ay, = Ayo X Aou (xy — yu) xz = aye — Ya, 
TRY — yen, ery — yx 
Aou = Arg X Aou Ce RS 3H, bases 8) can be 
(a? + y)xe = we Se 


2 


DEE + YY, 2% 


a 13 products —7 bases: 


Bi | — y*x — wy” — yxy 8 products — 2 bases: 
Po, X Ly = Bau = Bon Lan 3H, bases ß) can be 
Bau y> — ny — Ya? — XYx constructed 
Eyu = A9 X By, (x? + y?) x (@, y) = (a? + y?a, y® + wy), 
(2° + yay, y® + ayx), (0° + ay, y% + yx?) 
Di HA, (x? — y?, wy + yr) Xe = (82 — y?2, mye + Yue); 
r2 2 


(wen — yay, way + yeu), (en? —ey®, zy + 2yx) 


(*) Of the E bases only those containing -1 and B products are listed. 


those, while the coordinate products contained in B bases enter also H, bases 
and only those. 

Since the chosen z belongs to the representation A,, of Ds, and the chosen 
(æ, y) to the representation H,,, only one of the three bases «) for second 
order tensors can be constructed by multiplication of first order bases: 


Ax, X A = Ax, 5) XS — 2? 8 


The other two must be constructed with Wigner’s theorem: one obtains 
æ+y—A,, and sy—yax-A,,. The five coordinate products .2?, y?,,2?, 
ay, ye entering the three bases «) allow one to form only one E,, base f) 
(i.e. two linear combinations) which must be constructed with Wigner’s theo- 
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rem: E,- (0 — y?, xy+yx). The results for second order tensors in all 
groups considered are given in tables IV and V. 

The treatment and results for third and fourth order tensors are summa- 
rized respectively in tables VI, VII, VIII and in tables IX, X, XI. Table XII 
contains some results on photoelastic constants obtained by particularizing 
table X. 

To treat second, third and fourth order matter tensors in trigonal 
and 


TABLE VII. - Polar third order tensor. 


| Represen- | Corresponding | 
| Group | tation | representations Independent non-zero components 
| of interést of Der 
| Es = = = 
| DE A, see table V | CY = — Yue, 02Y = — Yeu, exy = — 2YD | 
| (Of. A, en Fb: a NT LL NEIN ee — CU 
Ce A | see table V | those of D, and C,, 
Dan A, Ayo: Bau a8 yu SF Ty? I yay | 
Og A’ Ayg, A3gs Baus Boy | À ya ay” yay | 
| | | 
| | y = — ay = — yx? = — aye 
D, A, see table V | those of D, and D | 
Os A, Are Asus By Bau 226022 VER OP, A RUF | 
| ye = — ay = — ya? = — aye | 
On on A see table V | those of D,, Og, and Oz, | 
| The tensor vanishes identically in the centrosymmetrical classes D C 
y À on» Yen» 
Stay ID at 


hexagonal systems by group theory one could also carry out the full reduction 
of the polar tensor representations (2), (3) and (4) of D,. The results for 
Do, are in order: 2A,,+ As, + 2H,, + En, 34x, + AA + By, + Bo, + 6H, + 
+ 32, 104,, + 94,, + 4B,, + 4B,, + 16H,,+ 11P,. Besides the bases 
that we have already met one should construct: two bases Z,, for second order 
tensors 


INGA SET ya IO 2X (a, y) = (24, 2%), (22, yz) ; 
three bases #,, for third order tensors 


ARE 2X (20, 29) = (22%, 22y) 5 (222, 2y2) , (022, ye?) 


PET 
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and finally four bases #,, for fourth order tensors 
A,,XH,,=Hy,, 2X (2m, zy) = (28a, y), (2202, y), (exe, 2y22), (02°, yes). 


All these new bases do not need to be used, since the corresponding polar 
and axial tensor components vanish in D,, and in all its subgroups con- 
sidered. 

The treatment that we have presented of second, third and fourth order 
tensors in trigonal and hexagonal groups has, of course, gained in simplicity 
from the gradual increase in tensor order. To treat directly a fourth order 


TABLE VIII. — Axial third order tensor (2%). 


| | 
| 3 
1 Represen- |  Corresponding | 
1 Group | tation | representations Independent non-zero components | 
I I | I | 
| of interest of Der | | 
: | | | 
ee el) Le tes + i | 
4 | | 
1 | 
| | | 
Dan A | Au Lye = — Yue, TRY = — Yeu, 20y = —2YX 
| | 
| | 
| | NE 
Con LE | AS ASE | “ye = — yar, rey = — Yen, LAY = —TYX 
| | DI VRR UYN, OG? QUE 
Se A, MA, Asus Pins Bo, | 22 = UT, 02) = —yen, 20y = 2x 
| CT 2 =— CYR, DER ey, 20° — CU? 
| | x = — yn = — ry? = — yay 
| | y3 ay = — yo? = EYE 
| | 
i Da Ae | Ae Br LYZ = —yuz, ney = — yeu, 20, = 2YX | 
| | D = — yx = — vy? = — yry | 


The tensor does not distinguish | 


| between ED ande, Oona; 
| C 


2 Cen » Rey ORE 
» S, » OER 
4 » Da » ID Ger: 


tensor without any previous knowledge on second and third order bases it 
would be necessary to use Wigner’s theorem extensively: indeed the 104,,, 
9A,,, 4B,, and 4B,, fourth order bases would come out of direct products of 
first order bases, which do not satisfy the condition that all representations 
in the products, except at most one, be one-dimensional since A, — Ad, 
AR, XE} As SA SOs elas, Die — Aon iia) and PB, A, X Hin 


143 lus 
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Tage IX. — Fourth order polar tensor representation of Dg, (+). 


Ay, = Au X Arg 2x2? = 2A 

(ae an y?) x (x? a y?) =: a + ys + a?y? + Hien 

xt t+ yt + xy?s + yory, wt + yt + ayay + yoyx 

(2? + y®) x 22 = we? + y222, 2202 + 2292, 41 products 
— 19 bases: 
11Z, bases By 
can be 
constructed 


202 + 2y?2, wera + yey, newe + yeyz, PART + yey 


Ane ANA, SEX (wy — yx) = aay —27yx, ty? — yar?, 

vey — y2?x, erye—eyne, weYe— Your, 2ARY — ZYRL 
(x? + y?) x (xy — yx) = ay — wyx — Px + yay, 
ay — yo? — yx + vy?, vey — vyx? —yPa + yoy? 


yr2y— a? zy —yueu—x yeu, yey?—Yyeu*—xeyx—xexy , 32 products 


zy? — 20°y — zyn® — enya — 8 bases: 
| | 12E, bases Bp 
: 2 can be 
ER N Re EE 
; | | ; ee ? : i constructed 


RC, (y8 — x y — yu? — EYE) X2 = y2z—ax*ye— YER— ryre, | 
x20 —Yy*2n—anyey—yury, RR —YELY —YRLYE—RRY? | 


ex? — gy?” — 2UY? — 2YXY 


Br br oui (0? + 0y?, y® + ya?) xe = (we + ay?2, y®e + yu?e), 
(az + yayz, Ye + ayaz), (ae + yrxe, Ye + x? yz), 
(0220 + yea, yrey + xPey), (wen + ayey, yrey + yuan), 
(vn + yazy, yrey + year), (wea + yery, yey + xzyx), 
(xz? + yeyx, yey? + xexy), (wen? + wey, yay? + yeu), 


(0° + z0y?, zy® + zyu), (22? + zyay, ey® + 2xyx), 


(2x? + 2y?xr, zy? + 20%) 
BBA (2? — y?, wy + ya) X 2? = (02? — y?2?, wye? + yxe?), 
( 2 y?, 22ary + 22yx), (exe — zy?r, 2xye + 2yxe), 
(we? — ye?y, yeèx + a2?y), (were — yeye, weye + Yyeue), 
( 


ZELE — eyey, emey + eyeu) 


= Baux Ein (a — ya — xy? — yay) x (@, y) = 
= (a — ya? — xx — yaya, ay — yay — xy? — yoy’) 


Byy X Eu] (@4 — y% wy + ya), (ny — yx, «yx + yPay), 
= (xy?a — yary, veya” + yay?), 
Bou X Eu) (ayy — yayı, ya? + wy?) 


(+) Of the E bases only those containing 4 and B products are listed. 
(*) To obtain simpler bases Wigner’s theorem is used to construct a*—y‘*, a?y?—y'a* and! 
wy + y’r: permutation yields the others. 


use: 


third order tensors distinguish. 


at = yt = Y'a? + aye + yayı, 

Yea? LUE X — Yay, Yaya = TYTY 
Bayer BR ey, RR — ZUR, 
LEN = YRYS, CULL = Lyey 


yz yurZ DYLR = — xyz, 
yey = — URYy = — YORE = — LYRT, 
YRY? = — neyn = — neny = — Yen’, 
zy? = — 2YL? = — 20Y% erty 


The tensor distinguishes those and only those symmetry classes that axial 


warn = yery, 
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TABLE X. — Polar fourth order tensor. 
7 | 
Represen- Correspond- | 
Group tation ing represen- Independent non-zero components 
of interest tations of Don | 
ee ae = =e | 
Don Ay, | Ay | 2 
| at = yt = port + aye + yoyo, | 
ya? = «7 y?, ay? = yay, yoy = ayay | 
De <2, 2: Rey, LE eyes, 0220 = yey, | 
| URE = YRYR, ZUR = ZYeY 
Con A, Ara» Avg at | 
we == YP Ye eye Hyaym, | 
yPa® = ay}, ry?a = yary, yayı = wyay 
SOL = ME, TR N N NEE 
| BERR— Yeye, exer = YU 
Lye*—= — yxe", 2° = — 22%, 20Y2 = — yak, 
YRL— — ney, neye = — Yeue, RARY = — ZYRS 
| ey = — pa = yay + ay? + yary?, | 
yeay = — ayn, oy = — YA, yxy? = — ayn? | 
Ss see table IV | see table IV | 24 | 
| | at = y = 202 | BEYER + YLYL, 
| yta? = ay, eye = yoy, yayı = yxy 
| (DERE I, BEE RN En, el | 
| | ene Yeyen 2ER — ZUR | 
| eye? = — ya2?, Roy = — yo, exye = — RYER, | 
| y2*n = —x2"y, xeye = — yere, RAY = — zyen 
wy = — ya yay + xy + yoy, 
| yxy = — nyr, ay? = — ya, yoy? = — vy? 
| | ype YyB22 Lyne = — “yz, 
y?2y = ay Yxee = — 2YRR, 
| | yay? = — eye LeLy — — Yeu", 
| | 2y?® ya: = — znyr = — 22?y - | 
| | aie = — zy?2 = — yayz = — y*xe, | 
02% yen = — DY2Yy = — yuey, | 
| 220% = — yeny = — yeya = — azy?, | 
| 2x8 zany? yxy = — 2y7a 
Desa see table IV | seetableIV| 24 
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TABLE XI. - Axial fourth order tensor (2°). 
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Represen- |Correspond- 
| Group tation ing represen- Independent non-zero components 
of interest |tationsofD,, 
D; see table V | see table V | 24 
Ct yk — yo + eye - yayn, Le ie D 
wy’n = yrry, yryx = TYTY 
LE Ye, NE Oe CU TR eee 
| URLS — Y2YS, RARE — ZYRY 
| 
CE see table V | see table V xyz” = — yaz?, 22xy = — 2?ya, zuy2 = — 2yxe, 
yen = — TRY, neye = — yeue, ZORY = — zyen 
ay = — pa = yPoy + xy? + yoy, yYay= —a’yR, 
| xy? = — ya, yay? = — xyx?, 
| 
| 
RCE | see table V | see table V | those of D, and O0, 
4 | 
| 
He De Ay Ans By | yz yure = — xy xz xyz, 
yray wey = — yuen = — LYRX, 
yey? = — wey = — LILLY = — yeu, 
2y? = — 2y0? = — 2Zayn = — 20Yy 
| D, | see table V | see table V | those of D, and Ds, 
| On 44 | Au Au» ye = ya?z BY — xyz, 
Dope ley yrzy a?ey YERE = — LYEX, 
| | yey? = — LRYR = — very = — YRL?, 
| | | BY — eyn? = — 200 = — 202% 
| | xz ay? yayz = — Ye, 
| | wen = — Yan = — ayey YERY, 
| 220° YaRy = — Yaya = — x2y?, 
| 203 = — 20y? = — eyxy = — zyx 
C> | see table V | see table V | xyz? = — yauz?, 22ay = — 2?yw, zaye = — zyxe, 
| yen = — ney, weye = — yee, RARY = — 2YRX 
| ay = — yea = yay + xy? + yayı, Pay= — ya, 
LYS = — yx’, yoy? = — yx? 
ae DUR = — yxye = — Y?xr, 
Dan = — y2n = — LYEY = — yuey, 
wee" yeny YRYL = — ney", 
2%° coy? = — 2yay = — 2y°% 
CO; see table V | see table V | those of D,, O,, and Oz, 
The tensor vanishes identically in the centrosymmetrical classes D,,, Ogn; 


| 


SD 
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TABLE XII. — Photoelastie constants in Ci, Cy, Oz, and S,, Cz. 


Photoelastic constants: fourth order polar tensor pq (p,q = 1, 2, 3, 4, 5, 6). 


Cons OF D On | S, > O3 | 

MATE 015 OO 16 ae A or PLAN 16 | 

no gel 1s OO lo nl le le 19 = 14 — 10 Io 

eat eine 835 T0 0 Cert) esi ER IN 
SONT HN Ar 0 A SALON LL 25 AG 

| 0 0 0 — 45 44 0 | —46 46 0 —45 .44 41 | 

oo 00a) 216. 16 0 15, 14-1, (11213) 


The schemes of POCKELS (17), recently reobtained by Bonp (14) and by 
Mason (38), are wrong for the ten classes O,, Su, Can C3, Sos Can Ce, Con, T 
and T,: attention to the incorrectness of Pockels’ constants for these classes 
| was called by Buacavanram (?°). A derivation of the correct schemes for sym- | 
| metries O,, S,, On, T and T, was given in ref. 16 a, section 3 b. 


4. — The basic elements of group theory are a useful tool for treating syste- 
matically tensor properties of symmetrical material systems. They allow one 
to obtain almost simultaneously the results in a symmetry group and in its 
subgroups; and to use results established in a given group for lower order 
tensors to study higher order ones, in the same group. Typical features of 
the treatment in each specific case are the clearcut distinction of vanishing 
and non-vanishing tensor components and the simplicity of the equations 
obtained among the latter: these features are particularly advantageous for 
high order tensors and for complicated groups. à 
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RIASSUNTO (*) 


Si dimostra che la teoria dei gruppi permette di ottenere agevolmente le compo- 
nenti cartesiane ortogonali indipendenti delle proprietà tensoriali dei sistemi materiali 
a simmetria finita. Si applica il metodo ai tensori polari e assiali del secondo, terzo e 
quarto ordine non dotati di simmetria intrinseca nei gruppi trigonali ed esagonali e si 
raccolgono in tabelle i risultati finali. 


(*) Traduzione a cura della Redazione. 


IL NUOVO CIMENTO MOL xe NGO 1 Settembre 1952 


Interazioni non localizzabili - Esempio dell’effetto Compton. (*) 


R. AScori 


Istituto di Fisica dell’ Universita - Torino 


(ricevuto il 9 Giugno 1952) 


Riassunto. — Si sviluppa il metodo proposto da G. WATAGHIN per intro- 
durre nella teoria dei campi il concetto della non localizzabilità delle 
interazioni. Determinati gli invarianti da cui dipendono i fattori di taglio 
che derivano dall’applicazione del metodo, si ricavano espressioni per 
calcolarli. Come esempio viene svolto lo studio dell’effetto Compton. 


1. — In una lettera dello scorso anno (!) G. WATAGHIN ha proposto una 
modificazione alla teoria dei campi, in cui si ammette che il contributo recato 
da ciascun punto dello spazio-tempo all’elemento di matrice corrispondente 
ad una interazione non dipende solo dai valori delle funzioni d’onda nel punto, 
ma da quelli in un dominio che abbraccia il punto; il dominio è assegnato in 
modo invariante e simmetrico rispetto alle particelle incidenti, mediante la 
sua definizione spaziale e temporale in un loro sistema baricentrico (2). Nel 
presente lavoro, supponendo valido il solito metodo delle perturbazioni, si 
mostra la via da seguire per l’applicazione della teoria allo studio di qualsiasi 
processo elettromagnetico, e si esamina in particolare come esempio il caso 
dell’etietto Compton. 


(*) Lavoro estratto dalla tesi di laurea in fisica, discussa il 12 dicembre 1951 presso 
l’Universita di Torino. 

(1) G. Waracuin: Nuovo Oimento, 8, 592 (1951). Vedi anche: Zeits. f. Phys., 88, 
92 (1934) e 92, 547 (1935). Sull’argomento vedi inoltre: J. RAyskı: Phil. Mag., 42, 
1289 (1951); C. MÖLLER e P. KRISTENSEN: Danske Mat. Phys. Acad. (in corso di pub- 
blicazione); C. MÖLLER: Rep. of Int. Conf. of elementary particles (Bombay, dicembre 1950) 
p. 163. 

(2) L’operatore di interazione dipende percid, attrayerso la definizione del sistema. 
baricentrico, dagli impulsi delle particelle interagenti. 
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Ci limiteremo a studiare l’interazione tra onde piane incidenti ed emergenti 
(matrice 8). Usando per semplicita la seconda quantizzazione solo per il campo 
elettromagnetico e riferendoci alla rappresentazione in cui le æ, sono diago- 
nali, chiamiamo y,(a,) le funzioni d’onda relativistiche degli elettroni e A,(®,) 
l’operatore corrispondente ai fotoni nell’immagine di Heisenberg. Ammette- 
remo che i prodotti tra tali funzioni d’onda e l’operatore A, che compaiono 
negli elementi di matrice hamiltoniana corrispondenti a ciascuna delle inte- 
razioni elementari (nodi dei grafici di Feynman) nelle quali si scompone qual- 
siasi processo, debbano essere calcolati con la regola (vedi (4): 


a) fat, | dn, DORE, + AE, nm) DOME, + tag) 0), 


se le particelle incidenti (assorbite) sono il fotone e l’elettrone di indice a 
(di energia positiva o negativa), ed à emesso l’elettrone b; con la regola: 


(1°) Jan. | din, Dia) Pelé, Sie JPN JAE, = Nu) Polé, “à Nu) ’ 


se l’unica particella incidente € l’elettrone a (di energia positiva o negativa), 
mentre sono emessi l’elettrone b e il fotone; D(n,) (e altrettanto D,(n,)) è una 


funzione universale del vettore n,, e tale che [ din,D(n,) ==.1. Poichè ab- 
biamo supposto determinati gli impulsi delle particelle, le particelle incidenti 
(sopra definite) determinano sistemi baricentrici, cioé tali che la parte spaziale 
dell’impulso totale vi é nulla; ammetteremo che in tali sistemi D(n,) abbia, 
rispetto all’origine, simmetria sferica nello spazio, e simmetria nel tempo, e si 
possa quindi scrivere D(n,) = D(n, |9|), dove n=|n|, e si indicano con 
sottolineatura le grandezze relative a un sistema baricentrico. Si suppone poi 
che D(n,) sia sensibilmente diversa da 0 solo quando 7 e n, sono in ordine di 
grandezza minori o uguali a una lunghezza universale l. 

Ponendo nella (1) A(x,) = aexp [—tk,m,], yıl©,) = uw exp [tp,,%,] si ot- 
tiene 


au, exp [— D,,, —#%, + DE] | dn,Din,)' 


“exp = UD z = k,)n,] ay, Di (Nan) ex} [Ps Mary! ® 
7 l I 


Posto 


(2) [am,Din, exp lin 6), | A, Daly) exp lip,n) = Gp), 


(5) Abbiamo posto y(x,) = p*(x,)ya © A(x,) = Aylazu)y„ dove y, sono le matrici di 
Dirac. Inoltre si usa un sistema di unita in cui # = 1, ce = 1 
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(3) Ua EXP [U(— Dy, — Ky, + Pen) E P u Eu) (Poy) = 
= PE) AE) E bar — KE) (Poy) - 


Analogamente nella (1'), posto A(w,)= a* exp [ik,a,], p(@,) = Ws exp [ip,,æ,] 
si ottiene 


(3') Ua * Uy Exp [i(—- Du k, + DER; Zn 
= PalE,A(E,) plE i Day) (Pon, — Ku) : 


Fer l’uso delle formule (1), (1') e (3), (3°) @ indispensabile osservare che 
uno stato intermedio 7, e quindi il suo vettore di propagazione Dir è diverso 
a seconda del sistema di riferimento S in cui si applica il metodo delle per- 
turbazioni. Converremo di scegliere lo stato di vettore p}, che si ottiene in 
un sistema baricentrico delle particelle incidenti, essendo questa Vunica pre- 
scrizione determinata in modo invariante e simmetrico rispetto alle particelle 
incidenti. 

Le (3), (3') sono le formule che useremo per il calcolo di elementi di ma- 
trice. Esse mostrano che effettuare il prodotto con le regole (1), (1’) equivale, 
nel caso delle onde piane, ad applicare i fattori di taglio Ge G,, dipendenti 
dagli impulsi, al risultato del prodotto ordinario. 


2. — La funzione G (e altrettanto la G,), definita dalla (2), dovendo essere 
invariante, conterrà gli impulsi sotto forma di invarianti. E utile vedere quali 
siano questi invarianti nell’ipotesi fatta che D(n,) abbia simmetria sferica 
nello spazio e sia simmetrica nel tempo rispetto all’origine. Si ha infatti pu 
cando la (2) in un sistema baricentrico: 


Gp.) == | din, exp [iPM] Din, |m|) = | de exp [— ip x}: 


- [ an.Den, |ml)exp [ipenol = | dr exp [— ip xn]Ftn, po) . 


dove la F(n,p,) € simmetrica in p, perchè trasformata di Fourier nel tempo 
di una funzione di Ino|, € dipende inoltre solo da 7. La @(p,), essendo tra- 
sformata di Fourier nello spazio di una funzione a simmetria sferica, avrà 
simmetria sferica rispetto all’origine delle p, e dipenderà quindi solo da 
PP) oltre che da |p|. Potremo quindi porre G(p,) = G(p, [Pol). Per 
avere un’espressione utilizzabile in qualsiasi sistema di riferimento, chiamata P, 
la somma dei vettori delle particelle incidenti, si pud osservare che per defi- 
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nizione si ha in un sistema baricentrico di tali particelle P = 0, e quiudi 


(4) |Po| = 


| PP, 


pe Se LR 
DE 2 are = Vp Ze ue 


E questa un’espressione generale degli invarianti da cui dipendono i fat- 
tori di taglio nella ipotesi fatta (4). 

Nel caso in cui aleuni degli impulsi che compaiono nelle (4) corrispondono 
a particelle virtuali, & comodo calcolare |p,| e |p| determinando per ogni 
particella anzichè il vettore di propagazione p;,, un vettore p,, che abbia di 
questo la stessa parte spaziale, e tale che per ogni interazione elementare si 
possa scrivere una relazione del tipo ÿ + Pi, = 0, dove il segno + corri- 
sponde ad emissione e il — ad assorbimento (5). 

Poichè p? = p,, la seconda delle (4) vale anche sostituendo ai vettori di 
propagazione i vettori p,, ora definiti. Per una particella di massa m si ha poi: 


ph = + VE + m3 = + /t ur — dd + M. 


Nel caso della (1) e (3) si ha P, == p,, + k,, e quindi, essendo p, = k, i 
fattori di taglio assumono la forma 


(6) GP at = G( |Pa N Ar k|, | Di = k,|) = 
((Pa = k (Pau jr k,))? P De (Ben Sr ku))? 
10 / E = ; ‘ Pau — Ku < ’ ; g { 
| (Pau Sr ku)? w i) | (Days + k,)? 
(Calter EN pe] 
5 | (Pau re ku)? 1 u 


== G (2 \/ (Paulin)? — Dives \/ (Paul )? = pike i m2 aie y- Dar = Pauley |. 
Pau aE Le de Dal Pas ST ku ar 2 Paul u : Dis IE ku IE 2Paukyu 


“Pap + M? + 


(7) Gy(P, 1) or G;(0, m) ’ 


dove vale il segno — se a è un elettrone, il segno + se & un positrone. 
Nel caso dalla (1’) o (3°) si hanno formule analoghe in cui l’indice a & 
sostituito da quello b. 


(4) Si noti che l’espressione generale ottenuta per i fattori di taglio contiene come 
caso particolare quello importante G(p,p,), usato, per esempio, recentemente da 
J. Rayski: Phil. Mag., 42, 1289 (1951). 

(5) Sono questi vettori che si usano correntemente nella teoria dei campi, per 
esempio nella formulazione di Feynman. 
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Si noti che affinchè le formule ottenute tendano a quelle ordinarie quando 
gli impulsi tendono à 0, & necessario che sia G(0, m) = 1, @,(0, m) =1 (°). 

Osserviamo inoltre che se si scegliessero le funzioni D(n,) e D,(n,) in modo 
che fossero diverse da 0 solo in un dominio spaziale, i fattori di taglio G(p,) 
e G,(p,) dipenderebbero solo dalla parte spaziale p di p,. Quindi nella (6) 
non sarebbe necessaria la distinzione tra il caso di elettroni e quello di posi- 
troni, e nella (7) sarebbe senz’altro G,(p,,,) = 1. Non sarebbe poi necessario 
precisare, per l’uso delle (1), (1’), (3), 
(3') il sistema di riferimento in cui 
applicare il metodo delle perturbazioni 
per ottenere i vettori corrispondenti 
agli stati intermedi, in quanto la parte 
spaziale di tali vettori non dipende 
dal sistema di riferimento suddetto. 


3. — A titolo di esempio appli- 
chiamo il procedimento esposto al caso 
particolarmente semplice dell’ effetto 
Compton. 

Contrassegnando con l’indice a le grandezze relative all’elettrone e al fotone 
incidenti, con b quelle relative all’elettrone e fotone emessi, con € e d rispet- 
tivamente quelle relative all’elettrone dello stato intermedio nelle due moda- 
lita del processo che corrispondono ai grafici di Feynman I e II di fig. 1, occorre 
calcolare l’elemento di matrice 


Fig. 1. - I e II diagramma di Feynman 
per l’effetto Compton. 


/ ; Das ka Al por <We|A 1k 
He H; | Hy = >. 2 = oe = EZ + 3 


PR 5 (Pas 1% |A | pa; ka, 1k) Ya 1%, 1k,| A|, 1k,» 
2 E,— Ba A £ 


dove i vettori di stato sono contrassegnati dalla funzione d’onda elettronica 
e dal numero dei quanti presenti di dato impulso. 

Aseguendo i prodotti con le regole (3), (3), si ottiene l’elemento di ma- 
trice H': 


H'= Hp, HP — Bou) + 
+ HD — Kou) (Pap — Kay) (Don) - 


Per applicare le (6), (7) scomponiano innanzitutto H,, nelle parti Hy, 
e H,,, corrispondenti rispettivamente à stati intermedi elettronici e positro- 


(6) Vedi J. Rayskı: Phil. Mag., 42, 1289 (1951). 
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nici. Si ha poi: py, = py, =m, ki, = k,, = 0; dalla legge di conservazione 
Pay + Kay == Pou + Mur Si ottiene inoltre pk, = Pyukour Poulin = Ponkau SE 
ha quindi: 

2Dau Kay m? 


(8) H'= Hg: Ben ——-- a (0, m) + 
; Vm? + 2DanKan Vim? + 2Daukau F 


je (Bir ce (Par, V(Payk Cou)” p+ en ar 
mM m ; 


», x à 
+ Hin 68 (Free V (Pasion)? + M + Pan | G?(0, m). 


m’ m 


Per avere risultati facilmente utilizzabili è necessario partire dalla cono- 
scenza degli impulsi in un sistema di riposo dell’elettrone incidente, in quanto 
questo coincide con un sistema di laboratorio quando si esegue l’esperienza 
nel solito modo. 


Si ha allora p, — 0, m = m, da cui pk, = Ma; Poules = M, quindi 
, 2mks m? 
H'= (mie: (a= = ÿ Aye 
V2mk, am m: V2mk, + m? 


+ Hi, G2(2k,, VIE + m? — k,) + Hi: @2(2k,, VIÈ + m2 + ky) G(0, m). 


Supposto @,(0, m) = 1, abbiamo per alti impulsi k, del quanto incidente 
(k > M): 


(9) H' = H/G*(V/2mk,, 0) + H,,- G2(2k,, 0) + H;,: @*(2k,, 2h). 


Si vede di qui che l’elemento di matrice H, nel quale l’emissione « segue » 
l’assorbimento & affetto da un fattore di taglio, dipendente dall’energia del 
fotone incidente, che si discosta sensibilmente da 1 solo quando sia vamk, S 
T 1/12140 MeV/c cioè k, S 40 GeV/c. 

Gli elementi di matrice H},. e Hj, nei quali l’emissione « precede » l’assorbi- 
mento, e che sono i soli importanti nel caso di assorbimento di fotoni di alta 
energia con emissione di fotoni di bassa energia, sono invece affetti da fattori 
di taglio dipendenti dall’energia del fotone emesso; non si ha quindi «taglio » 
sensibile nel caso di emissione di fotoni di bassa energia, qualunque sia l’energia 
del fotone incidente purchè sufficientemente grande, mentre si avrebbe «taglio » 
se l’energia del fotone emesso diventasse dell’ordine dei 140 MeV. 


4, — Le regole poste a base del procedimento proposto forniscono un metodo 
applicabile allo studio di qualsiasi processo elettromagnetico. E chiaro che 
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pur di scegliere opportunamente la funzione D(n,) è possibile ottenere la « con- 
vergenza » di tutti i processi del calcolo delle perturbazioni. Il metodo proposto 
presenta il vantaggio di ottenere il risultato in modo sistematico e di essere 
basato su un motivo fisico anzichè su artifici analitici. Si differenzia inoltre 
dagli altri metodi perché il «taglio » viene introdotto in ogni «nodo » dei gra- 
fici di Feynman, anziche ai soli «lati» che danno luogo a termini divergenti; 
segue da cid che esso opera su ogni processo al quale prendano parte particelle 
di sufficiente energia. 

L’esempio dell’effetto Compton mostra che i fattori di taglio possono pro- 
durre risultati assai diversi sugli elementi di matrice corrispondenti alle diverse 
« modalità » con cui pud avvenire uno stesso processo. E notevole in parti- 
colare che la probabilità di emissione di particelle di bassa energia possa non 
essere influenzata sensibilmente dai fattori di taglio, comunque grande sia 
l’energia della particella incidente. 

Nel caso dell’effetto Compton i fattori di taglio non portano evidentemente 
aleuna modificazione nel campo di energie nel quale é possibile l’esperienza: 
lo studio di altri fenomeni poträ permettere un confronto con dati sperimentali. 


Ringrazio il prof. G. WATAGHIN per il costante interessamento e per le 
discussioni avute sull’argomento. 


SUMMARY 


The author develops the method proposed by G. WATAGHIN to introduce the idea 
of non-local interactions in field theory. The invariants have been determined, on 
which the cut-off factors deriving from the application of the method, depend. Express- 
ions have been deduced to calculate such invariants. For an example the method 
has been applied to Compton effect. 
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Summary. — A theory of the polarization effects in the energy loss of a 
charged particle is given, in which the mechanism of the loss is different 
from that of the Fermi theory. This theory leads to a distribution of 
the loss between ionization-excitation and emission of Cerenkov radiation 
which differs considerably from that of the theories on Fermi’s lines, 
although the stopping power is not essentially changed, being somewhat 
less than in those theories. This treatment leads to an increase of the 
direct ionization at impact parameters larger than atomic dimensions 
after the relativistic minimum, the increase being large in the case of gases. 
The Fermi saturation effect of the loss at distances larger than atomic 
dimensions does also exist in the present theory, although due to a dif- 
ferent mechanism. The relativistic increase of the loss in distant inter- 
actions is due largely to an increase of the radii of action for ionization 
and excitation (Bohr-Williams mechanism), but there is also a contri- 
bution of the Cerenkov radiation (Fermi mechanism). The saturation 
arises from a limitation of the increase of the radii of action due to the 
polarization of the medium and to the saturation of the emission of Ce- 
renkov radiation. This theory leads to a modification of the formula of 
Frank and Tamm for the rate of emission of Cerenkov radiation, the 
emission of high energy Cerenkov radiation being considerably reduced. 


Introduction. 


1. — The first analysis of the polarization effects in the theory of ionization 
is due to FERMI (!). FERMI’s work was extended and improved by Wick (?), 


() E. Fermi: Phys. Rev., 57, 485 (1940). 
(2) G. C. Wick: Rie. Scient., 11, 273 (1940); 12, 858 (1941), Nuovo Cimento, 1, 
302 (1943). 
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HALPERN and HALL (?), A. Bour (?) and SCHÖNBERG (5°). With the exception 
of A. Bour’s paper, all the above papers follow the line of FERMI’s analysis. 
The work of FERMI, Wick and HALPERN and HALL had mainly in view the 
saturation of the rate of loss due to distant interactions, although in FERMrS 
paper it is indicated that the Cerenkov radiation is included in the loss. The 
importance of the Cerenkov radiation was first emphasized by A. Bour (4), 
who showed by an intuitive argument, and also by using FERMT’s results, 
that the whole increase of the loss at relativistie energies was accounted for 
by the emission of Cerenkov radiation. At that time this conclusion did not 
attract much attention, but, later on, OCCHIALINI’s measurements (7) of the 
tracks of relativistic particles in photographic emulsions showed that there 
was no increase of the grain density with energy within the experimental 
errors (10%). These measurements have stimulated theoretical work on ioni- 
zation and led Messen and RITSON (8) and SCHONBERG (5) to the same con- 
clusions as A. BOHR, using FERMI’s methods. SCHÖNBERG (5) has shown that 
the conclusion regarding the increase of the loss in the relativistic region was 
largely independent of the actual expression of the dielectric constant used 
in the Maxwell equations for the dispersive medium, but he remarked that such 
a conclusion might not hold in a rigorous quantum treatment and also that 
the spectrum of the Cerenkov radiation obtained from the theory of FRANK 
aud TAMM (?) and FERMI (!) required important quantum corrections. 

The conclusions of the analysis of the polarization effects on Fermi’s lines 
seemed sufficient to explain the experimental data of OCCHIALINI and the 
Bristol group (1°), as well as those of BowEn and ROSER (1"). However, re- 
cently, experimental data were obtained by Voyvopic (12) and GosH, JONES 
and WILSON (1%) which cannot be explained by the theory on Fermi’s lines, 
because it is necessary to assume an exceedingly strong absorption of the high 


(*) 
(*) 
(8) 
(°) 


HALPERN and H. HALL: Phys. Rev., 57, 459 (1940); 73, 477 (1948). 
Boxr: Det. Kgl. Dans. Vid. Sels., 24, n. 19 (1948). 
. SCHÖNBERG: Bull. Cent. Phys. Nucl. Brux., n. 20 (1950). 
. SCHÖNBERG: Nuovo Cimento, 8, 159 (1951). 

(7) P. S. OccHIALINI: Como Congress (Nuovo Cimento, Suppl., 6, 377 (1949)). 

(8) H. Messer and D. M. Ritson: Phil. Mag., 41, 1129 (1950). 

(9) I. FRANK and I. Tamm: Compt. Rend. Ac. Sei. USSR, 14, 109 (1937); I. Tamm: 
‚Journ. of Phys. USSR, 1, 439 (1939). 

(20) P. H. FowLer: Phil. Mag., 41, 169 (1950); U. CAMERINI, P. H. FOowLEr, W. O. 
Lock and H. MUIRHEAD: Phil. Mag., 41, 413 (1950). 

(1) T. BowEn and F. X. Roser: Phys. Rev., 83, 689 (1951). 

(2) L. Voyvovıc: Bristol Conference (1951); E. Pıckup and L. Voyvopic: Phys. 
Rev., 80, 89 (1950). 

(3) S. G. Gosx, G. M. D. B. Jones and J. G. Wırson: Proc. Phys. Soc., 65, 
68 (1952). 
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energy Cerenkov radiation to get the experimental values of the increase of 
the ionization after the minimum. SCHÖNBERG pointed out that the emission of 
the high energy Cerenkov radiation would be improbable, according to quantum 
theory (4), and proposed a modification of Fermi’s methods (5), which re- 
duces considerably the intensity of the high energy Cerenkov bands and gives 
a larger amount of ionization. In this paper we shall work out in detail this 
modification of the Fermi theory and improve the formulas of reference (15). 

In the existing theory, the polarization effects are taken into account by 
assuming that the icnizing particle moves in a classical homogeneous dispersive 
medium and using the solution of the Maxwell equations for such a medium 
given by FRANK and Tamm (?) and Fermi (1). Actually the solution of the 
Maxweil equations is used to compute the energy loss due to interactions at 
distances larger than 10-8 cm, the loss at smaller distances being taken from 
the quantum theory of BETHE and BLOCEH ('). This procedure, originally due 
to FERMI (1), was never satisfactorily justified. It seems that the validity of 
the theorem of WILLIAMS (1?) and BLocH (16), concerning the equivalence of 
the classical and quantum methods of computing the average loss due to: 
interactions at distances larger than atomic dimensions, has been extrapolated 
to the theory including polarization effects. Such an extrapolation is actually 
not justified, specially in the determination of the relative contributions of 
the ionization-excitation and the Cerenkov emission to the total loss. The 
relative contributions depend critically on the way in which the interactions 
between different atoms influence the resonance and non resonance transfers. 
It is shown by the formulas given in references (5) and (6) that the resonances 
are shifted from the oscillator frequencies w, to the zeros w,; of the dielectric 
constant e(w); thus a range of frequencies preceding immediately an oscillator 
frequency wm; becomes a non resonance region in which coherent emission of 
radiation is important. On the other hand, the shift of the resonance from 
w, to &, and the fact that V1 — P? is replaced by V1 — P2e(o), in the Fermi 
theory, cancel completely the relativistic increase of the ionization-excitation 
loss, since V1 — Be(@;) — 1. The shift of the resonance frequency @; is 
approximately (4xne?/3m)(f;lw;), n being the number of electrons per unit vo- 
lume, e and m the charge and the rest mass of the electron and f; the oscillator 
percentual (5 f;= 1). In the neighbourhood of the path, the polarization 


(4) M. SCHÖNBERG: Nuovo Oimento, 9, 210 (1952). 

(15) M. SCHÖNBERG: Nuovo Cimento, 9, 372 (1952). 

(6) H. BETHE: Zeits. f. Phys.; 76, 293 (1932); Ann. d. Phys, 5, 325 (1930); 
F. Boca: Ann. d. Phys., 16, 285 (1933); Zeits. f. Phys., 81, 363 (1933); C. MÖLLER: 
Anm. d. Phys., 14, 531 (1932); E. J. WILLIAMS: Proc. Roy. Soc., A 135, 108 (1932). 

(17) E. J. WILLIAMS: Proc. Roy. Soc, A 139, 163 (1933). 
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effects are small and the total field is nearly the Lorentz transformed Coulomb 
field, hence the times of collision are essentially the same as in the Williams 
theory. The short times of collision do not allow such a sharpness of definition 
of the resonances as the classical theory requires, in order to separate the 
shifted resonance @, from the peak of the corresponding Cerenkov band, be- 
cause of the uncertainty principle. At distances of the order of R=V me2/4ane? 
the polarization effects are strong and the times of collision become much 
larger than in the Williams theory, so that at such distances the classical 
theory may be applied. It will be shown that the amount of high energy 
Cerenkov radiation which remains, when polarization effects are neglected up 
to distances of the order of R, is rather small. 

It results from the preceding considerations that, at distances from the 
path much smaller than R, the interaction between different atoms may be 
completely neglected, because the polarization is not large and the clean cut 
separation of the frequency ranges relevant for ionization-excitation and Ce- 
renkov emission on both sides of an oscillator frequency @, is not allowed by 
the uncertainty principle. The adequate treatment in this interval of distances 
is therefore the Bethe-Bloch-Möller-Williams theory, which is actually equi- 
valent to a classical theory with the oscillator model, in what regards the 
value of the average loss at distances larger than 1078 cm. 

Following strictly the lines of reference (15), we should apply the Bethe- 
Bloch-Möller-Williams theory up to the distance R, and a modified form of 
the Fermi theory at distances larger than R. This procedure is not very 
accurate, because the polarization effects are already considerable at distances 
somewhat less than R. It is, however, possible to introduce a correction to 
the results of the Bethe-Bloch-Möller-William theory for distances less than R, 
in order to take into account the non negligible polarization effects. Instead 
of solving the Maxwell equations for a dispersive medium by assuming the 
dispersive medium to occupy the entire space, as was done by FERMI, we shall 
first assume that there are no electrons up to the distance R and determine 
the solution of the Maxwell equations for a dispersive non homogeneous me- 
dium in which the dielectric constant has the value 1, at distances less than R, 
and the Lorentz value e(w) at distances larger than R from the path. Thus 
we neglect the effects of the polarization at distances less than R in the com- 
putation of the field at distances larger than R. Our solution of the Maxwell 
equations allows us to determine with sufficient accuracy the effect of the 
polarization at distances larger than R on the field at distances less than R, 
and thus to correct the rate of loss obtained with the Bethe-Bloch-Möller- 
Williams methods. 

In the theory developed in this paper, most of the relativistic increase of 
the loss comes from distances less than R, in which the polarization effects 
play a minor part, and it is mainly due to an increase of ionization and exci- 
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tation. The increase of loss due to interactions at distances larger than 10-8 cm 
is essentially due to the increase of the radii of action, i.e. to the Bohr-Williams- 
mechanism limited to distances of the order of À. The upper limitation of 
the radii of action leads to saturation, whereas the indefinite increase of the 
radii, that appears when polarization effects are neglected, leads to the inde- 
finite logarithmie increase of the Bethe-Bloch-Möller-Williams theory. The 
contribution of the distances larger than R includes the loss due to Cerenkov 
radiation and is small. This part of the loss also gets saturated, by the same 
mechanism as in the theory on Fermi’s lines. The total rate of loss is often 
smaller in the present theory than in the corresponding one on Fermi’s lines, 
there is also a very different distribution of the loss between ionization-exci- 
tation and Cerenkov radiation. The results of the present theory agree satis- 
factorily with the existing experimental data. 

The solution of the Maxwell equations for the field created by the ionizing 
particle and the electrons lying at distances larger than R from the path is 
derived in section 2. The value of the loss due to transfers at distances larger 
than R is computed in section 3. The modification of the formula of Frank 
and Tamm is given in section 4. The contribution of the different oscillator 
frequencies to the ionization is discussed in sections 5 and 6. 

The modifications of the Bethe-Bloch-Möller-Williams formulas due to the 
limitation of the impact parameters to R are given in section 7. The cor- 
rections due to the action of the field of the atoms lying at distances larger 
than R on the atoms at distances less than R are given in section 8. The 
quantum effects connected with the short times of collision are discussed in 
sections 9 and 10. 

The application of the theory to the analysis of the experimental data. 
of GosH, WILSON and JONES and of Voyvopic is given in section 12. 


Modifieation of the Fermi solution of the Maxwell equations. 


2. — In Fermi’s theory, the energy transferred to the medium at distances: 
larger than atomic dimensions is computed with the field determined by the 
Maxwell equations for a dispersive non magnetic medium of dielectric con- 
stant &(w): 


19 1 0H 

— — &pE = rot H— 4zj, - — =—rotE, 
(1) ce dt e Ot 

div gop == 4j , divH=0, 


w denotes a circular frequency and the operator &op is: 


, 0) 
(2) fon = & ( à ; 
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E and H are the electric and magnetic field vectors, respectively. j, and ej 
are respectively the charge and current densities corresponding to the ionizing 
particle of charge z,e moving with constant velocity v along the x axis: 


(3) jy = Aed(æ — vt) (y) da) , fash 


We shall use, instead of Fermi’s fields E” and H®, fields which satisfy 
the Maxwell equations for a non dispersive medium at distances from the 
path less than R 


L Ee ne Ne (= Vy+#<R), 
(4) € ot ; ce dt = 
7e \/ mc 
div E* = 47), , div H* — 0 | = a) ; 


and the Maxwell equations for a dispersive medium at distances from the path 
larger than R: 


1 oH" 
es TO, ae Tot | 

(5) ce dt C dt (e Ss. R) E 
div &pE" —(, dvH”—0, | 


We shall denote by C, the cylinder of radius o with axis on the path. We 
must fit the solutions of (4) and (5) at the cylinder C,. It is convenient to 
use cylindrical coordinates (0,9, %), @ denoting the azimuthal angle around 
the path. By taking into account the well known formulas 


I 2) 10U, DU. 
: je ees "ne Q ec, 
(6) div U 2 9 (20e) Boop ea 
1 - U 
Got); =a eee, 
0 dx OL 
(7) 
OU, Oe ‘ MT ou ns dU, 
(rot U), Er de , (rot trees (oU,) Mp ’ 


we get from (4) and (5) the following boundary conditions for the fields: 


(8) Hi Ende, Bh = ER, B= EY, HH". (0 = R)- 
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1 
(9) E 1 _ grad 1 H — rot À. 


We shall denote by (Aj, A’) the potentials inside C, and by (45 A") the 
potentials outside C,. The equations for the potentials are: 


1» + 124 ep 
(10) ER Ab, A} er Sar Fin 40, 
Sa EP: 2 
en (rt, marne 
(11) The jt AY, 4 0, hab! 


Because of the geometrical symmetry of the problem, the vector potentials 
are parallel to the path 


(12) A, = A,=0, 
and the Lorentz conditions can be satisfied by taking 


(13) At == 45, A == GAP, 


since the potentials depend on # and ¢ only through the combination «#— vi. 
The boundary conditions for the field (8) are satisfied, by taking: 


0A, dA; 
(14) (1 — B)eop Az = — Peon) Az ; RE TEE, 


v 
— — Y — R, Seh é 
do  d f ) 


In the case of an absorbent medium, the condition that the fields vanish 
at infinity shows that we must take: 
5. {ae 
io(* — ) do, 
a Ne 


(16) kw) = 7 (1 — Be()) , k(— ©) = k*w), — 


+00 
(15) Pe = i B(w)K,(ok(w)) exp 


TT 
9 


<arg kw) < 3 


wla 


for w>0., 


B(w) is determined by the boundary conditions (14). A} must be of the 


IONIZATION AT RELATIVISTIC ENERGIES AND POLARIZATION EFFECTS 771 


form 
+00 
17) A — = | oo) exp jie (? = ) dw + 
: +co : 
ve | O(e)Lo(oko(w)) exp lia 2) dak (te) = lel fi #) 


The K and / are modified eylindrical functions: 
(18) K,(u) = = exp [ion /2]H (iu) , I,(u) = exp [— inv/2]J,(in) . 


C(@) is also determined by the boundary conditions (14). The first term in 
the right hand side of (17) is the x-component of the vector potential of the 
field created by the ionized particle moving in vacuum. 

We get from (14): 


i 1= 
(19)  B(o)K,(Rk(w)) = Ku = > Olw)I(Rk, (0) ' Ta)’ 


(20) B(o)k(w)K,(Rk(o)) = LH (Ro) = x Co) (Ro) kolo) . 


By taking into account the well known relation 


(21) Kuh u) + Ki (hu) =, 
we get from (19) and (20): 
(22) B(o) = R|K,(Rk(w))Lo( Rko(@)) kw) + 
1 — P?e(o) 1 
OEL ER) , 
2 | 1— PPe(o) = 
(23) C(w) = in U Fc) K,(Rk(w)) B(w) — Ko(Rko(w))| { To Rko(w)) Y* | 


By taking into account the expressions 


U 2 1 
(24) L(u)æ1, Liu) & 9) Kılu) = jo K;(u) Sy 
(24a) Gul< I, log y = 0,577.:.); 
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and treating R as a variable parameter, we get: 


(25) lim B(w) =1. 


R— 0 
We shall denote by k.(w) the limit of k(w) as v—c. It is easily seen that: 


(26) lim Bio) = LE (Rh (Oo) _ RO) Ho)e o)| | | 


v—>C 


The fields outside C, are: 


+00 
1216 1 . (a 

27) Hy; =——,| B —— — 82) K,(ok ef iby = m, 
( ) : | (@) Es ß | o(@ (®)) exp |ıw É ) ms (60) Bop KY 

+00 
(28) [Dp = = B(w)e"(w) k(@) Ky (ok(a)) exp liw 2 i) do = Bopk®, 

+co 
(29) A? = AY = B(o)k(w)K,(ok(o)) exp \io (2 en i) do = BopH®. 
—6o 
: À 

(30) Bo=ıB (i 4 Ë 


The fields (E", H'') may be obtained from the Fermi fields by the application 
of the operator Bop- 

The fields inside C, are: 
(31) E=E®-E, H = H° + H'. 


E© and H are the fields created by the ionizing particle in vacuum, they can 
be obtained from the Coulomb field by Lorentz transformation: 


po = pe 
(32) BY = —=* a — pt) x 
+ co 
ee Be a CI _ __ &e(1 — PP) — vi) 
x | Hotton exp |io F od — fe 7) + @— ope? 


eS ER 
(33) m VISE 


+ co 


x i Kylek(w)) exp 


| mee tle pe 
io(® old [ed Pe) FU)? 


(34) HO HO BE. | 
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(E', H’) are the corrections to the field of the particle moving in vacuum due 
to the polarization outside C,: 


+00 


(35) E, = = (1 — ps | C(o)L(ok,(o)) exp ie E = i odo, 

: => 2 
(36) E, = — = VAE | C(o)i(oko(®)) exp |iw (= ) lo|do , 
(37) H'=H,=ßR,. 


The field E’ is finite at the position of the partiele at any instant and allows 
us to get immediately the rate of loss per unit length of path due to distant 
interactions, since it is the difference between the total field and the Lorentz 
transformed Coulomb field: 


+ oo 

1W 210 2 

(38) 48 = — z (v-E'(part.)) = — ae (1 — B?) IK (o)o do . 
2e) 


By distant interactions we mean always those at distances larger than R, 
whereas usually the distant interactions are taken at distances larger than 
10e 


3. — We shall now compute explicitly the value of the integral in (38) by 
taking v — c. It results from (23) that: 


dW, 22e? D'OR »\Ko( Rk(w)) b 
3 eurer ; B(o)io d 
(39) dx nv? x | Cas p en (w)io Oy 


(R — real part). 
We shall consider only the limit for » +c, because for low velocities the 


contribution of the interactions at distances larger than R to the loss are small. 
It results from the second equation (24) and (26) that: 


| aw, dn Kl (Rk.(o)) ED d 
(40) lim (— aa ae a [| +4 (Rk«o)) (Rk.(®)) 7) eee 


(41) - k.(0) = lim kw) — © V1— ea), (00 5 <a y <3). 
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The dielectric constant is given by the Lorentz formula: 


FRI ares 
; = Bi De Pa 
(42) A) = DER rs, : (> Br ya 
(see I Be 2a m 
3 7 vi — v? — ig;v 


The »; are the oscillator frequencies measured in units « and the f,; the cor- 
responding oscillator percentuals. The damping constants g are positive and 
very small, when there are no conduction electrons, as we assume. It is 
shown in Appendix I that in the case of a complex e(w) the quantity within 
the square bracket in (40) has zeros in the neighbourhood of the frequencies w, 
and lying in the lower part of the complex w-plane, but has no zeros or singu- 
larities within the quadrant of the positive real and positive imaginary axes. 
We shall deform the integration path and replace it by the positive imaginary 
axis and a quarter of circle at infinity, between the positive imaginary and 
the positive real axes. We shall then neglect the damping constants. The quan- 
tity inside the square bracket in (40) is real on the imaginary axis, so that 
the only contribution comes from the quarter of circle at infinity. The asymp- 
totic expressions of the dielectric constant e and of k,(w) are: 


Agra? Agro? > 
TO age Tr? OR 


(43) e(o) ~ 1 — 


, C 
0) w? @ 


Finally we get: 


(44) lim 


v—c 


( oe = a Le ul wer | 2rnèet 


dæ R? 7 DR mer 


Modification of the formula of Frank and Tamm. 


4. — We shall now examine more in detail the nature of the loss due to 
distant interactions, in order to separate the ionization-exeitation from the 
loss due to emission of Cerenkov radiation. The energy transmitted to the 
medium at distances larger than o, o being larger than R, per unit length of 
path is: 

il AE r 
(45) aA = zal [E" x H"], dS — > | pue ER AS. 


Co 6, 


The integral is taken over the surface of the cylinder ©, and n denotes the 
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normal component. It is easily seen that: 


+ 69 


aw À 
TE 


+ oo 


1 
B(w)B(— o) (= — pe) k(— w)K(ok(—o))K,(ok(o))o dw— 


| 
à 
aie 
vw. 
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CO 


z er a]. 
ES BC B)RK (ok*o) Kahle) io de | 
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The corresponding formula in the Fermi theory 


awY 22e? 1 
pe DO MAIRE ——_ — 2) k*(w) K,(ok*(w)) Ko(ok(w))ia d 
(47) a ae 2| (; a ge) (o)K;,(ok*(o))K(ok(w))io do , 
can be obtained from (46) by taking B — 1. 
It was shown in references (5) and (6) that in the case of à non absorbent 
medium (real &(w)): 


aw? aw Deze? _ lolo læ | o lo |e 
4 2 = = + Ya,o = K 
(#8) dæ da v2 as Cher ; | v | | v ; 
dw” Een U oe i 
49 4 - Blade , 2: 
(49) EP: Ar | (= fa 0) (a > 0) 
1-P?e<o 


where — dW%/dx is the rate of emission of Cerenkov radiation per em. For- 
mula (49) is due to FRANK and Tamm (?). The a; are the residues of &!(w) at 
the poles @. These poles are all real and simple, in the case of a real &(w): 


(50) 


The contribution of the residues in (48) gives the rate of loss per cm due to 
direct ionization and excitation. Equation (48) shows clearly that the reso- 
nances occur at the frequencies w, as it was said in the discussion of section 1. 

We may expect to get a formula analogous to (48) in the present theory. 
There are now some new circumstances due to the presence of the factor 
| B(w) |? in (46). The w are no more poles of the quantity under the integral 
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in (46), but there are poles introduced by the factor | B(w)|?, as we shall now 
prove. It was shown in references (5) and (6) that there is a pole @, of &(w) 
very close to each frequency w; (@;<«@,). We may write: 


(51) | B(w) |?k*(w)K, (ok*()) Ko(ok(w)) = B(w)B(— o)k*(w)K,(ok*())K,(ok(o)— 


K,(ok*(w))Ko(ok(@)) 7. ee 
ON NA Ce CAE 


K,(Rk(o))k 
(52) D(o) — I(Rkylo)) 4 ne 


e-Uo)L,(Rk,(œ)) . 
In the case of a real e(@), the function D is real and monotonic in the interval 
©; == ©; . Since 


(53) D(@;) TA T,(Rk(@;)) >0, Do; — 0) =-— oo, 


there is a zero ©, of D in that interval. ©, is a double pole of the quantity 
under the integral in (46). When the damping of the electronic vibrations is 
not neglected, e(w) is complex and its zeros ©; — in, and poles &; — in, have 
negative imaginary parts which tend to zero together with the damping con- 
stants g in the Lorentz formula (42). In the case of very small damping 
constants, the phase of D(w) varies approximately by —z as w varies from @; 
to @;. This shows that the zeros of D(w) corresponding to the ©, have also 
negative imaginary parts—7;. The corresponding zeros of D(—q) have posi- 
tive imaginary parts. In the limit of zero damping constants, D(w) has no 
zeros in the Cerenkov bands, because the ratio of the A functions in (52) 
cannot take imaginary values inside a Cerenkov band. (These points are 
discussed in Appendix I). 

The transition to the case of zero damping, in (46), is more complicated 
than in the case of the theory on Fermi’s lines, because the quantity under 
the integral has poles on both sides of the real w-axis. In order to obtain a 
formula analogous to (48), we shall consider equations (27) and (28) as obtained 
from the inversion of a two sided Laplace transformation, instead of a Fourier 
transformation, and the integration paths will be taken parallel to the real 
@-axis at a distance 6, in the upper part of the w-plane. The well known 
formula for the Laplace transform of a product of two functions gives | 


+09 
Ser 5 | Ey Eop Hg exp |— 216 i de = 
—co 
+cot+id 
ee | en, (= 6) K,(ok(o)) K,(ek(2id - @))k2id-—- wo dw . 
zw? &(@) 


—cot+id 
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By taking a sufficiently small 6, we may keep it constant and make the damping 
constants tend to zero, and then let 6 tend to zero. Thus we get: 


a | II . 2 
(55) lim (— = — Jim in Aa 
g—>0 dx ô—>0|9—0 TV 
+69 .16 
: 1 A re 
x i B(w)B(2i0 — w) Fr e) Ko(ok(w))K,(ok(2iö -w))k(2iö -w)ode| . 
— Go +76 


The important thing is that, in the limit, we do not get the Cauchy prinei- 
pal value definition of the integral of B(w)B(— @)[(1/e(@)) — P?]K,(ok(w)) x 
x K,(ok(— w))k(— w)w along the real w-axis, but the integral taken along the 
real axis indented by infinitesimal semi-circles centred at the ©, and lying in 
the upper part of the w-plane. We shall denote this integration path by (ol 
and drop the lim in the left hand side of (55): 


(56) en. = m | B(w)B(— 0 (= =: pe) k(— @)K,(ok(-w))K,(ok(o))w do. 


It was shown in references (5) and (6) that: 
(57) Rliok(— w)K,(ok(— w))K,(ok(w))] = arg k(w) (wo > 0). 


Therefore the only contributions to the right hand side of (56) are those of 
the Cerenkov frequencies and the residues at the @;: 


dx v? ©) 
ee) 
DA sees ~ ~ a 
ar > %,(0\8,0hl@,)Kılok(@,))Koloko;)) , 


Residue Bo) en Br ko) Æ(ohk(o))Ko(ok(o))o 
= a E(w) 
(58a) I a = 


(jk(@;)K, (ok(@;)) Ko(ok(as)) 


The flux of the Poynting vector through the cylinder at infinity gives the rate 
of emission of Cerenkoy radiation: 


AW Cer 


dx 


(59) 


II 202 [7 
er oe lee ~ B*)| Blo) to den = — 
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This formula differs from that of FRANK and TAMM (49) by the presence of 
the factor |B(w)|? which reduces considerably the intensity of the high energy 
radiation. We may write: 


dW? AW  2ie: x R 
(60) nn : a4 ae 2 35,00 \,0k(o,)K,(ok(0,))Ko(ok(o,)). 


The second term in the right hand side of (60) gives the rate of loss per em 
due to direct ionization and excitation at distances larger than o> R. 


5. — We can compute the rate of loss of energy due to distant interactions 
by using the Poynting vector of the field inside C;: 


(61) 2 18 = => [E'xH'),.ds = --z/ EE, AS = 
Oo Co 
+ co 
--8[ [EDER + EYE, + ELEY + E,Boldw. 


It is easily seen that the terms ASE and ELE, do not give any contribution 
o (61). By working out the contributions of the other two terms, we find 
again the formula (38). We shall now separate the contributions of the direct 
ionization-excitation and the Cerenkov radiation in (39), by indenting the posi- 
tive real axis at the &, corresponding to zero damping and making the damping 
constants tend to zero. Thus we get: 


22e? \ 
(62) — dW., = 224 @| K,(Rk(o)) | 1 pe) B(o)io do — 


daz Tv? Rk,(@)) \e(w) 
Cer bands 
22 ER a . 
— I, (Rk(o; ©, Residue [D(o)]! 
v2 2| 1 )) oV1 pi 2, [ )] 


It seems Dee, that the real part of the integral over the Cerenkov bands is 
the rate of emission of Cerenkov radiation per em and the contributions of 
the residues the rate of loss due to direct ionization and excitation at distances 
larger than R. It is shown in Appendix II that the first term in the right 
hand side of (62) coincides with — dW%/ds. Thus we get a more convenient 
form of the loss by ionization and excitation. 

The residue in (62) can be computed by the well known formula: 


(63) Residue [D(w)}-* = ER 


03 dw (o= @; 
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In many cases of practical interest w,R/v>1 and we can determine with 
sufficient accuracy ©, by the equation 


‘ Er: BETT an Lo(RKo(&;)) 
(64) Let): V1 — Bre(,) = — V1 —B LR (@))* 
and take 

‘3 api ff, Be \ide * 

(65) os (: 2(1 os al dé L(Rk(@)) en 


For extremely relativistic particles and w;R/v > 1 we get 


(66) Ge (ét), Fee 


20 de m ©; \ 2% 


by using the approximate expression 


à = 
Ww; — w? 


(67) e(o) = (w= w;), 


for values of w very close to &,. The contribution of the residue at w,; to 
—dW,,/dz is approximately 


Agne2e4  [w,R\-4 wok w,R  — 
38 LS > j i 5 il, ed 4/1. 221. ; 
En mov? fh | 2 | | v À v pes 


A general approximate expression of the contributions of the residues can be 
obtained from the formula (9) of Appendix III. It results from (68) that the 
contribution of the large w, is negligible, at distances larger than R. It is 
interesting to remark that their contribution is nevertheless larger than in 
the theory on Fermi’s lines. Indeed, since 


(69) Os, = — O0_G, = 27%07f;, 


it follows from (48) that the contribution of a large frequency w,, in the theory 
on Fermi’s lines, is: 


(70) 


mnie ref; oe | | (> i 
mu? 2 v v 


When w,R/v<1, we may use the approximate expressions (24), so that: 


= 1 /&,R\? OR dD E de 
71 wu —— K A -— = |- — 5 
(71), OD Al © | at v | IE € dw)w= 6) 
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It follows from (62) and (63) that the contribution of w, to —dW,,/d@ is 
approximately ; 


= Bee NE (2 2 ok _ 
(72) Te GK (0e) CRE En ’ 


and does differ little from the corresponding value in the theory on Fermi’s 
lines: 


22262 _ &,R\/de\-! = e 
(72a) IR (2) 5) (a 


v2 do/w= 3) 


Discussion of the energy loss due to distant interactions. 


6. — We shall now compare the results of the present theory and that on 
Fermi’s lines. We shall first assume that 


. R = 
(73) oe ep RTE 
on 


for all the frequencies and we shall replace the 7 and X functions by their 
asymptotic expressions 


eu 


rn 7 = a 
(74) I(u) = Va Ku) = js exp [— u] ((ul>1), 


in equations (22) and (52). Thus we get e(w,) = — 1 for B<1 and using for- 
mula (9) of Appendix III 


2826? UD inne! „ wih aw ke 
75) — —— [(1,(Rk,(@,))Rky(@;)|'@ D + f “re ’ 
(75) yp [ ı( o(@;)) o(@;)| ©; an aes = mw? es : 5 
Hence the results of the present theory do not differ practically from that on 
Fermi’s lines for non relativistic velocities and distances larger than R. 

The total loss at distances larger than R is given approximately by (44), 
for extremely relativistie partieles. It is easily seen that the formulas given 


~ 


in references (5) and (6) lead to the following result 


(76) lim 


w—>e 


aqaw® A4anze4 1154 3 tune It; 
oo > al log 1 < 


da na en me Ira," 


(77) i= 
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I, being the mean ionization potential of the atoms of the medium (atomic 
number Z). The value (76) is, in most usual cases, considerably larger than 
(44), particularly for gases. The total amount of Cerenkov radiation is, the- 
refore, much smaller in the present theory than in that on Fermi’s lines, in 
most cases. We shall now prove that this is due to the strong reduction of 
the high energy Cerenkov bands, i.e. those in which mR/v>1. |Rk(w)| is 
equal to zero at the lower limits of those bands, but has large values on a 
large part of the bands and becomes infinite at the upper limits. The term with 
the K, function in (22) is equal to zero at both the limits of a band and is 
not very important, so that: 


(78) B(w) = [To(Rko(w)) Ky(Rk(w)) RE(w) | - 


In the part of a band that gives the strong intensity in the theory of Frank 
and Tamm, we may use the asymptotic expression (74) of Æ,, and since Æ(w) 
is imaginary: 


a = 


(79) | Ba) |? = = [LL (RE (o))V R|k(o)|] * > (e(@) — 1) 


The factor |B(w)|? in formula (59) cuts down strongly the intensity when 
lo| R/v>1. 

It is interesting to remark that the factor | B(w)|? in (59) does not change 
appreciably the intensity of the low energy Cerenkov radiation. When 
|w|R/v<1 we get from (22): 


R A 
(80) B(w)=1, (els < 1; lw| A a) ; 


The preceding discussion shows that the essential difference between the 
present theory and that on Fermi’s lines comes from the drastie reduetion of 
the high energy Cerenkov radiation. Whenever the theory on Fermi’s lines 
does not lead to a large contribution of the Cerenkov radiation to the total 
loss, the results of both theories are not very different. It results from (76) 
that this may only happen when 


IR 


RI, 
(81) log 


i.e. for condensed media of low atomie number. It will be shown later that 
the loss at distances less than R is usually less than in the theory on Fermi’s 
lines. 
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The loss of energy in close collisions. 


7. — Quantum corrections become important in the computation of the ave- 
rage loss at distances less than R, as it was shown in the Introduction. It 
results from those considerations that the uncertainty principle leads to à 
broadening of the frequency ranges in which resonance occurs and the emission 
.of coherent radiation is not possible in these ranges. Since the main pola- 
rization effect at such distances consists in cutting down the direct ionization- 
excitation and replacing it by the emission of coherent radiation, it is satis- 
factory to neglect altogether the polarization effects for 0o< R and to apply 
the methods of Bethe-Bloch-Möller-Williams. At distances close to R the 
polarization effects become considerable, as it was shown by A. BoHR (®). 

In reference (#), one of us applied the Bethe-Bloch-Möller-Williams me- 
thods up to distances Rk’, R’ being close to our À, and the modified form of 
the Fermi theory to distances larger than R’. Such a treatment is satisfactory, 
when a high accuracy is not required, and it has the advantage of simplicity. 
It can be improved by computing the energy transfers at distances less than R 
with the field (E', H'), instead of the fields (E, HO) used in the quantum theory 
of Bethe-Bloch-Méller Williams. 

BLocH (1) and WILLIAMS (!?) have shown, by neglecting polarization effects, 
that, for distances less than R, >10" cm, it is necessary to use quantum 
mechanics. But this is no more necessary for distances larger than R,, if only 
the average loss is required, provided the atom is assimilated to a system of 
classical harmonic oscillators with frequencies w, corresponding to the quantum 
jumps, there being Zf, oscillators of circular frequency w;. We shall sepa- 
raté the contributions of distances less than R, and distances between R, 


and R: 
As AWes «dW, Wh, 
. da da da ’ 


—dW;,/dæ being the contribution of the atoms between the cylinders Cp, and Oz. 
We shall derive —dW_„/dx from the Bethe-Bloch quantum treatment and 
compute —dW,,/dx classically as the rate of work per em done by the field 
E' on the electrons between C,, and C,. If we neglect completely the pola- 
rization field E’, —dW,;,/dw is replaced by the rate of work done by the 
field E®. It is shown in references (5) and (6) that this rate of work 
— dW,/dæ is: 


€ 2 
dWz,  2rneet 
dz mov? 


(83) > fiko(w,;)| U(w;, R,)— U(w;, R)] ; 
j 


(84) Ulm, 0) = = 0? |B?K?(ok,(w)) — 2K 3(eky(@)) + 
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For (ae?/hv) << 1 the total loss is given by the Bethe formula 


Dm? hi 2 

ME max DCE 

wi |, (32 <a), 
10 


Ne 741 —)? 


mw? 


| =) 2runzzer | 
Bet. 5 


when polarization effects are neglected, Emax being the maximum energy trans- 
ferable to an electron by the ionizing particle. In order to obtain —dW_,,/da, 
we must subtract from (85) the contribution of the impact parameters larger 
than R;: 


aw dw / aw 2 2 
6 É Skt = Ban Ry ; = 2 1 : 
(26) da | dz i | dx i ET ho = 
daw Q7cne? 74 2mv°E ace: 
ST) == <A, __ 1 > “max 9 Be y;) UIs Ra & eZ Ih 
(27 da mw? os © 12(1 — B?) i 2 a Us, #1) ko ~ 


A more general expression can be obtained by replacing the Bethe formula 
by that of Bloch, which is valid for any value of 2,¢?/hv 


(88) 


AW.,, _2menzje‘ 8: 2mv’ 2mv’E max 
dr mv? 272 #3 


_. 2B% + y) 2e y ( ae = )}- = EH (w;)U(w,, R,)| , 


+ denoting the logarithmic derivative of the gamma function: y(u) = 
= d/dw log Iu). It is important to remark that (87) and (88) do depend 
on the f; and , individually and not only through the mean ionization po- 
tential /,. This is a very unpleasant circumstance, because these quantities 
are in general unknown. However, if we limit the energy transfers to a 
value E’ and consider the limit of very high energies, we get (*): 


da 


Zrınzie‘ my?’ RE 12.6? 
= — — 2 1+ —— ; 
me? jos 2h3 Sr oma Er 


(B’) 
(89) lim ve ee) — 
\ 


This limit value does not depend any more on the f, and w;. 


(*) This formula is valid when E' > Vme’%o, For small values of EB’, E'< 
< V’me?%o, the term — 1 inside the bracket in the right hand side is replaced by 0. 
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A rough approximate value of —dW_,/dz can be obtained by adding 
to (88) the value (83), with R replaced by R'Z 0.37R, in order to take into 
account the reduction of the loss by the polarization effects at distances close 
to R: 


(90) AWir = 2rneiet [2mv? Emax 
5 de ” mv? |T(1—Bß?) 
age ar pr 4 8 | ~ 2 filet U(w;, K')| . 


In the next section we shall derive à more accurate expression and justify the 
choice of R'. 


8. — In order to determine the density of polarization P' induced by E', 
we shall write 


+00 
(91) Fahre €'(w) exp io? SE 4] do , 
—co F 
hence 
ave 
(92) P' besoin | e(w) — 1)E'(w) exp (51) do. 
4x dn v 


—dW;,/dxv is equal to the rate of work per em done by E* on the electrons 
between (, and Cp: 


R 27 + CO 
ARS ON i oP" 
93 pe Be soa Eire 
> da ‚fe 2 | av | ( ya 
Ry 0 —co 
a +00 
| 2x | ode R | (1 — e(w)) (E(w) -E*( ))io do. 
| u 
! Rı 0 | 


It results from (32)-(33) and (35)-(36) that: 


R 
aw 
a + 2a | 2o)e de, 


Ry 


dws _ 
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de 
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195) 0%) = 5 1 6 
ef | ie(w)O(w)O(— w) {(1 — P?)I2(ok,(o)) + T( ok (w))} dw + | dolce) + 


| + O(—o)) {a — B?) Ko(oko(@))Lo(oko(@)) — Kıloko(w))Iı(o ko(c))} «0° de : 


We get from (23) and (52): 


w? dia? Ky( RB ko(w)) 
9 I do) Oo) + C er ' ae 
(96) 7 e(w)[ O(a) + C(— w)] = &(®) T (RE (0) + 
_ | km) K,(Rk(o)) , ew) k(—o) K(RKÇ—-w))| | ia a 4 
7 rem | &(— ©) DI-w) K,(Rk(— | Un ann) 
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 |Diw) Kı(Rkio))  &(-@) D(— a) Kı(Rk(—- o))| I:(Rko(o)) | 4 p Fe 


nil k(w)k(-- @) Ko( Rk(w))Ko( Rk(— w)) | m? @® 
' eo) D(»)D(— o) K,(Rk(o))K,(Rk(— o)) | 4 v? 


(1 — f)2R?12(Rho(@)) ie 


We must start with a complex e(m) and then make the damping constants g 
tend to 0. In the limit of zero damping constants, the only contributions 
to 2(o) come from the poles of the quantities under the integrals in (95), since 
these quantities are imaginary in the limit. The first term in the right hand 
side of (96) and the first term in the right hand side of (97) have the same 
poles as e(m). In the corresponding integrals the integration paths will be 
indented with infinitesimal semi circles centred at the poles of e(w) and lying 
in the upper part of the w-plane, because these poles have negative imaginary 
parts when there is damping. The second group of terms in (96) and the 
second group of terms in (97) have simple poles due either to D(w) or D(—@). 
In the corresponding integrals, the integration paths must be indented with 
infinitesimal semi-circles centred at the ©, and lying above the real w-axis 
in the case of D(w), and below that axis in the case of D(—w). The last term 
in the right hand side of (97) has double poles at the @,, in the limit of zero 
damping. The corresponding integral can be computed by the same method 
applied to (46). 
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We shall denote by Q2,(o) the part of Q(g) due to the poles of &(w) and by 
Q;(o) the part of Q2(o) due to the poles at the w,. It is easily seen that: 


2 ue ge 
(98)  Q(e) = zur 1) Eva | — B2)Ko(gh(@,))Lo(okel@,)) — 
= RR ¥ = 
— K,(0k,(@;))L,(oko(@;)) 5 N - B2)I2(oko(@,;)) + Ti(oko(@;))}| , 
— KÂER(G) 2. a MER) 
OUT Liebe) de no 


The contribution of the poles @; is small at relativistic energies, since: 


(100) lim 2,(0) — 0. 

The contributions of each pair of terms within the square brackets in (96) 
and (97) cancel, so that the only contribution to 2,(@) comes from the last 
term in (97): 


ae? 
(101) 2;(0) = vip 
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By using the well know Lommel formulas 
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The integration od 0,(0) can also be done by using (102) amd (103): 


(105) 2x | Qiao do = 2n [er do — 
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It results from the equation (18) of the Appendix II that: 
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The residue of the product of a function F,(u), regular at u, by a funetion 
F,(u) having a second order pole at wu is given by the formula: 


(107) Residue [ F(u)F(u)| = Fy(u) Residue F,(u) + 4 
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It follows from (105), (106) and (107) that: 
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The expression (66) of (dD/dw)w=3#) is not sufficient, because it is only appli- 
cable when |e(@,)|>$1. A more general formula will be derived in Ap- 
pendix III: 


(109) (2). Arene?f,; 1—B'e(&;)  .(@,) 


do)}w=ay MO (RG) 2 — Pre(@;) (1 — es)? 
The approximate expression of e(w) given in Appendix III shows that: 


B?mo; ea — £(©,))? 


110 m : 
ar wi Anneïf, 1— Belo) 


d 
me mo) 


Finally we get: 


(111) 2a | Codes" Zi €(@;) | Be Oo; y. 
j 


mw: 1— e(@,)) \2 — P?e(@; 


Ry 


x [B?{ Ii( Rk(@,;)) — 12(Rkg(@;)) } + 2(Rko(@;)) Lo Rko(@;)) Ls (Rko(@;))] x 
X [1o( Rko(@;))L,(Rko(@;)) 7? : 
This correction remains finite as v tends to c: 


Rk 


; Annee O;R e(@;) \? 
27 | 20 dd D = 21; (: 2 oe) Vee) | RE 


Rı 


(112) lim 


v—c 


For non relativistic energies the Rk,(®) have high values and the corrections 
due to due to ©, and Q, are negligible, because of the asymptotic behaviour 
of the 7 functions: 


et 1 — 4»? 
113 Ts ae ee Le 
a ARTE IR | 
We get from (111) and (113): 
n = 
(114a) 2x | 0, (0)o do = 
16n?neizes e(@;)(1 — B2e(@,) 8 5 z 
= 0 oe Se Rk j) C2 — 2Rk j 
en OB) Slag ciao pret )| O9 exp (28 (5), 
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It results from (104) that: 


R 
4a? nere4 w;R\? exp|— 2Rk,(@;)] 
Elston | Ox (o\p do AS a 
| 1(e)e de = mw? ir a; | G Rko(w;) 
Ry 


with (V1— p?~1). 


* [Oxygen 


For extremely relativistic 
energies the correction (111) is 
roughly equal to — 4rınz?e!/me?, 
as a consequence of (112). The- 


g refore we may take in (90): 
à 
& 
= = Il, 
5 Enc oe il, 
“dl = 4 6 8 10 2 4 6 6 10 2 R 


(115) 
Fig. 1. — Rate of ionization at distances larger | RES OST. 
than R, = 1078 cm. 


The saturation value of the 
total loss per em, with transfers limited to #’, is roughly: 


i 1W Dance my?2R2£" N 
(116) lim (— Se | zen ox be {y ( ae } . 
ME 


oe 
Zu Mare 
Cae 


da: 


2.2 u 
| nee 
2 2 4 6 810 2 4 6 
Fig. 2. — Rate of ionization at distances larger than R, = 1078 cm. 


This value is smaller than that given by the theory on Fermi’s lines: 


5 AW)\®  Qaneres | 4mR?E' f 12,62 
LH? = a Dey DR tee : 
yew) im | de | me oe on: a el "ke }}l 


The variation of the Bethe loss at distances larger than R, (— dW.,,/dt)Bet.» 
of — 4W,,/dæ and of — 4dW;,/dæ with the energy of the ionizing particle are 
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represented in fig. 1 for oxygen (normal conditions) and in fig. 2 for AgBr 
crystals. In fig. 2 the variation of the rate of ionisation and excitation at 
distances larger than R is also included (*). 


Diseussion of the energy loss due to close collisions. 


9. — In the computation of the energy loss at distances less than R, we 
took into account the interactions of the electrons of this region with the ionizing 
particle and with the electrons outside C;, but we neglected the interactions 
between electrons lying inside C,. The contribution of the electrons at di- 
stances less than R to the field outside C, has also been neglected in the com- 
putation of the loss at distances larger than R. If we would have taken into 
account the contributions of the electrons inside C, to the field, both within 
and outside O,, by applying classical mechanics and electrodynamics, we 
would have, of course, been led to the theory on Fermi’s lines. In our treat- 
ment, the loss at distances between R, and R consists of ionization and exci- 
tation, whereas in the theory on Fermi’s lines the main effect of the electrons 
between R, and R at high relativistic energies is, in most cases, to give rise 
to the high energy Cerenkov radiation, as we shall now prove. It was shown 
in references (5) and (6) that the rate of work done by the field on the electrons 
at distances between o and o + do, in the theory on Fermi’s lines, is 


’ 24 Fe Ke 
(118) __ dwr(o) do en o do 3 10% Ke (22) + K? (22) : 
dav MU SES v v 
with 
(119) 2a; — ap, (fi =f) 


the notations being the same as in formula (48). The corresponding quantity 
in the Bethe-Bloch-Möller-Williams theory is 
_ 4nnaïe® 


dwwz(e) 
da 


(120) — do ode > f;k2(w,) | Kr(oku(w;)) + (1—B?) K2(oke(,)) |, 
ij 


mv? 
and in ours 


dw(o) dw,(o) nr 
(121) ha de = — ran do + 27 Q2(o)o do. 


(*) In fig. 2 the value obtained with the approximate formula of Halpern and Hall 
is indicated as Fermi approximate. 
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The total loss by ionization and excitation in the interval of distances R, — R, 
in the theory on Fermi’s lines, is: 


R 


(122) IE DANONE: ue ra be cal BE 


da e mw? v 


Ry 


@;R K o,R K oR Anne, 2hw 
i Ta Ur : = me "yRl, 


At extremely relativistic energies the rate of emission of Cerenkov radiation 
per cm is given by (76) and is larger than the value (122) for most elements 
in a condensed state and even more in the gaseous state. The Cerenkov 
radiation is due to the interaction of the ionizing particle with electrons at 
distances larger than atomic dimensions and the high energy bands, which 
account for the largest contribution in that theory, arise largely from inter- 
actions with electrons at distances less than R. (This is shown by the fact 
that, in our treatment, in which the electrons at distances less than R are not 
taken into account in this connection, the high energy bands are practically 
suppressed). Therefore, we may conclude that the main effect of the electrons 
in the interval of distance R, — R, at high relativistic energies, in the theory 
on Fermi’s lines, is to give rise to the high energy Cerenkoy radiation, in 
most cases. 

In our treatment the quantity corresponding to (122) is —dW;,/dæ. It 
results from (83) and (115) that: 


: dw, Annee R 
123 a NS SE Me ln 
( ) im | Ae | = [1 R, | 


d—>C 


This value is larger than that given by (122) in all cases, excepting condensed 
media of very low atomic number. 

We have already considered in the Introduction reasons for, not taking 
into account the field of the electrons in the interval of distances R— R,, in 
the computation of the loss. We shall now examine this point in more detail. 
The influence of the field created by the electrons in the loss of energy is not 
of the same nature as that of the Coulomb force of the nuclei. It is an effect 
of the field due to the acceleration of the electrons and the coherence of the 
contributions of the different electrons plays an essential part. In order to 
see whether the interference of the fields of the electrons belonging to different 
atoms may be of importance, we must consider the times of collision. We 
(shall denote by T.(o) 


(124) Te) = = Vale BR 
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the time of collision defined, as usually, with the Lorentz transformed Coulomb 
field of the ionizing particle. The important contributions to the coherent 
fields of the electrons are due to particles whose phase differences correspond 
to time lags 7 of the or ii of T,(o). It results from the uncertainty principle 
that 


(125) AtAw >1, 


and, since we need 7 of the order of 7,(o), we must have 


(126) At<T.(o), 
hence: 

Ù 
oa, 20V1 —p? 


Important coherence effects can only occur in frequency intervals of the order 
of ©@;—@,. It was shown in references (5) and (6) that: 


Agra f; he? hi 


12 Ee E = o, 
(128) un 30); 3R? ho, 


In order to have important coherence effects the uncertainty Aw must satisfy 
the condition 


he? 
(129) Aw< RI,’ 
which requires that: 
(130) Sd a, (us 
he VT — p 13741 NEE - 2 R, 


Since the minimum of ionization corresponds to von p? = 1/3, the condition 
(130) cannot be satisfied, after the minimum, at distances less than R. The 
condition (126) is certainly too stringent, since the polarization effects lead 
to an increase of the times of collision. The preceding considerations should 
therefore not be valid when the polarization effects are important, i.e. at 
distances larger than R. Moreover the small frequencies become very important. 
In most of the interval R,— R the considerations leading to (130) are valid 
and show that the contribution of the electrons of that interval to the field 
may be neglected. The applicability of the classical treatment “of sections 
2, 3 and 4 to the atoms at distances larger than R depends essentially on the 
modification of the times of collision due to the strong polarization effects. 
We shall examine this point in the following section. 
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10. — The value (124) of the time of collision 7,(o) follows immediately 
from the formula: 


CON LES 2,e(1 — Pe : 
(131) Ey a [o?(1 — p?) au (a — vt)?]3/2 = 


We cannot apply the same method to the case of a dispersive medium, because 
we do not have any more a formula corresponding to (131). There is, how- 
ever, another method of estimation of T,(o) based on the Fourier integral (33). 
It is well known that the duration of an electromagnetic pulse, whose spectrum 
covers an interval Aw, is of the order of 2r/Aw. Equation (33) shows that 
the width Aw at the distance o is determined by the condition 


(132) ok(Aw) Z1, Nee — 


which leads to (124). In the case of the Fermi field, the function Æ,(ok,(w)) 
is replaced by K,(ok(w)): 


+00 
(133) Bu | e—'(w)k(w) K,(ok(@)) exp |iw (? = i) do . 
= He 
65) 


When is much larger than the oscillator frequencies w;, we have: 


470700? 02 
(134) e(o) = 1 Pre 1 ask?’ (> w,). 
Hence 
—— 
(135) ok(w) = vg, + o°k5(o), (o>), 


and at distances larger than R the contributions of the large frequencies to ESP 
become small, since 


(136) ok(w) > 1 (wo >o;, p > R). 


The preceding considerations show that the distance R behaves as à radius 
of action in the theory on Fermi’s lines, as it was shown in an elementary 
way by A. BOHR (*). 

It results from (33) that: 


7 v 2210. 
(137) lim BY = a jo» lo ( = ‘| do = = Ô(æ — ct). 
d =: CO C | © 


T—+>c 
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This formula shows clearly that the collision time 7, vanishes in the limit 
of » = c; the frequency range that gives important contributions to the field 
is now infinite. In the case of the Fermi field 


oo 
(138) AU | e.m)k.(w)K,(ok.(w)) exp io (=-. ) da, 
dc 7 TTC b (a 
— co 


and the transverse component of the electric field is not delta-like in the limit 
of v—6. This could also be seen by taking into account that the range of 
frequencies giving important contributions to the field does not tend to ox, 
as à consequence of (136). 

In the case of E”, (133) is replaced by (28). It results from (22) that: 


(139) B(w)= 1, (o > @,). 


The contribution of the very high frequencies is not larger than in the case 
of the Fermi field. 


Summary of the formulas. 


11. — We shall now summarize the main formulas. We shall not take into 
account the corrections corresponding to the case in which the ionizing 
particle is an electron, as we did already in section 7 (*). The total rate of 
loss per cm — dW/dx is 


aw dW, AW, 
dpe 2 a2 das” 


—dW_,/dx and —dW.,/dx being the contributions of distances smaller than R 
and larger than R, respectively, with 
(m: mass of the electron) 


2 
(141) wear 
Anne? (n: number of electrons per em®). 


The loss at distances less than R is the sum of two terms: 


(142) ie, Sue te 


=i Orsi . 
da da da (A LOUER 


(*) Our formulas are nevertheless applicable to fast electrons when high energy 
transfers may be disregarded (see MOLLER (16)). 
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— dW.,/dæ is given by a modified form of the Bethe formula: 


; dW Zrengze‘ a =a 
Cea en log 7, == de = ZL fskelo,)T(@,, R)|, 


(1434) — <1, 


(v = velocity of the ionizing particle) (z;e = charge of the ionizing particle). 
The w and f are the oscillator frequencies (circular) and percentuals, respe- 
ctively. J, is the mean ionization potential of the atoms of the ionized medium. 


The function U is: 


(144) U(@, e) = | oko( w)) — K?(ok(co))} + Ko( oko(@))K, (ok (@))| - 


The K are modified Hankel functions and 


(144a) ko(w) AR 18, 


— dW;,/dx is the contribution of the distances between A, and Rk: 


P R 
dw; 2renzze: 
(145) nn = D Iklw)| Ul@s, Rı)— Uo,, R)] + ET de. 
Ry 


The last term in the right hand side of (145) is the correction to the loss at 
distances less than R due to the polarization at distances larger than R: 


(146) Q(e) = Ale) + (0), 
(147) 2n | 2 (e)e de = — er a) Z © Kı(Bko@,)) I(Rko@,)) + 
: + K,(Rk(@;))1,(Rko@;)) — 2(B2Rky(@;))Ko(Rko(@,))L(Rko(@;)) + 
+ { Rkg(@,)Lo( Rk@,)) (RG ))] en DR Res lo). 


e(w) is the dielectric constant of the medium for waves of circular frequency w= 


ung Inne: 
(147a) @;=0;, (®,; = pole of e(w)), Res e(&) = — fr 


©; 
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The I are modified Bessel funetions 


(147) lim Q,(0) =0. 


v—>o 


The correction due to 2, is less important than that due to Q,. 


(148) 2x | Q,(0)o de = lc e(@;) JS &(@;) a € 


mw? 1 (Co Be(@; 


x [B? {Ti Rko(@;)) — B( Rky(@,))} + 2(Rkg(G;))-Lo( Bkio(@;)) Ls (Rkeo(@;)) | x 
x [1o(Rko(@;)) I, (Bko(@;)) |-? 


(@;-pole of B(w) close to w,). 


(1484) [B(o)] = R|K,(Rk(o))I,(Rk(o))kko) + 


Klo) re , 
\ e(w)ko(@) K(BK))T (Hl) ’ 
(1485) k?(w) = = — P?e(o)) , (-3< arg k(w) <3, aS 0) 


—dW_,/da# corresponds to direct ionization and excitation. The loss due to 
the emission of Cerenkov radiation is included in —dW,,/da: 


dw. dWe 22,€? Ok. |e dD 
149 ser SE I T,(Rk,( lt 
{ ) da da v? 2 ate (05 IE vip — pe oe da (= w,) ; 
AW cor 2e? / ( EU x 
149 - —. Ale ar 
(149a) a Fe | B(w) ea 8 Jo do, 
Cer 
ao>o0 
(1495) D(w) = | K,(Rk(w))Rk(w) Blo) | 
dD\-1 Arne?f, 1 — Pe; ) &(@;) 
(149¢ nerd 2 a Z mee Ma 
oe Gale no; o(Rko(o )) 2 2 — Be(@;) (1 — e(o;))? 
; UWE u 


All the preceding formulas correspond to the case of a real dielectric cons- 
tant and are not valid when there are conduction electrons. The case of a 
medium with conduction electrons presents special difficulties. Nevertheless 
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some of the formulas derived in the preceding sections are still valid in this 
case. For instance, we get from (46) an expression of — dW.,/dæ which is 
general: 


A ay, 
(151) de ea er [| Boil (aes 52 | k*(w)K,( Rk*(w)) K(Rk(w))io do. 


The expression (95) of Q(o) is also general and allows us to get a general 
R 


formula for [are do: 
Ry 


co 


22 02 RR? , 
1° a | ee | O(@) |2{ T(Rko(o) —I2(Rk(w)) + 


2v2c? 


A —f)R 


a 


(152) 2a | 2(o)o de = 


Ry 


+ 2(P?Rk,(@))Io(Rk(o))I(Rk,(w)) ho do — 


-| e(w)(Olw) + C(— w)){ Ky(Rko(w))Lo(Rk(@)) + K,(Rko(w))1,(Rky(@)) — 


= 2(BRk(o)) 1 Ky(Rko(w)) 1, (Rko(w)) 12% + ; 


(152a) Cee Br 


in (EC) e '(w)B(w)K,( Rk(w)) — K(k, (ao DIET L(Rk(@)) F1 


Diseussion of some experimental results. 


12. - We shall examine only some of the existing experimental results in 
this paper, a more systematic discussion will be published later. We shall 
not consider here the experimental results regarding only the effects of the 
polarization on the stopping power, since our main interest is in the 
relative contributions of ionization-exeitation and Cerenkov radiation to the 
stopping power. 

The most important experimental results, from the point of view we are 
assuming, are those of GosH, JoNES and WILSON (1) for mesons in oxygen. 
In this case the increase of the ionization after the minimum is considerable, 
and the difficulty of explaining it with the theory on Fermi’s lines is obvious. 
Indeed, all the relativistic increase in that theory is due to the emission of 
Cerenkov radiation and the absorption of the radiation within about 2 mm 
is required to explain the experimental value of the increase. The present 
theory leads to an increase of ionization of about 30%, by neglecting the 
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ionization due to the absorption of Cerenkov radiation. (The loss at distances 
larger than R is almost entirely due to Cerenkov radiation, as shown by the 
discussion of section 5). This value agrees satisfactorily with the experimental 
one. It is important to remark that the theoretical value of the total increase 
is obtained by comparing the loss at the minimum with the loss at saturation, 
and is not affected by the uncertainty on the values of the oscillator fre- 
quencies w, and oscillator percentuals f,, since in the region of the minimum 
of ionization the rate of loss 

in a gas depends only on the 


iz mean ionization potential 


and the value at saturation 
does not depend on the a; 
and f;. The mean ioni- 
BO EBEN ee [Es zation potential was taken 
an as 89.2 eV for oxygen, i.e. 
the experimental value ob- 
tained by LIVINGSTON and 


Nb, of ion-pairs per cm 


45 = 
} E/nge? BETHE (12) for air multi- 

ohh 2 ATEN 9100 2 al NOONAN, Pe eS To plied by 87.22. 
Fig. 3. — Total rate of ionization (transfers less The variation of the rate 
than 1000 eV). of ionization for oxygen at 


normal conditions is shown 
in fig. 3, which corresponds to the values of the v, and f,, given in table I. 
There is some uncertainty in these values, but the effect on the rate of loss 
is not considerable. The agreement with the experimental data is satisfact- 
ory, by taking as 34.5 eV the average energy necessary to create an ion pair. 
The transfers in elose collisions were cut at 1000 eV. 

In the case of silver and bromine, the results were found not to be very 
sensitive to the choice of the w, and f;. The two sets of values A and B 
of table I led to results differing always by less than 0.5%. The mean ioni- 
zation potentials of Ag and Br were taken as 470 eV and 367.5 eV respectively, 
by interpolating the data of MATHER and SEGRE (1°) and taking into account 
those of BAKKER and SEGRE (2°). The energy transfers in close collisions were cut 
at 5000 eV (see JANSSENS and HUYBRECHTS (?1)). The results are insensitive 
to the choice of this cut-off value. By neglecting the absorption of Cerenkov 
radiation, the increase in ionization from the minimum to the plateau in AgBr 


(18) M. S. Livineston and H. BETHE: Rev. Mod. Phys., 9, 246 (1937). 
(9) R. MATHER and E. SEGRÉ: Phys. Rev., 84, 191 (1951). 
(2°) C. J. BAKKER and E. SEGRE: Phys. Rev., 81, 489 (1951). 
( 


es 
21) P. JANSSENS and M. HuUYBRECHTS: Bull. Cent. Phys. Nucl. Brux., n. 27 (1951). 
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crystals was found to be 3.8%. The total rate of loss, including Cerenkov 
radiation, with transfers in elose collisions eut at 5000 eV is of 5.96 MeV/em 
at saturation, the corresponding value in the theory on Fermi’s lines being 
6.42 MeV/em. The value 3.8% seems too low to explain the observed increase 
of grain density of about 7%. Moreover the saturation of the grain density 
occurs faster than that of the total rate of ionization. The discrepancy bet- 
ween the theoretical results and experimental data can be considerably atte- 


Silver Bromide 


} = A: +4 

0.95 Experimental values _| _ 
of Morrish 

0.90 


4 2 4 6 8 10 2 4 6 8 10 


Fig. 4. — Total rate of ionization (transfers less than 5000 eV). Curve I: Contribu- 

tions of all the shells. Curve II: Contribution of the outer shell (j = 1) neglected. 

Curve III: Contributions of the outer shell and one third of that of the second shell 
(j = 2) neglected. 


nuated by taking into account that small energy transfers are probably not 
sufficient to develop a grain of AgBr. By neglecting the contribution ‘of the 
electrons of the outer shell, the increase from minimum to saturation be- 
comes 4.8%, and an increase of 6% is obtained by neglecting entirely the 
contribution of the outer shell and 33% of the contribution of the next 
shell; this procedure leads also to a faster saturation. The theoretical results 
and the experimental values of MORRISH (*) are compared in fig. 4 (*). The 
agreement seems satisfactory, when the small energy transfers are neglec- 
ted, although the saturation of the grain density is still faster than 
that of the reduced theoretical loss. The eventual existence of an instru- 
mental eut-off of the very small transfers in the photographie emulsion can 


(2) A. H. Morrısm: Phil. Mag., 43, 533 (1952). 
(*) The normalizations of the three curves in fig. 4 are of course different, in order 
to get the same saturation plateau. 
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TABLE I. 
Silver Bromide Oxygen 
7 ee à a = : = pr = 
à Shells 837, tj ue) = Shells | 8f, v; 
1 N Br 8.4 1.025 1.502 L2p 4.2 54.31 
2 N Ag | 20.4 3.33 shes L 2s 2.0 108.62 
3 M Br 18.0 | 7.18 6.00 K 1.8 727.10 
4 M Ag 18.0 10.62 10.62 
5 IR AB | 7.2 33.50 | 33.50 | 
6 L Ag 7.2 72.45 72.45 
m ike Br | eA: 277.80 277.80 
8 RACE 1.4 526.60 526.60 
a 11.76 - 1015 8-1 | 0.1861-1015 s-1 
| R | 4.065- 1077 em 2.502 -10752cm 
| Vi Ag 470 eV Br 367.5 eV 89.2 eV 
0; 
vy; = ‘ 
27a 


be put in evidence by accurate measurements in the region before the 
minimum of ionization, where the polarization effects are very small. Such 
an effect would make the effective mean ionization potential for the develop- 
ing of grains in the photographic emulsion higher than the true mean ioni- 
zation potential. The effective mean ionization potential for the developing 
of grains may eventually depend on the average size of the grains and other 
conditions affecting their sensitivity, as well as on the processing of the plate. 

Unpublished results of Voyvonıc (*) for ~-mesons and electrons in the 
photographic plate show also an increase of about 6-7% and the agreement 
with the present theory seems quite satisfactory by assuming the existence 
of the instrumental effect discussed in connection with the MORRISH data. 
The results of L. JAUNEAU and F. HUG-BousSER (2?) show the existence of 
an increase of the grain density after the minimum in the photographie emul- 
sion, but the scatter of the results in different series of measurements appears. 
to be considerable and a precise evaluation of the increase after the minimum 
is difficult. 

The results for the photographic emulsion can probably be improved by 


(*) Private communication. à 
(3) L. JAUNEAU and F. HuG-Bousser: Journ. de Phys. (in press). 
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using a more accurate expression of the dielectric constant, in which the effect 
of a shell is not assimilated to that of a single oscillators (see WICK (?)). 


Note added in proof. — Dr. A. H. Morrisx has called our attention to the paper 
of Mc Drarmip (Phys. Rev., 84, 851 (1951)), whose results agree rather better with 
ours with respect to the beginning of the saturation in photographic emulsions. 


APPENDIX I. 


Let us denote by F(w) the quantity within the square brackets in (40): 


(1) F(o) 


a 2e(w) K,(Rk.(o)) LME 

Rk.(o) K(Rk.(@)) 
It was shown in references (5) and (6) that the w corresponding to real values 
or poles of the function e(w) defined by (42) have negative imaginary parts, 
when the damping constants are not 0. (qm) has therefore no other singu- 
larities in the quadrant of the positive real and positive imaginary axis 
besides that at the point m — 0. In order to sbow that #{(w) has no zeros. 
in this region, we shall consider the closed contour J’ formed by an infinite- 
simal quarter of circle centred at the origin, the positive real and positive 
imaginary axes and a quarter of circle centred at the origin with infinite radius. 
F(w) is real and positive on the positive imaginary axis and on the infinite 
quarter of circle, and the variation of its argument along the infinitesimal 
quarter of circle and the positive real axis is equal to zero. Since the variation 
of the argument of F(w) along the closed contour J’ is equal to zero and F(w) 
has no poles, either on the contour or inside it, when the damping constants. 
in (42) are not 0, it follows from the argument principle that it has no zeros 
in that region. 

F(o) is related to the function D(w) defined by (52): 


: . 2e(o) KılRk.(o)) 
(2) F(o) + Rk.(o) K,(Rk.(o)) Ze | 


: Rk,(o) K,(Rk, 
eh lem ean er Ta K > a @ 


We have shown that D(w) has zeros in the neighbourhoods of the w,, whose 
imaginary parts are negative when the damping constants are finite. It results 
from (2) and (3) that, in this case, F(w) has also zeros in the neighbourhoods 
of the w;, with negative imaginary parts. 

We shall now prove that, in the case of a real e(w), the function D(w) has 
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no zeros in the Cerenkov bands. It follows from (52) that, at a zero @ of D(w): 


K(RK@)) (RO) hol) ~ 
) 


Æ(RKG) ae 


(4) ) I,(Rkg(@)) k(0) 


If @ would lie in a Cerenkov band, we would have: 


(6) 


Since J (u) and J,(u) do not vanish simultaneously for any real value of u 
and the same happens with NV,(u) and N,(w), equation (6) will only be satisfied 
when k(@) — 0. It is therefore not compatible with (4). 


AMPIPEINIDILX HE 


It results from (22) that: 


1 
(1) (= _— pi )z (co) Kq(Rk(o)) (= = 6 B(—w)K,(Rk(— )) = 
1 1 
— B(w)B(-w) 2 ©) (= ) K,(Rk(@)) — B-(w) > = ») KulRk-0) = 
= RI, (Rk,(®))B(®)B(— w\ (Gs — B ET K,(Rk(w))K,ı(Rk(— w)k(— w) — 
= = —#) K,(Rk(— OR RK) 7 
e(— @) 
hence 
2) & Es =p ) Blo) Ko( Rk(w))ic ) = 
= RI,(Rk,(w))| B(o) |? R Fr x ul (Rk(o))K,(Rk— w))k(— w)iw| . 


This equation shows that (39) is equivalent to (46) for o = R. 
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In the Cerenkov bands arg k(wm) = — 7/2, therefore: 
(3) k(—- wm) = exp [ix]k(w) , (w > 0, 1 — Bew) < 0). 
By taking into account the well known formulas 
(4) K,(exp [ixlu) = K,(u) — ixli(u), K, (exp [ir]u) = — K,(u) —ixl,(u), 


and (21), we can simplify (2) in the Cerenkov bands: 


| 1 | 
5) & (a Se | B(o)K,(Rk(o))ia Le 


= 50 (=. e) L(Rk,(©))| B(w)|?, (© > 0, 1 — B?e(w) < 0). 


It results from (5) that the first term in the right hand side of (62) is equal to 
the right hand side of (59). 

When the damping constants are finite, we can obtain regular analytic 
functions in a strip around the real positive w-axis starting from the functions 
k(m) and B(w) defined on that axis by (16) and (22), with 


(6) kl) = —V1— 6, 


k(— w) and B(—«), as defined by (16) and (22), are not obtained by analy- 
tical continuation of the regular analytic functions k(w) and B(w). We shall 
introduce the analytic functions k(w) and B(w) defined by the equations: 


(7) k(w) = k*(w),  B(w) = B*(w), (m > 0). 


It results from the well known symmetry principle of Riemann and Schwarz (*) 
that, with a suitable choice of the branches: 


(8) ko) = (k(o*))*,  B(@) = (B(&*))* 


By introducing the function G(w) 


(9) Go) =w (<5 = ) B(w)K,(Rk(o))[l — RI,(Rk,(0))B(®)K,(Rk(o))k(o)] : 


equation (1) can be written as: 


(10) FmGo)—=0: for. -w>0. 


(*) See, for instance, Z. NEHARI: Conformal Mapping (New York, 1952) p. 184, 
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Hence we have 

(11) Hw*) = (Hw))* , 

as a consequence of the symmetry principle. It results from (11) that: 


(12) Res G(w) = — {Res G( w) }*. 


wj+in; wj-in; 
This equation can be obtained by comparing the values of Î G(o) do and 


| G(o)dw, I” being a circle with sufficiently small radius centred at the point 
rx 

©, + in, and J* the circle with the same radius centred at w,—in;, both 
circles being described in the same sense. By introducing the expression (9) 
of G into (12), we get: 


(13) _Res Fe — pt) B(w) B(w)1o(Rko(o)) (Hb) (Eko) Ho) | + 
oj +ing 
1 N 
+{ Re le +) BoB I (Eko) Ko( Bk()) K( RE(c0)) (oo Ie 
Ra. (a o)K „(Kk(o))o | 


In the limit of zero damping constants, equation (13) becomes: 


(14) Res 


d; 


(in +) #4 ©) ER OA (Ko) (0) RH = 


a & = — p? )x o(RK(@;))@; Res B(«) . 


&(@; 8; 
It is easily seen that for any analytic symmetric function H(k(a), k(w); ©) 
(15) H(k(o), k(o); ©) = H(k(w), k(®); ®) , 


regular in neighbourhood of the poles ©; + in, of G(w) and taking real values 
for positive w, we have: 


(16) Res H(k(o), k(o); o)B(o)K,(Rk(w)) (= ~ po . 


win; 


= Res L(w) + Res L(o), 


Oj+ ini @j—inj 
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x B(®)B(w)T,(Rk,(w))K,(Rk(w))K,(Rk(w))Rk(o)o £ 


In the limit of zero damping constants, we get: 


(18) H(k(@;), k(@;); ©) K(Rk(&;) ler >), Res B(w) = 
= Res H(k(o), k(w); ¢ MM — ß? ) B*(w)Io( Rko(w)) Ko( Rk(w)).K, (Rk(w)) Rk(@)o 
dj 


This equation shows how it is possible to transform the expression of the loss 
at distances larger than R by direct ionization-excitation obtained from the 
equation (58) of section 4 into that given by (62) of section 5, which involves 
only the residue of [D(w)]". 


APPENDIX ITI. 


We shall now prove formula (109) of section 8. In order to get an appro- 
ximate value of dD/dw at the zero ©,, we shall start from the formula (52) 
of section 4 and we shall neglect the variation of the functions J, as well as the 
variation of the ratio of the two Æ functions. We may write: 


(1) D(w) = I,(Rky(o)) + en — = ee T,(Rko(@)) . 
a= é(@) 


With the aforementioned approximations we get: 


| aD. Bulk) LCA) Ar nn 
(2) Te Leo) V1 — P2e(w)] . 
Since 
: r _V1i— peo) 2— Brew) _, de 
(3) x Le 1(o)V1 — pre oi] SES Tee az TETE (0) Te 
and 

si WEEE) vx A BHGe)) (REG) 
(4) . e'(@,)V1— Bea) = — VIP 7 KORG) T,(Bko(@,)) ’ 
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dw “a 


(= ®,) 


(5) ae Ca 


1,7 m KılRk@,)) Lo(Rko(@;)) 2 — Area) ~ (dE 
NE RE) LE) œ al! le 


It results from (2) and (5) that: 


(6) 


aD Le 2 — pre(B,) 
ona) a2 De 


= de A 
AR I,(Rko(®;)) a (es, Ee@;). 


It results from formula (50) of section 4, combined with formula (119) of 
section 9, that in the neighbourhood of @,; we have: 


‘ 47 2-/2f. 
(7) ae 
Oo Os 
We get from (7) 
de d œ(e(æ) — 1)? 
= RTS RE ue, SE aa I = 
(8) me Nas ser: Zu 27?0?f;e(w) 
hence 
(9) ER Nr dovne*f; 1 — Be(@s) _&(@s) a 
do a mo;To(Rko(® ;)) 2 — P?e(@;,) (e(@;) = a): + 


This is formula (109) of section 8. 

It is interesting to remark that formula (66) of section 5 can be obtained 
immediately from (9). This is important because (66), for extremely rela- 
tivistic energies, gives a good approximation when @, differs little from on 
whereas (7) is only satisfactory when « differs little from @,. Formula (9) 
gives therefore a good approximation when &, is close to either ©; or @;, i.e. for 


an entire interval which contains &,, at least for extremely relativistic energies. 


RIASSUNTO (*) 


Si formula una teoria della perdita d’energia di una particella carica, in cui il mec- 
canismo della perdita & differente da quello della teoria di Fermi. La nostra teoria 
conduce a una distribuzione della perdita tra l’eceitazione per ionizzazione e l’emis- 


(*) Traduzione a cura della Redazione. 
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sione di radiazione Cerenkov che differisce notevolmente da quelle delle teorie del tipo 
di quella di Fermi, per quanto il potere di frenamento non sia essenzialmente diverso, 
essendo soltanto nella nostra teoria alquanto minore che in quelle. Questo trattamento 
conduce a un aumento della ionizzazione diretta per parametri d’urto maggiori delle di- 
mensioni atomiche dopo il minimo relativistico, l’aumento essendo considerevole nel 
caso dei gas. L’effetto di saturazione di Fermi della perdita a distanze superiori alle 
dimensioni atomiche esiste anche nella nostra teoria, per quanto in virtü di un diffe- 
rente meccanismo. L’aumento relativistico della’ perdita nelle interazioni a distanza è 
dovuto largamente a un aumento dei raggi d’azione per la ionizzazione e l’eccitazione 
(meceanismo di Bohr-Williams), ma & dovuto anche a un contributo della radiazione 
di Cerenkov (meccanismo di Fermi). La saturazione € dovuta a una limitazione del- 
l’'aumento dei raggi d’azione dovuta alla polarizzazione del mezzo e alla saturazione 
dell’emissione di radiazione Cerenkov. La presente teoria conduce a una modificazione 
della formula di Frank e Tamm per il tasso d’emissione della radiazione Cerenkov, 
l’emissione di radiazione Cerenkov di elevata energia essendo considerevolmente ridotta. 


IL NUOVO CIMENTO Vor. IX, N. 9 1 Settembre 1952 


Sulla corrispondenza tra elettrodinamica classica e quantistica. 


G. MORPURGO 


Istituto di Fisica dell’ Universita, Centro di Studio per la Fisica Nucleare del O.N.R. - Roma 


(ricevuto il 2 Luglio 1952) 


Riassunto. — Si mostra come nel caso non relativistico e di dipolo 
esiste una completa corrispondenza tra l’equazione di Dirac-Eliezer 
classica dell’elettrone puntiforme e l’equazione che in elettrodinamica 
quantistica regola la variazione col tempo del valor medio quantistico 
in uno stato (expectation-value) della posizione dell’elettrone stesso. Nel 
caso di un oscillatore armonico la corrispondenza & cosı stretta che il 
valor medio quantistico della posizione dell’oscillatore soddisfa proprio 
all’equazione di Dirac-Eliezer. Come conseguenza di questo si ha che, 
per lo meno nel caso dell’oscillatore armonico, Vintroduzione della co- 
stante di Planck non toglie le difficoltd delle soluzioni fisiche e non fisiche 
dell’equazione di Dirac-Eliezer. Ed inoltre, nella variazione col tempo 
del sopramenzionato valor medio, si ha, sia pure in un caso particolare, 
un esempio di risoluzione di un problema di elettrodinamica quantistica 
completamente indipendente dalla teoria delle perturbazioni. 


Introduzione. 


Nel presente lavoro vogliamo mostrare come, almeno nel caso non rela- 
tivistico (n.r.) e nell’approssimazione di dipolo (d), esista una interessante cor- 
rispondenza tra l’equazione di DIRAC-ELIEZER (">%*) (D-E) classica dell’elet- 
trone puntiforme e l’equazione la quale in elettrodinamica quantistica regola 
la variazione col tempo del valor medio in uno stato (expectation-value) della 
posizione dell’elettrone stesso. Sempre nelle approssimazioni sopra accennate 


(7?) P. A. M. Dirac: Proc. Roy. Soc., A 167, 148 (1938). 
(2) J. ELIEZER: Rev. Mod. Phys., 19, 147 (1947). 

(3) M. Markov: J. Phys. USSR, 10, 159 (1946). 

(4) D. Boum e M. WEINSTEIN: Phys. Rev., 74, 1789 (1948). 
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e nel caso particolare di un elettrone soggetto a una forza elastica di richiamo 
(oscillatore armonico) la corrispondenza & cosi stretta che il valor medio quanti- 
stico della posizione dell’oseillatore soddisfa proprio all’equizione di D-E. 

Per stabilire la ora detta corrispondenza, nel $ 1, nelle approssimazioni 
sopra riferite, dedurremo V’equazione (non relativistica) di D-E dal formalismo 
hamiltoniano; nel $ 2 mostreremo come la stessa deduzione possa essere svolta, 
riferendosi ai valori medi, nel caso quantistico e specializzeremo le formule 
al caso accennato in cui la corrispondenza & più stretta, cioè al caso dell’oscil- 
latore armonico; nel $ 3, per mostrare che l’hamiltoniana usata nelle appros- 
simazioni n.r. e d. & abbastanza generale, ed oltre a cid per completare un 
nostro precedente lavoro (5), dedurremo brevissimamente da essa, sempre per 
l’oscillatore armonico, il Lamb-shift. Seguirà nel $ 4 una breve discussione. 
Aleuni ealcoli saranno svolti in appendice. 


1. - L’hamiltoniana dalla quale partiremo (nella quale gli oscillatori longi- 
tudinali e scalari sono stati eliminati) sara dunque la seguente: 


2 


e e 
= 2 V ey an 7” — Nv ee OU 
on © + V(x) mc 2 (ex P)dx sis AMC? fe (er'ex)gidu 4 


K, ke 


GIE 


1 
+5 SCORE + M). 
AR 


In (1) le notazioni sono quelle stesse del WENTZEL (°). L’indice e in alto & stato 
inserito per ricordare il fatto che in questo $ le variabili che compaiono in (1) 
sono variabili dinamiche elassiche (numeri e). Nel prossimo $, in cui faremo 
uso della stessa hamiltoniana essa sarà invece affissa di un indice q à deno- 
tare che le grandezze che in essa compaiono sono degli operatori (numeri q). 
V(x) in (1) & l’energia potenziale di un campo esterno di forze, m, la massa 
meccanica dell’elettrone; l’indiee di polarizzazione & sempre sottinteso. 

Notiamo che l’hamiltoniana (1) contiene (nel quarto termine) la diffusione 
della radiazione (formula di Thompson). 

Le equazioni di Hamilton associate all’hamiltoniana (1) sono: 


1 e 1 
= BEN 

PT me Vi A 

: dV 

PR esse 
(2) ce 

de = Pr; 

: Gal Cael 
für = k?qx + moe ye (p'°ex) == Moe? 17 > (er ex) 


(5) G. MORPURGO: Nuovo Cimento, 7, 313 (1950). 
(6) G. WENTZEL: Quantum Theory of fields. (Interscience Publishers, 1949). 
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dalle quali si ricava: 


% DE CRE 5 

(3) MoX = a . yır > Cndn » 
De ec = 

(4) Qe + c7kq, = AT (rene 


Tl secondo termine a secondo membro della (3) à semplicemente il campo elet- 
trico nella posizione dell’elettrone: si tenga presente che è: 


À, Am D eng exp [ik'x]. 
R 


1 
yılz 


(5) Gut) = s - Al (e,°%(s)) sin ke(t — s)ds + 


Pig twig 
+ ro COS ke(t — ty) + xo + sin ke(t — to), 


dove do © Geo rappresentano i valori iniziali (al tempo { —t,) di q;,(t) e della 
sua derivata. 

La questione di quali debbono essere tali valori iniziali è probabilmente in 
relazione con il problema delle soluzioni fisiche e non fisiche dell’equazione 
di D-E. Noi non vogliamo qui discutere questa questione e porremo percid, 
come si fa (sottintendendo la cosa) di solito, quando l’equazione di D-E viene 
dedotta ammettendo che il campo nella posizione dell’elettrone sia quello 
dovuto ai soli potenziali ritardati: 


(6) Iro = xo 0€ 


Ci limiteremo inoltre a discutere soluzioni della-(3) in cui x(t,)= x, = 0. Se 
si impone la (6) & ragionevole imporre questa condizione. 

Sostituendo la (5) (in cui sia stato fatto qu = do — 0) in (3), con suc- 
cessive integrazioni per parti e facendo uso dell’integrale di Dirichlet si giunge 
rapidamente (v. appendice I) alla equazione di D-E n.r.: 


= 
= 
R: 
u 
| 
| 


dove si € posto 
e? 
M Mg Ar Mem + Mem dk ’ 


ossia si € rinormalizzata la massa. 
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Come abbiamo detto non vogliamo qui discutere la (7); essa è gia stata 
oggetto di numerose discussioni; il problema delle sue soluzioni fisiche e non 
fisiche resta ancora un problema aperto. L’unica osservazione (7) che vogliamo 
fare è che, se si impongono le condizioni iniziali sul campo, la condizione ini- 
ziale sull’accelerazione non & imponibile arbitrariamente ma & determinata 
dalla (3). Se come abbiamo supposto à Qro = Yo = 0 risulta 


(8) ee ee 


Si noti che nella condizione iniziale (8), che come abbiamo detto deve essere 
imposta se le condizioni iniziali sul campo sono le (6), compare la sola massa: 
meccanica m, (che non ha significato fisico) e non la massa m = mo + Men 
che compare nell’equazione di D-E e che é la sola ad aver significato fisico. 
Cid fa pensare che le (6) non siano le corrette condizioni iniziali da imporsi 
al campo. 


2. — Se x, P, 9x, Pr Vengono in (1) interpretati come operatori soddisfa- 
centi alle usuali regole di commutazione, (1) diventa l’hamiltoniana dell’elettro- 
dinamica quantistica e come tale sara indicata con «H. La generica funzione 
d’onda del nostro sistema sara indicata con W(x, g,,t) (essendo sottinteso 
che & funzione di tutte le q,) ed il valor medio all’istante ¢ di un generico 
operatore O nello stato W(x, q,, t) sara indicato con < W*(t)|O|Y(t)> o anche 
pit brevemente con O(t) o infine anche con O semplicemente: 

öl) = | P(x, Qu OW, qu, t) xT] du - 
J k 
B allora ovvio che le equazioni che regolano la variazione col tempo dei valori 
medi degli operatori x, P, dx, Pr, hanno esattamente la stessa forma delle 
equazioni (2) che regolano la variazione col tempo delle variabili classiche 
x P, Qe» Pr. Ossia hanno la forma: 


a 10e e 1 Su 
Fm me vg Ale 
Ran 
((2') dx 
dh CD 
x e le er | = 
DZ u. ya) nn er ee 


x 


(?) Questa osservazione & dovuta a MARKOF (?); indipendentemente, il prof. B. 
FERRETTI mi ha fatto rilevare la cosa durante una discussione per la quale gli sono 
grato. 
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La cosa è un’ovvia conseguenza del fatto che: 


dyes eee 
= — = a 
(9) 3, 0 =; FAO CE), 


e che 


Procedendo a partire dalle (6’) in maniera completamente identica a quella 
del § 1 e imponendo le condizioni 


eo = Gio = 0, 
si ricava l’equazione: 
A SV Bt 
10 mx lern: 
= On: (se) Grec? ” 


La (10) € una equazione che lega la derivata seconda e la derivata terza di x 
con il valor medio di dV/dx nello stato che si è preso in eonsiderazione (che é 
del tutto generico). Ora il valor medio di dV/dx in generale nen pud essere 
espresso mediante x. Tuttavia un caso importante & quello in cui si abbia 
V = 1} ox? (oscillatore armonico di costante di elasticità o). In tal caso 


dV/dx — ox e quindi la (10) si riscrive: 


e? 


Er ea 
6e? 


(11) mx + ox — 


Questa equazione & una equazione differenziale in x(t) ed & perfettamente iden- 
tica all’equazione in x(t) che si ha (per un oscillatore armonico) classicamente. 
La sua soluzione generale è ben nota ed anzi in questo caso notoriamente esiste 
anche la soluzione fisica (8). 

Notiamo che, se si ammette che l’hamiltoniana 27 (in cui si ponga !/:0x? 
al posto di V(x)) abbia (degli 0) uno stato stazionario, in tale stato la (11) é 
identicamente soddisfatta, nel senso che tutti i suoi termini sono identicamente 
nulli per ogni valore di ¢ (°). 


(8) Vedi, ad esempio, A. LOINGER: Nuovo Cimento, 6, 361 (1949). 

(°) La cosa puö anche essere controllata direttamente osservando che x, x sono 
evidentemente nulli e che x & pure nullo; quest’ultima eircostanza dipende dal fatto 
che 7H & invariante rispetto allo scambio di x in — x e q; in — q,, e conseguentemente 
le autofunzioni appartenenti a stati stazionari sono o pari o dispari rispetto allo scambio 
simultaneo di x, gx in — x, —q,- Percid il valor medio di x, che & una funzione dispari, 


x 


in uno stato stazionario & nullo. 
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Nel caso generale la (11) permette di determinare esattamente, cioé a tutti 
gli ordini in e, il moto del baricentro di un pacchetto di probabilità dell’elet- 
trone, qualora per lo meno, cosa che puö essere dubbia, ci si fidi della solu- 
zione fisica. Discuteremo più a fondo la (11) nel $ 4. 


3. — In questo numero, per mostrare che l’hamiltoniana che abbiamo usato 
nella deduzione delle (7) (10), sebbene abbia validità limitata, è tuttavia tale 
che molte delle caratteristiche dell’hamiltoniana esatta sono in essa rispec- 
chiate, dedurremo, facendo uso di essa, il Lamb-shift di un oscillatore armonico. 
La cosa à immediata. La trasformazione canonica di funzione generatrice 8, 
tale che (1): 


(12) Sie ee [| | = parentesi di Poisson 
dove 
1 2 1 n2 1 x »2m2 ! le? 2 
(13) Ay = amt VE D DE 3 2 (cp + koi), 
e 1 
(14) BF === > (er P)Qx , 


me VIE À 

elimina l’interazione del primo ordine H, dando luogo a un’interazione sup- 
plementare del secondo ordine pari a: !/2[S,, H,]. Il Lamb-shift dell’oseilla- 
tore armonico in uno stato y, stazionario imperturbato con zero fotoni vale: 


il 
(15) <Oyr |5 [Sı, Hall wn 0 > rinormalizzato. 


Se si sceglie $, pari a: 


| me V2 £ (o/m,) — k?e? 


(16) A e 1 zart: AR C(p ex) Px 


la (12) e soddisfatta. 
Con la (16) risulta: 


. i GF tl erahnen! c:(e.'p)? 
[S,, H] = > u Sean 


1 == TASSE = ae SES > u . 
ae 2 mee? V po k?e2—(a/m) 2 mee v2 k2c?2 — (o/m,) 


bol rR 


Il secondo termine a secondo membro della (17) è identico a quello che ab- 
biamo ricavato in un precedente lavoro. Il primo termine naturalmente non 
compariva in tale lavoro in quanto classicamente esso è nullo quando ci si 


(1°) Nelle (12), (13), (14) le variabili possono. essere indifferente mente numeri € 
o numeri q. 
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metta in condizioni in cui l’energia del campo elettromagnetico sia nulla come 
ci si era messi allora; infatti se l’energia del campo elettromagnetico & nulla 
tutti i g, Sono nulli e quindi tale primo termine scompare. 

Ma quantisticamente non possiamo evidentemente annullare i qg,. Sussi- 
stono sempre le fluttuazioni del campo elettromagnetico che fanno si che il 
primo termine sia in ogni caso non nullo. Ora, come è stato visto nel lavoro 
citato il secondo termine a secondo membro della (17) vale semplicemente 
(m,„/2ms)p? ed & quindi nient’altro che il termine di rinormalizzazione della 
massa del secondo ordine. B facile controllare che il primo termine & proprio 


il responsabile del Lamb-shift. Basta porre in esso 4% =y® (a, + dy) € 


fare il valor medio in uno stato di data energia imperturbata e con zero fotoni. 
Il calcolo e immediato (14) e lo accenniamo in appendice II. Quello che ei 
interessa & di aveı mostrato come la hamiltoniana “H da noi considerata con- 
tenga il Lamb-shift. Ed inoltre ei interessa di aver confermato il risultato 
del nostro precedente lavoro (5); in tale lavoro dal fatto che se il campo di 
radiazione viene trattato classicamente non si ha il Lamb-shift deducemmo 
indirettamente che tale Lamb-shift doveva essere una conseguenza delle flut- 
tuazioni del campo elettromagnetico in accordo con il lavoro di WELTON (?2). 
Qui, mediante la trasformazione canonica (12), abbiamo mostrato diretta- 
mente come il Lamb-shift sia una conseguenza delle fluttuazioni. Come ab- 
biamo gia detto il primo termine a secondo membro della (17), responsabile 
del Lamb-shift, non compariva nel nostro calcolo in quanto, essendo il ealcolo 
classico avevamo potuto porre i gq, = 0 fin dal principio. 


4. — Conclusione. 


Facendo uso di una hamiltoniana che presenta certe restrizioni ma che 
tuttavia nel caso non relativistico contiene le principali caratteristiche del- 
Vhamiltoniana corretta, abbiamo mostrato come sussista tra le derivate se- 
conda e terza del valore medio di x e il valor medio di 3V/dx una relazione 
che & la corrispondente dell’equazione di Dirac-Eliezer. Nel caso particolare 
dell’oscillatore armonico il valor medio di x obbedisce proprio alla (ossia varia 
col tempo proprio secondo la) equazione di D-E. Quindi possiamo enunciare 
il risultato: la soluzione esatta (cioè indipendente dall’uso della teoria delle 
perturbazioni) del problema di elettrodinamica quantistica : determinare 


(1) T. A. Werron: Phys. Rev., 74, 1157 (1948). 

(2) Si tagli naturalmente un a un certo k,,,,. Si tenga anche presente che m, nel 
secondo membro della (17) puö essere sostituito con m a meno di termini di ordine 
superiore in e?. 
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la variazione col tempo del valor medio della posizione di un oscillatore armo- 
nico in interazione col proprio campo, si riconduce alla risoluzione esatta (cioè 
senza far uso della teoria delle perturbazioni) dell’equazione di Dirac-Eliezer. 

Quindi, per lo meno nel caso di un oscillatore armonico, l’idea (2?) che la 
introduzione della costante di Planck possa togliere le difficoltà dell’elettro- 
dinamica classica non sembra corretta. 

Una indagine approfondita del problema delle soluzioni fisiche e non fisiche 
dell’equazione di D-E diventa assolutamente necessaria. Nel caso dell’oscilla- 
tore armonico esiste, come & ben noto, una soluzione fisica dell’equazione di D-E. 
Questa soluzione fisica & probabilmente la soluzione alla quale tende 1a solu- 
zione perturbativa. Comunque, per quanto sappiamo, ci sembra il presente 
il primo caso nel quale, sia pure partendo da una hamiltoniana di validità 
limitata, si conosce la soluzione esatta (cioè indipendente della teoria delle 
perturbazioni) di un problema di elettrodinamica quantistica. 

Si noti a questo proposito che l’hamiltoniana di Bloch-Nordsiek & assai 
più limitata di quella di eui abbiamo fatto uso in quanto non contiene il ter- 
mine in A? e non consente quindi il fenomeno della diffusione. 


APPENDICE I. 


Dalla (5) integrando per parti si ha: 


t 
a 


(e.'x(s)) cos ke(s — t)ds = 


t 
L z 

== DE (e.'x(s)) sin ke(s — t) I a] (ex: x(s)) sin ke(s — 1) ds = 
lt, 


ty 


t e il \t 


e , : REN aS 
us (ex x(s)) sin ke(s —t) RT a; (ex: &(s)) cos ke(s — t) À 
: t 
— DE | (ex: x(s)) cos ke(s —t) ds, 


to 
conseguentemente l’espressione — (e/e)(1 /V¥?) N e.g, risulta somma di tre ter- 


mini; il primo si annulla se poniamo x(t)) = ary il secondo vale: 


il e? 


2 x 1 2 ; 
7 2 ade e,(e,°X(t)) = M.m&(t) E > — Om Ja e|, 
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dato che il contributo dell’estremo inferiore è nullo come si vede facilmente; 
il terzo é: 


t 
1 e? i wee 
(Zu + y > js O | (e,-x(s)) cos ke(s —t)ds = 
ty 
co t t 
oe ei ies e? : SIDE Cet 
= I») dk | x(s) cos ke (s —t)ds =—, lim | x(s) ve vale == 
37°C? | 370°C? ko 09 | c(s —t) 
0 to to 
9 \ 
25 hi of sin kyu aan Ce 
~ 8268 > coude \e Y= Gne ” 
—c(t—t,) 


Onde segue la (7). Si noti che se (¢—t,) + co questo ultimo passaggio € per- 
messo solo se x come funzione di ¢ tende a zero in modo conveniente per 


t + co. Questa osservazione pud essere utile nella discussione del problema. 
delle soluzioni fisiche. 


APPENDICE II. 


Il valor medio del primo termine della (17) in uno stato in cui l’oscillatore 
armonico ha energia determinata e non si hanno fotoni vale: 


Ckimax 
LÉ A o he We BENT, Ne Ck 
3 AB = a — = — — = dk= 
(5) 2 m°e? V 2 k2c? — o/m 2k 127? a | ee — (o/m) 
0 


: de Pas 
— 192 \re)\me) ° "Avolm” 


come si vede AB & eguale per tutti gli stati dell’oscillatore armonico. Lo stesso 
risultato si ha naturalmente dalla formula (5) del lavoro di BETHE (1). 
Indicato con n il numero quantico di un generico stato dell’oscillatore armo- 


4 — 
nico, facendo uso delle notazioni di Bethe e in più posto b— a/omh?, tenuto 


(3) H. A. BETHE: Phys. Rev., 72, 339 (1947). 
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presente un fattore 47 dipendente dalle diverse unità, la formula (5) di 
Bethe dà: 


RCK ax 


En es, & [rer 3 > (nel — En) _ 


7H Seer a: 


Ackmax 


Er hrs [arc | — b?nh Voim bin +1 1% Vo/m |. 
= Voim Im + fick h Vom + hek 


~ 12n®heim: 
0 


ckmax 


h ck dk ; 
i of |) = aie + termine che — 0 al tendere di 
me c?k® — (o/m) ke At 

0 max @ ° 


La (4), coincide con la (3),. 


SUMMARY 


We show that in the non-relativistic case and in dipole approximation a complete 
correspondence exists between the classical Dirac-Eliezer equation for a point electron 
and the equation which, in quantum electrodynamics governs the variation with time 
of the expectation value of the position of the electron. For an harmonic oscillator 
the correspondence is so close that the expectation value of the position of the oscil- 
lator satisfies just the Dirac-Eliezer equation. As a consequence of this we have the 
result that, at least for an harmonic oscillator, the introduction of Planck’s constant 
cannot solve the difficulties of the physical and non physical solutions of the Dirac- 
Eliezer equation. Moreover, in the variation with time of the mentioned expectation 
value, we have (naturally under the mentioned limitations) an example of resolution 
of a problem of quantum electrodynamies completely independent from perturbation- 


theory. 
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Interferometria ultrasonora nei gas (III). 
Velocita e assorbimento di ultrasuoni nell’anidride solforosa. 


S. PETRALIA 


Istituto di Fisica del? Universita - Bologna 


(ricevuto il 4 Luglio 1952) 


Riassunto. — Si sono misurati la velocita di propagazione e l’assorbi- 
mento di ultrasuoni, di varie frequenze, nell’anidride solforosa. Vengono 
calcolati il calore specifico molecolare ©, dovuto alle oscillazioni interne 
delle molecole, il tempo di rilassamento + dell’energia di oscillazione e 
la probabilità P,, di diseceitazione degli stati di oscillazione. Risulta 
0, = 15476..cal/mole, 22 1, S5.1 0222876077, = 0.041052. 


1, — Mentre esistono numerose ricerche sulla dispersione della velocità di 
propagazione e sull’assorbimento del suono, di origine molecolare, nell’anidride 
carbonica, poche e incomplete ricerche sono state fatte, in tal senso, sull’ani- 
dride solforosa. Per questa sostanza sono stati determinati aleuni puuti della 
curva di dispersione da parte di W. T. RICHARDS e J. A. REDD (!) e dei punti 
della curva di assorbimento da parte di E. GROSSMANN (?) e più recentemente 
da E. F. Fricke (°). Si tratta in ogni caso di determinazioni fatte a frequenze 
sonore basse. Col metodo cosidetto delle impurità, ideato da KNUDSEN, FRICKE 
prevede un massimo di assorbimento per l’anidride solforosa a una frequenza 
ultrasonora di 1040 kHz, per cui la ragione di dispersione e d’assorbimento 
molecolare per questa sostanza deve estendersi verso le frequenze ultrasonore 
elevate. 

Pud quindi essere interessante completare la curva di dispersione della 
velocita per valori della frequenza ultrasonora, o del rapporto f/p tra la fre- 


(1) W. T. RrcxarDs e J. A. Rew: Journ. Chem. Phys., 2, 193 (1934). 
(2) E. Grossmann: Ann. d. Phys., 15, 681 (1932). 
(3) E. F. Fricke: Journ. Acoust. Soc. Amer., 12, 245 (1940). 
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quenza f e la pressione p del gas, tali da raggiungere il valore limite della velo- 
cita V_, per frequenza infinita, e allo stesso tempo confermare l’andamento 
teorico dalla eurva di assorbimento costruita da FRICKE. Ho eseguito pertanto 
misure di velocità e di assorbimento su questa sostanza, servendomi di un 
interferometro ultrasonoro, già da me (4) descritto. 

Come sorgenti ultrasonore ho usato dei cristalli di quarzo oscillanti sulle 
frequenze di 166 e 1400 kHz. Il gas veniva preparato per azione di acido sol- 
forico su mercurio a tem- 
peratura elevata, veniva 
quindi condotto su ani- 
dride fosforica, poi lique- 


v3 (ns 


x Schweikert 


fatto e distillato più volte. 50000; “Richards e Reid ( 
+Aurore 


Dal refrigerante veniva im- 
messo nell’apparecchio a 
pressione di una atmosfera. 
Nelle misure la pressione 
del gas nell’interferometro 
poteva essere ridotta con- 
venientemente, onde otte- 
nere diversi valori del rap- 
porto f/p. La temperatura 
de] gas veniva tenuta co- : 
stante, a 20 °C. 470001. 2 

Aleune delle numerose 
determinazioni eseguite di 
velocità di propagazione 
sono raccolte nella tabella I, mentre nella fig. 1 è segnata la curva V(log f/p). 
Per quanto si sia curata la preparazione del gas, onde ottenerlo puro in ogni 
esperimento, pure da una serie di misure all’altra si ha una notevole variazione 
nei risultati, ei valori di V? determinati con la sorgente di frequenza più bassa 
sembrano stare su una curva diversa da quella che potrebbe riunire i punti 
che danno V2, ottenuti con la sorgente di frequenza mageiore. Nel grafico sono 
anche segnati dei punti dedotti dalle misure di RICHARDS e REID e due punti 
dedotti dalle determinazioni fatte da G. SCHWEIKERT () alle frequenze di 
5920 e 13062 Hz. Quest’ultimo autore da per la velocita del suono nell’ani- 
dride solforosa a 0 °C il valore 209,4 m/s. Da tale valore si & ottenuta la velo- 
cità a 20 °C, tenendo conto che la velocita di propagazione del suono in un 
gas reale pud essere espressa, in funzione della temperatura, da un’equazione 


49000; 


480007 


Ig, VA 


(4) S. Perratia: Nuovo Cimento, 9, 351 (1952). 
(6) G. Scuwexerr: Ann. d. Phys., 48, 593 (1915). 


54 - Il Nuovo Cimento. 
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del tipo 


RTq Rd 
2- — feos, 
(1) Vie (1+), 


dove g e d sono fattori di correzione dipendenti dalla temperatura e dalle co- 
stanti di van der Waals, a e b. Tale equazione é quella proposta dalla scuola 
di W. H. PIELEMEIER (%). Trascurando l’influenza della temperatura sul ter- 
mine di correzione d, risulta per la velocità V,, a bassa frequenza e per la 


jn 


O15- 


» Fricke 
o Grossmann 
4 Autore 


0,05 


10 100 1000 4000 F (kHz) 
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temperatura di 20°C, V, = 217,3 m/s, mentre sembra possa assumersi, se- 
condo le mie misure, come velocità limite per frequenze elevate e per la tem- 
peratura di 20 °C il valore V = 224,3 m/s. 

Il coefficiente di assorbimento dell’energia ultrasonora per lunghezza 
d’onda, u, è stato determinato seguendo il procedimento indicato da PIELE- 
MEIER (7). I valori ottenuti mostrano scarti notevoli da una serie di misure 
all’altra; si è ritenuto di conservare, come abbastanza probativi i risultati di 
misure fatte su un campione di gas, per f/p uguale a 1400 e 2560 kHz/atm. 
E iso = 0,132 € 560 = 0,094. Tali valori si inquadrano abbastanza bene 
sulla curva di Fricke, riprodotta in fig. 2, insieme con le determinazioni di 
questo Autore e quelle di GROSSMANN. 

Secondo un computo recente di J. J. MARKHAM, R. T. BEYER e R. B. 


(6) W. H. PIELEMEIER, H. L. Saxton e D. TEILFAIR: Journ. Chem. Phys. 8 
106 (1940). 
(7) W. H. PIELEMEIER: Phys. Rev., 34, 1184 (1929); 41, 833 (1932). 
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LINDSAY (8), il coefficiente d’assorbimento classico dell’ampiezza sonora nel- 
Vanidride solforosa a 20 °C e a pressione atmosferica vale « — 1,37 :10-1%f2 em! 
e quindi per il coefficiente d’assorbimento classico dell’energia sonora per lun- 
ghezza d’onda, y, Si ricava pe = 204 = 2,74fV :107%, 


“ns log f/p V (m/s) Poe log f/p V (m/s) 
ite 4 = Wr Be eat EN eed ee 

166 | 5,213 DES 1400 | 6,143 220,3 
15299 217.7 | 6,324 221,5 

| 5.398 218,7 6,470 221,1 

5428 219,1 6,547 RR 

| 5,665 220,0 6,595 221,9 
5,710 220,6 6,700 223,4 

5,747 220,7 | 6,801 223,9 

| | ee re! 
| | 224,3 


Per f=1400 kHz risulta 1,4906 0,008 e per f=2 560 kHz risulta 455001 = 0,015, 
valori questi melto inferiori a quelli determinati sperimentalmente; l’assorbi- 
mento & percid qui prevalentemente di origine molecolare. 


2. — Come & ben noto, l’esistenza di una banda di dispersione ultrasonora, 
che per l’anidride solforosa risulta sia daile misure di velocità che dalle misure 
di assorbimento, & conseguenza della lentezza con cui avvengono gli scambi 
di energia tra i gradi di libertà di traslazione e i gradi di libertà interni di 
oscillazione delle molecole del gas. Secondo le recenti tabelle di A. EUCKEN 
e K. H. HELLWEGE (?) le oscillazioni fondamentali delle molecole dell’ani- 
dride solforosa hanno i numeri d’onda 1151,2 em; 1361 cm71; 520 em. 
Tl corrispondente calore specifico a volume costante per grammimolecola alla 
temperatura di 20 °C, nell’ipotesi che tutte e tre le oscillazioni siano eccitate 
a questa temperatura, à date da 


» pS p (2855:8\ 9 (19908) (7479) | 
2) GE Nott |, ET EE 


ove F(x) & la funzione di Einstein. Il contributo più importante a C; è certa- 
mente dato dal terzo termine della (2), Si trova 0, = 1,525 cal/mole. Ora C; 
non influisce più sul valore del calore specifico per grammimolecola della so- 


(8) J. J. MARKHAM, R. T. BEYER e R. B. Linpsay: Rev. Mod. Phys., 28, 353 (1951). 
(9) A. Eucken e K. H. HELLwEGE: Atom. und Molekularphysik, 2, Molekeln I. 
(Bertin, 1951), p. 335. 
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stanza, quando questo viene determinato da misure di velocità sonora fatte 
a frequenze opportunamente elevate. Introducendo nella (1) i valori della 
velocità di propagazione V, e V_, si ottiene per i calori specifici molecolari 
a volume costante C, e C, 


Co = 7,611 cal/mole, C = 6,135 cal/mole, 


e perciö C; = OC, —C, =1,476 cal/mole, che è prossimo al valore determi- 
nato per via spettroscopica. FRICKE, basandosi sui dati spettroscopiei di BADGER 
e BONNER per l’anidride solforosa, calcola per la temperatura di 23 °C, 
Co = 7,524 cal/mole, C, = 6,108 cal/mole e ©; — 1,416 cal/mole; viene fatta 
l’ipotesi che solo le due oseillazioni fondamentali più basse della molecola SO, 
siano attive à quella temperatura. 

Per il calcolo del tempo di rilassamento 7 caratteristico dei processi di ecei- 
tazione e diseccitazione per urto degli stati di oscillazione delle molecole, la 
teoria fornisce varie espressioni, per vero dedotte nell’ipotesi che il gas si possa 
considerare perfetto. Si ha 


1 Cs 1 | (Cy. + RC, 
+ RO,’ 


= 7 Sn(fipho Ca 2alp)m (6, = 


dove (f/p), rappresenta l’ascissa per la quale la curva V?= V?{f/p) presenta 
un punto di flesso, mentre (f/p), rappresenta l’ascissa per la quale la curva 
che da il coefficiente d’assorbimento molecolare per lunghezza d’onda, 
u = u(fjp), ha un massimo. Le frequenze ultrasonore corrispondenti sono 
legate dalla relazione 


y, 
(4) Im = Vv 1 . 


La determinazione di (f/p), nel grafico di fig. 1 soffre naturalmente del- 
l’incertezza con la quale tale grafico & stato segnato. Si pud assumere 
(f/p),, =1070 kHz/atm; dalla (4) risulta allora per la frequenza di massimo 
assorbimento (f/p)m = 1036 kHz/atm, che praticamente coincide con la fre- 
quenza prevista da FRICKE, e cioé 1040 kHz. Secondo la prima espressione 
di (3) il tempo di rilassamento vale 7 —1,84:10-7s, mentre per la seconda 
espressione di (3) alla frequenza di massimo assorbimento & r = 1,86 107 s. 
FRICKE dä per r il valore 1,81:10-7s. Pertanto l’accordo tra le determinazioni 
di FRICKR e quelle da me fatte pud dirsi soddisfacente. 

Secondo la teoria della dispersione ultrasonora, il tempo di rilassamento 7 


(20) P. Vicourrux: Ultrasonics (London, 1950), p. 40. 
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pu essere. ancora calcolato,.per ogni frequenza ultrasonora, dal coefficiente 
di assorbimento per lunghezza d’onda, determinato per quella frequenza. E 


(5) Umol = Usper — Mel = ug + Q?n? ; 
con 
ILE = 
bie Ü 
Q? ==, n= 2nft = 
1+ R Co 
er 


Il coefficiente d’assorbimento per lunghezza d’onda molecolare w,,,, Si ottiene 
manifestamente dal valore di per determinato sperimentalmente dal quale si 
sottragga il valore di w,, classicamente calcolato. Per la frequenza di 1400 kHz 
dalla (5) si ottiene += 2,610” s; alla frequenza di 138 kHz, coi dati di 
FRICKE, quali si possono ricavare dalla curva di fig. 2, risulta 7 = 2,3 108. 

Il coefficiente di massimo assorbimento molecolare 


(6) ft 


am 
mol 


vale 0,152, mentre secondo FRICKE € Ummoı = 9,148. 

Possiamo ancora calcolare il numero di collisioni Z,, che devono aver luogo 
perché un quanto di oscillazione si tramuti in energia di traslazione moleco- 
lare, e la probabilità P,, che una molecola perda per urto la sua energia di 

’ 10 Ï 1 Le) 
oscillazione. Secondo la definizione di r, à : 

0 
Ent 


dove K,ı € Ky. stanno a indicare, rispettivamente, il numero di volte che nel- 
Punitä di tempo si eccita un quanto di energia di oscillazione a spese del- 
Venergia di traslazione o viceversa il quanto di oscillazione si tramuta in 
energia di traslazione, e la temperatura 9 & definita dalla relazione RO = hey, 
con k costante di Boltzmann, À costante di Planck, ce velocita della luce e 
y numero d’onda. 

Dalle (7) si ottiene 


7 


ji KGa VE 
RT 


? = 5 = 
Kur Ko Kyo 


(7) C= 


| — exp 


(à) | ET (: + exp 


4) 


Per’ la frequenza delV’oscillazione fondamentale di deformazione della mole- 
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cola di SO,, di numero d’onda 520 em-!, si ha allora K,, = 0,5-107. Ora à 


(9) T0 = 5 = 


e qui N esprime il numero di urti che una molecola subisce in media in un 
secondo. Prendendo per diametro delle molecole il valore dato da FRICKE, 
cioè 3,3:10-® em, risulta a 20°C e a pressione normale, N= 9,17-10° urti 
per cm? e per secondo. Corrispondentemente abbiamo 


Las ass 0510 


L’Autore ringrazia il prof. G. VALLE per l’interesse con cui ha seguito questa 
ricerca. 


SUMMARY (*) 


The propagation velocity and the absorption of supersonic waves of different 
frequencies in sulphur dioxide have been measured. The molecular heat capacity ©, 
due to the internal vibrations of the molecules, the relaxation time + of the vibration 
energy and the probability P,, of disexcitation of the vibration states are calculated. 
The A. finds ©; = 1.476 cal/mole, 7 = 1,85-10-7s and P,, = 0.54-10°. 


(*) Traduzione a cura della Redazione. 


IL NUOVO CIMENTO Vor. IX, N. 9 1 Settembre 1952 


Sulla probabilità di localizzazione degli elettroni 
negli sciami estesi. 


G. PoIANI e C. VILLI 


Istituto di Fisica dell’ Universita - Trieste 


(ricevutc il 10 Luglio 1952) 


Riassunto. — Viene studiata la distribuzione degli sciami estesi nell’aria 
applicando la teoria della cascata della componente elettrofotonica svi- 
luppata dal gruppo di Padova. Si trova che la probabilitä che un elet- 
trone venga localizzato nell’anello x’, x'+ da’ è D(x') dx’ = 0,153(dx'/x"8), 
in cui x’ & misurato in unità Molière. Questo risultato viene confrontato 
con quelli previsti da altre teorie. 


1.- E noto che MoLrèRe (!) ha ricavato una espressione del tipo 
D(x') da’ = 0,1(da'/ax’?) per la probabilit che un elettrone venga rivelato alla 
distanza (x', dx’) dall’asse dello sciame, essendo x’ misurato in unita ©’=(E./P)@. 
Sulla base di questa relazione, valida per distanze relativamente grandi, & 
possibile prevedere, in disereto accordo con l’esperienza, il numero di coinci- 
denze che si ottengono con contatori posti a distanze variabili tra venti e 
centocinquanta metri, mentre per distanze superiori, aggirantisi attorno ai 
trecento metri, il numero di coincidenze previsto & inferiore a quello sperimen- 
talmente registrato. 

Questo fatto denota che il calcolo del coefficiente numerico della anzidetta 
probabilitä differenziale & stato fondato su di una schematizzazione fisica del 
problema che limita lo sviluppo laterale della cascata e pud essere, almeno par- 
zialmente, imputato all’approssimazione di Arley, usata da MoLIERE, la quale, 
trascurando i processi radiativi sotto l’energia critica, conduce ad un numero 
troppo piccolo di elettroni di bassa energia, che sono principalmente respon- 
sabili dell’allargamento dello sciame. 


(1) G. Morrère: Cosmic Radiation (Heisemberg, 1946), p. 26. 
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Un calcolo più dettagliato à stato successivamente eseguito da EYGEs (?) 
supponendo che la distribuzione radiale, mediata su tutto lo sviluppo longitu- 
dinale dello sciame, costituisca una accettabile approssimazione della distri- 
buzione radiale in ogni punto, tranne che in prossimita dell’origine e della 
coda. Il calcolo, fondato sulle espressioni della traccia elettronica e della lar- 
ghezza media quadratica date da BELENKY (?) e valide anche al di sotto del- 
Venergia critica, conduce per il coefficiente di Moliere al valore 0,162. Questo 
risultato deve, perd, essere accolto con qualche riserva, perchè il procedimento 
adottato da EYGES tiene conto della variabilita delle sezioni d’urto relative ai 
processi Compton e radiativi non in modo sistematico, bensi come correzione 
sui risultati numerici ottenuti usando la sezione d’urto asintotica, il che com- 
porta una valutazione troppo approssimata degli effetti, opposti sui risultati 
finali del calcolo, di detti processi nella zona delle energie molto basse. 

Scopo di questo lavoro é il calcolo del coefficiente di Moliére applicando la 
teoria della cascata elettrofotonica sviluppata dal gruppo di Padova (458), la 
quale consente di impiegare, nel modo pit semplice e diretto, gli usuali proce- 
dimenti matematici fondati sulla trasformazione di Laplace. 

Il calcolo verra eseguito sulla base delle due seguenti ipotesi: 

a) che le perdite di energia al di sopra di € —= 161 MeV (limite inferiore 
di validità della sezione d’urto asintotica relativa al processo di bremsstrahlung) 
siano dovute esclusivamente ad irraggiamento con sezione d’urto costante e 
al di sotto esclusivamente alla ionizzazione; 

b) che il processo di produzione delle coppie si arresti all’energia 
M — 631 MeV (limite inferiore di validità della sezione d’urto asintotica rela- 
tiva al processo di formazione di coppie). 

Queste due ipotesi sono accettabili sin tanto che si considera l’aria come 
sede dello sciame a causa della elevata perdita per ionizzazione (8 = 98 MeV), 
di poco inferiore al valore limite e’. In particolare, l’ipotesi b) equivale ad am- 
mettere che alla quota variabile ¢’ i primari elettronici 0 posseggono una energia 
sufficiente per produrre particelle rivelabili alla quota t, ed allora vengono 
diffusi entro ’angolo massimo 6,,,,— (66 10%/E), oppure, in caso contrario, ven- 
gono assorbiti per ionizzazione nelle regioni periferiche dello sciame senza 
contribuire in modo rilevante al processo moltiplicativo della cascata nella 
successiva unità di radiazione. 

E possibile in tale modo ricopdurre il calcolo del coefficiente di Molière 
a quello dello sparpagliamento medio quadratico, mediato sullo sviluppo longi- 


tudinale dello sciame, dei //(E,,0,t) elettroni che supporremo rappresentati 


(?) L. EygEs: Phys. Rev., 74, 1801 (1948). 

(8) S. BELENKY: Journ. Phys. USSR, 8, 305 (1944). 

(4) N. DALLAPORTA e E. CLEMENTEL: Nuovo Cimento, 3, 235 (1946). 
(°) 
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dello spettro integrale di CLEMENTEL) (*) e generati da (Evo, H, t) primari di 
energia H >», (spettro elettronico differenziale di ROSSI-GREISEN) (7). 

Tenendo conto dell’ipotesi 6), si pud prevedere che il valore numerico del 
coefficiente della probabilita differenziale, calcolato seconde questi presupposti, 
risulterà approssimato per difetto rispetto al valore vero. 


2. — La probabilita che un elettrone di energia H venga diffuso entro l’an- 
golo solido 27 sin dd nel percorrere la distanza dt! espressa in unità di radia- 
zione, è, tenendo conto dell’ipotesi b) (sin @~ 0; 1 — cos 0 ~ 62/4) 


E.\? . 46 
dé = «(Z) (5) tae 


in cui «(Z) = [2-In(183-Z-"/%)}-! e H, — 21 MeV & Venergia caratteristica di 
scattering. Seguendo il procedimento di NORDHEIM (8), si trova che il numero 
totale di elettroni di energia maggiore di zero, generati da primari di energia 
E>n e localizzati alla distanza (x, dx) dall’asse dello sciame é: 


Le ‘QE! f we 
(1) N(E,, 0) = «(ZyR® fa | ae | MB, E’, t')I7(B', 0, t— t'\(¢—t')2 dt". 


Siccome il numero totale di essi che possiede energia Æ > 0 é espresso dalla 
traccia integrale: 


F(E,, 0) =| II(E,, 0, t) at, 


0 


—~ 
bo 
— 


risulta che la probabilità differenziale che un elettrone venga localizzato nel- 
Vanello circolare (x, dx) 6: 


N(E,, 0) de 


“ HA HE = M 5 
(3) D(a) da F(E,, 0) = 
in eui 
a CON Æj t 
a(Z) EB "dE Ty, 1 1 / 

== == at) == I 19, AKI, On i SEO, 

(4) M rn | PE [x x II I ) 
0 No 0 


(6) E. CLEMENTEL e L. FABBRICHESI: Nuovo Cimento, 5, 78 (1947). 
(7) B. Rosst e K. GREISEN: Rev. Mod. Phys., 13, 240 (1941). 
(8) J. NisHIMURA e K. Kamara: Progr. of Theor. Phys., 6, 628 (1951). 
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Tenendo conto che lo sparpagliamento quadratico medio ¢é: 


(5) I (By, 0, t) R20, > = aff Es, E', DIT", 0, 1—t')(t—t')? di, 


No 0 


dalla (4) si ricava, dopo aver integrato la (5) sulla quota, 
(6) M= az) <R?(0)) , 


in cui < R20) 0)> & lo sparpagliamento medio quadratico di tutti gli elettroni di 
energia H >0, mediato sullo sviluppo longitudinale dello sciame. 


3. — Applicando il teorema della Faltung, la Z-trasformata della (5) risulta 
co 


(7) [exp [— At]< R2(0, t)) 17 (BE, 0, t) dt = 


ele latin ea ee 
= 2 “u un, 3, |, 


"o 


in cui 2 e À sono, rispettivamente, il sımbolo operatorio e il parametro di 
trasformazione, L, e Ij sono le L-trasformate dello spettro differenziale di 
Rossi-Greisen di alta energia e di quello integrale di Clementel. Per 2 = 0 (s=1) 
il primo membro della (7) si identifica con il primo membro della (5) integrato 
sulla quota e, quindi, essendo il prodotto dei due operatori 2 e 2+ uguale 
all’operatore unita, la (6) diviene: 


&(2) Ze 5 
FES | La(Ep, E, 0) E 


Ne 


(8) M — 


d° ’In(E', 0, 2] dh’ 
d? 1=0 BE? 


Per il calcolo della (8) & necessario considerare il fatto che, quando gli elet- 
troni che alimentano il processo moltiplicativo alle quote inferiori hanno alla 
quota intermedia ¢’ un’energia che & uguale o prossima al limite inferiore 7, 
Vapprossimazione implicita nell’aver sviluppato la teoria verticale della ca- 
seata in base es Be differenziali incompleti di ee. «A» 
[L,(Eo, E, 4) =— (H(s)E/A (JE) e L(E, W, À) = L(s)Æ/s-À;(s): Ws*4)], 
diviene molto soie perchè le loro espressioni nn 


Go + À Bol as. 
LA (BE, E, A) = ul ans i 


1 Fi C(s) E,\s 
LH Ww Wee eae Laan REN 
fa RR am Re Ces n 
6-00 
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ede dei contributi non solo dei poli sull’asse reale positivo, proporzio- 
nali a (E,/E)’, ma anche di quelli esistenti nel Semipiano negativo, dato che, 
a causa del processo moltiplicativo rappresentato dal secondo membro della (at), 
la condizione E,>E, che rende trascurabile l’effetto dei secondi rispetto a 
quello dovuto ai primi, non & più verificata. 

Di questa circostanza analitica, strettamente connessa al fatto fisico che 
col progredire del processo moltiplicativo verso le basse energie e le quote 
inferiori, gli elettroni generatori coesistono in bande energetiche che sono 
sempre più contigue à quelle degli elettroni generati, si pud tener conto risol- 
vendo le espressioni della teoria verticale (GE 


E/no 


[A + A(s)]L_(Bo, B, 2) = 2 | 1,( Bs, 


Wile 
AoW) + XW) (Re, W, 2) = | Tal Be 


0 


Vax 1 
1e 1) Dv) , 
® ® 


dopo aver sostituito le (9) in luogo delle L, e L, incomplete. Conseguente- 
mente si trova che le Z-trasformate dei quattro termini dello spettro di 
Clementel sono: 


RER sun 
Zu (9,04) ile ek RL Ok | jan 
d-i00 42 | Wks 
Dn ey OA) a eo Als ni Erz et sane) R 
34100 4 
a wy ae IPC 
Ae 
ADI, i = EIER one uae 
ao) P(E) 
IB, 0, De ol] ar hi aon rial a 


= iCO 


E 
ds 
oe he „jean? ) ; 


ly too 
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in cui nella quarta delle precedenti relazioni la funzione q(0, A, s) tiene conte 
dell’effetto dovuto alla ionizzazione. 

Il calcolo dell’integrale che compare nella (8) si scinde in quattro integra- 
zioni distinte del tipo: 


+ico co 
H(1) E,£2 0? 1 RD AE) | dE' 
11 TER EE ~ dr | ——— BZ D). 
( ) k À 460) Iri d 2 N at A(s) | Vp i E'st3=r 10 ( ) 
— 109 No 


in cui si è posto per L,(B,, H, À) l’espressione approssimata di Rossi-Greisen 
dato che il termine 1/Æ'? contribuisce a rendere trascurabile l’effetto dei poli 
del semipiano negativo. La funzione ®,(A,7) è una qualsiasi delle funzioni 
M, N, P, Q, che compaiono nelle (10) e y si identifica, a seconda dei casi, con 
7» oppure €’. Siccome quando E,>E le sezioni d’urto assumono il loro valore 
asintotico costante, à da ritenere che porre Z,— oo nell’estremo superiore di 
integrazione non alteri in modo rilevante il calcolo; con questa posizione @ 
possibile risolvere l’integrale in campo complesso applicando il teorema dei 
residui nel polo s +2 —r—0; fissato dall’integrazione sull’energia. 

Esprimendo la (3) in unità w= (E,/ß)«, la (11) diviene nei quattro casi 
richiesti dalle (10) 


ARE: Ze B\? > [M(A, s +2) 
al nod, vole 24? Fe, + A(s) FE 
I B \? 0? [N(4,s +2) 
FE Sagal al | d4?| À + Als) ke 
(12) 
Bl Se an 
er en e'] 2042| A + Ads) Jao’ 
_ AZ) ByE (1) (B\ 2 (QA, 8 +2) 
Da te SE dA2| A + Als) [3-07 


in cui per l’aria 1/&(7,22) = 9,06. La (6) diviene: 
(13) M= F1 0 > jis); 


ove la traccia elettronica integrale & nel nostro caso: 


7 Ev C(1) 209 p(0,0,1) — g(0,0,1) 
14) P(B,, ( EOP oe iy a à 
CR En D An USE er 2’ | 


Dalle relazioni (12) e (14) risulta che il coefficiente di Moliére & indipen- 
dente dall’energia del primario elettronico. 
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4, — Per calcolare le espressioni I, (k = 1, 2, 3, 4) e la traccia integrale ab- 
biamo applicato alle equazioni di diffusione pubblicate in (% 5) un procedimento 
risolutivo approssimato consistente nell’eguagliare le Z-trasformate degli spettri 
elettronico e fotonico di bassa energia [L-(H), E/(1—v), t)/L-(B,, E,t) e, 
rispettivamente, L_(E,, W/v, t)/L-(E,, W,t)] ai corrispondenti rapporti delle 
L-trasformate di alta energia. In tale modo il calcolo numerico si semplifica 
notevolmente perchè il parametro J che nei due citati lavori era stato definito 
in funzione della quota (individuata dalla condizione di sella relativa agli 
spettri non corretti per effetto Compton) e del rapporto fra l’energia del pri- 
mario e del secondario, pud essere identificato con il parametro s caratteri- 
stico della teoria degli sciami. 

Nelle menzionate equazioni di diffusione € stato inoltre introdotto il ter- 
mine 


No 1 
jrs E du 
Im, W,t)k(W, W— E)dw= | (m. _ ) Mi, 


i ) 
B IDE 

che esprime il numero, di elettroni originariamente di energia inferiore ad Æ 
1 quali, per collisione con un fotone di energia inclusa nell’intervallo (W, W+ dW) 
eutrano nell’intervallo energetico (Z, E + d#). Questi elettroni trascurabili nella 
schematizzazione unidimensionale della cascata, essendo per lo più lenti, con- 
tribuiscono fortemente allo sviluppo laterale dello sciame e sono responsabili 
in modo piuttosto rilevante delle coincidenze registrate con contatori posti à 
distanze superiori a 150 metri. 

L’introduzione degli elettroni Compton modifica, come si pud facilmente 
rilevare, la funzione &(H, À, n) che compare nella teoria verticale di DALLA- 
PORTA, nel senso che nella funzione integranda alla sezione d’urto y(W) rela- 
tiva al processo di produzione di coppie si aggiunge la sezione d’urto Compton 


x(W) = | ko) a. 


ujaw 


L’ influenza sullo spettro differenziale elettronico dei fotoni secondari diffusi 


per effetto Compton: 
1 


* dw 
y, W/w, t)k(w) We? 


wie‘ 


i quali tendono ad aumentare la componente fotonica alle bassissime energie, 
& pure contenuta nella funzione &(H, À, n) e si manifesta nella correzione 


2Gu s + 1 


W (s +1)(8 +3)’ 


WW) = x(W) — 
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sulla sezione d’urto totale Compton, cosicche nel denominatore della funzione 
integranda compare l’espressione A + o(W) + x(W), in cui il parametro A 
deve intendersi variabile, come è richiesto dalle argomentazioni relative alla 
influenza dei poli dovuta al processo moltiplicativo. La dipendenza della 
E(E, A, n) e della O(H, À, n), pure definita in (4), dal parametro A & molto sen- 
sibile soprattutto nella regione delle basse energie; per questa ragione l’appli- 
cazione del teorema del prodotto integrale ha richiesto la valutazione di tutti 
i termini che definiscono le derivate seconde rispetto a 2 delle funzioni p(0, À, s) 
GOs ARO 
Esprimendo le energie in unità me?, si trova: 


B E BE 
I, = 0,0082 ee EEE 
Jo No Mo No 


E 
F(E,, 0) = 9,5012 — , 
No 


e quindi il coefficiente della probabilità differenziale D(x’)dx’, espresso in unitä 
Moliere, risulta: 


M=— 0,153(4). 


Un confronto abbastanza significativo fra le varie teorie della cascata: 
elettrofotonica proposte per le basse energie pud essere eseguito in base al 
valore di M che consentono di calcolare, dato che tale valore dipende in modo 
molto critico della maggiore o minore adeguatezza della schematizzazione uni- 
dimensionale della cascata. 

Con un calcolo piuttosto laborioso si ricava, in base alla teoria di BHABHA- 
CHAKRABARTY (°), considerata nella forma più rigorosa di SNYDER (?°), che il 
coefficiente di Moliere ha il valore M= 0,093, mentre, fondando il calcolo pre- 
cedente sullo spettro integrale di« Approssimazione B » di Rossi-Greisen invece 
che su quello di Clementel, si trova: 


H,(1) NENG KGS, Ss) NEC) ASE Sue 
ie oe ee 1 Le RE te 
i 34/(1) In (1832-18) — 73(3) — A3(3) areas! 


Infine, in base al valore di {Ry calcolato da ROBERG-NORDHEIM, corretto 
in modo da tener conto di tutti gli elettroni di energia maggiore di zero e non 
solo di quelli che possiedono energia maggiore di 0,058, si trova M= 0,1412. 

In conclusione, il valore trovato M= 0,153 sembra essere abbastanza plau- 
sibile, perché il procedimento delineato, in quanto consente di valutare con 
una certa accuratezza il numero di elettroni generati sotto l’energia critica 
e tiene conto dell’espressione corretta della traccia integrale, non possiede, 


SULLA PROBABILITA DI LOCALIZZAZIONE DEGLI ELETTRONI ECC. 833: 
oltre all’approssimazione b) del n. 1, quelle approssimazioni per difetto e rispet- 
tivamente per eccesso che sono implicite nei caleoli di MOLIÄRE e di Evans. 


Ringraziamo i proff. N. DALLAPORTA e E. CLEMENTEL per le utili discus- 
sioni e la dott. S. CHERSI che ha eseguito i calcoli numerici. 


SUMMARY 


The radial distribution of electrons in air showers is studied on the basis of the 
cascade theory of the photoelectronie component worked out by the Padua group. 
It is found that the probability of a singly scattered electron to be localized in the an- 
nular ring (x’, dx’) (#' in Molière units) is D(a’) dx’ = 0,153(dx//x°). This 


result is 
compared with the predictions of other theories. 
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Su una nuova equazione ondulatoria 
per una particella a spin 1/2. 


P. CALDIROLA e P. GULMANELLI 


Istituto di Scienze Fisiche dell Universita - Milano 


(ricevuto il 15 Luglio 1951) 


Riassunto. — Si confrontano i risultati che si ottengono da una nuova 
equazione ondulatoria per una particella a spin 1/2, recentemente pro- 
posta da H. C. CorBEN in una teoria unitaria, con quelli che derivano 
dall’equazione di Dirac. In particolare si tratta il problema della strut- 
tura fina dei livelli dell’atomo idrogenoide e si esamina il significato che 
nella teoria di Corben assume il momento magnetico (normale e anomalo) 
di una particella. 


1. — L’equazione di Corben. 


H. ©. CoRBEN (!) ha recentemente sviluppato una riformulazione della 
teoria quantistica dei campi ricorrendo all’uso di un formalismo in uno spazio 
a cinque dimensioni secondo il quale la carica elettrica di una particella é intro- 
dotta come operatore anziche come numero ordinario. La teoria di Corben 
risulta sostanzialmente equivalente all’ordinaria teoria dei campi; l’unica dif- 
ferenza di rilievo consiste nel fatto che l’equazione ondulatoria per una parti- 
cella di spin 1/2 non coincide esattamente con quella ormai tradizionale di 
Dirac. La nuova equazione di Corben per una particella sottoposta all’azione 
di un campo elettromagnetico pud essere scritta nel seguente modo: 


[ ie ie 
if 311% SAN + 4 =), 


(7) H. C. CorBen: Nuovo Cimento, 9, 580 (1952);. 
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dove i y; sono i consueti operatori di Neumann 


Ve = — Par (K=1,2,3), ya—=B, = Yıyıyaya = ir = — 0, ; 
e nella quale 


e 
(2) b= m,c? ’ 
essendo m, la massa che la particella viene a possedere pel fatto di essere dotata 
di carica elettrica e, mentre 
76 = 


MC 
h ’ 
essendo m, la massa della particella «nuda ». 

Si vede chiaramente, tenendo presente la (2), come l’equazione (1) diffe- 
risca da quella di Dirac unicamente per il termine in y,. In particolare, qua- 
lora si interpretasse tutta la massa osservabile di una particella dotata di carica 
come dovuta a questa, l’equazione di Corben diverrebbe: 


(3) in a + eV — a x (ihe grad — eA) + o,m,c?| y =0, 
da confrontarsi con quella corrispondente di Dirac: 
in = ol eV — ax (ihe grad — eA) + pane 20% 


In questo caso la differenza fra le due equazioni consiste dunque nella sosti- 
tuzione nella prima del termine o,m,c? al termine fm,c?. 

Si pone pertanto il problema di studiare se una siffatta modifica della equa- 
zione di Dirac porti o no a qualche conseguenza suscettibile di verifica speri- 
mentale. In relazione a questo scopo ci & parso utile studiare due problemi 
particolari: a) il moto di una particella in un campo centrale e, in particolare, 
coulombiauo, b) il momento magnetico di una particella e le sue correzioni 
relativistiche. 


2. - Moto di una particella in un campo centrale. Atomo idrogenoide. 


Ci proponiamo anzitutto di calcolare gli autovalori e le autofunzioni per 
una particella a spin 1/2 che soddisfi alla (1) e che si muova in un campo cen- 
trale caratterizzato dalla funzione potenziale U(r) = eV(r). Per amore di gene- 
ralità supporremo che siano presenti entrambi i contributi m, e m, alla massa 
della particella. 


55 - Il Nuovo Cimento. 
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L’equazione degli stati stazionari si scrive pertanto: 


(4 


(4) Wu — lea Xp + Pe? + Ex + eV(r) 


u, 


essendo W l’energia totale relativistica della particella. 
La (4) equivale al seguente sistema di equazioni differenziali: 


© 


4 mc + U—W en FIRE dus @ 5 0 
h € +" (Da, Jeo te 008 ay 
I none Se (i W d ee dur e 
: SN thee 4 Y= 
h ce Marl & ee dæ, Abe * 3 
à — moe + U—W D u, e 

= A Seh 0 
h ¢ et 5 : : Vz ann dx, hbe “a 4 
i me + U—W RAT RE NT dl ae 
h € Re Re de 


Si pud facilmente verificare (?) che il sistema ammette soluzioni di due 
tipi in corrispondenza di ogni singolo autovalore W— Win, l,j) caratterizzato 
dai tre numeri quantici n, |, j. 

Soluzioni di Tipo I: 


u =iFY", +il+m+1HYr, w= iF Y41—i(l—m)HY™, 


Ug =(l+m+1)GYr + KY", , Ug = | CS) G Mi Kaye 


Soluzion di Tipo IL: 


u = i(l + m)F'Y", +iH'Yr, u, = —i(l— m — 1)F'Yr4} +iH Yes, 
us = G'YP +(l+m)E'Yr,, Uy = @ Ir — ((—m—1)E'Yr4. 


Ci limiteremo ad esaminare esplicitamente solo le soluzioni di Tipo I, po- 
tendosi verificare immediatamente che i due tipi di soluzioni conducono & 
risultati equivalenti (si passa dalle soluzioni di Tipo I a quelle di Tipo II sosti- 
tuendo J con — (1 +1)). 

Le Y, che compaiono nelle espressioni precedenti, sono le ben note fun- 
zioni sferiche, mentre F, G, H, K sono funzioni di r che soddisfano al seguente 


(?) G. PerraAu: Journ. de Phys., 12, 810 (1951) ha usato un metodo analogo a 
quello da noi seguito per integrare le equazioni di Dirac per una particella sottoposta 
all’azione di un campo pseudoscalare. 
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sistema di equazioni: 


1 me +U—W ù l e 

Er N Co a Fr AL Te Là 
ob mea US W d 112 e 
= a | D het 
1 mee U == 4 d !+2 e 
annee né Rire ee 


1 — moo? + ane | ec ip 
h € 


Tale sistema, ponendo: 
ls VEN; I = Hie Kein 
ei 
k = vof, h = v9 , 


puö essere ricondotto alla seguente forma piu semplice: 


d Uses 1 moc? + U—W e = 
dr r “| h € Male ; 


EI 1 me + U—W N 
: ar r ee cc 7 AL oS 
(5) leat 
à +1 1 f—mc? + U—W e 
ù +1 il RP UW e 
a 3 al € il alt =o 


Le condizioni di compatibilità impongono che le funzioni Uy € v, Verifichino il 
seguente sistema di equazioni algebriche: 


Up + Vo = 0, 


W 
i (L + Uo) + M¢(Uy) — Vo) + vs (to + %) —0, 


da cui risulta: 


: MoC2b [mc2b\2 
(6) Uy = — Vo, =" = | ad 


e 


Per di più una discussione sui segni, che non riportiamo, porta alla scelta 
del segno positivo davanti al radicale. 
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Con cid il sistema (5) si riduce al seguente: 


d =, Me ei e | = 
(= r Jo % e a LE 
(7) 
d i+1 1 /— mc? + U— W an 
tr = Cc a 
Nel caso che il campo sia columbiano (V(r) = — Ze/r), per eseguire l’inte- 


grazione si pud seguire, come nel caso corrispondente delle equazioni di Dirac, 
il ben noto metodo di Darwin. Si arriva cosi alla seguente espressione per gli 
autovalori: 


Z?o? —1/2 


à (n! + Vj? — Zot? : 


2 


(8) W= (met) + > 


ove j=1+1/2 ex è la costante di struttura fina 

Come si vede l’espressione degli autovalori coincide completamente con 
quella che si ottiene risolvendo per lo stesso problema l’equazione di Dirac, 
pur di interpretare l’espressione: 


SS 
m = | me + pa 
come massa effettiva della particella. 

Si osservi esplicitamente che anche nel caso che la massa della particella 
fosse completamente legata alla carica elettrica, onde m, = 0 em = e/be? = m,, 
l’espressione degli autovalori per l’atomo idrogenoide è identica a quella della 
teoria di Dirac. Va perd notato che, come risulta dalle (6), le autofunzioni del- 
l’equazione di Corben differiscono da quelle della corıispondente equazione di 
Dirac; nel caso particolare di una massa tutta legata alla carica, si dovrà porre 
nelle (6) m, = 0 onde risulta », = 1 anzichè v, = 0 come nella teoria di Dirae: 


3. — Il momento magnetico della particella a spin 1/2. 


Come secondo problema ei siamo proposti di calcolare il momento magne- 
tico di una particella che soddisfi all’equazione di Corben, considerando sia il 
caso che essa si muova sotto la sola azione diun campo magnetico costante, 
sia il caso che essa si muova sotto l’azione di un campo di forze centrali (non 
necessariamente elettromagnetiche) a cui si Sovrappone quella di un campo 
magnetico costante. 
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Particella libera in un campo magnetico costante. — Supporremo il campo 
magnevico di intensitä costante H, diretto come l’asse z ma orientato in senso 
opposto. Il calcolo é analogo a quello eseguito da RABI (#) per un elettrone 
di Dirac e da uno di noi (‘) per una particella dotata anche di momento 
magnetico supplementare (termine di Pauli). 

L’equazione di Dirac-Pauli per il problema in considerazione si scrive: 


W € CA, 
(9) u a Xp + Po + be & a x en ee 
con: 
S=—fpo e o=—iamen. 


Ricorrendo a coordinate polari le equazioni esplicite, equivalenti alla equa- 
zione spinoriale (9), si scrivono: 


h € ORTE 26 h be de 
le me +) at exp Eig 2 + Ger) pot te = 0 
a i = lu +exn in is ger) ; Uy er 


Questo sistema ammette una soluzione particolare, che € quella di interesse 
fisico, del tipo: 


Un Er n)lexp ? (Da 22982) 


. i 
inp + 7 pe 


dove p & una costante qualsiasi e rappresenta fisicamente l’impulso della parti- 
cella parallelo al campo, e gli n, numeri interi che risultano legati dalle se- 
guenti relazioni: 


Ny = N,—1, Ns — "Na 1, Ne = Ny « 


(3) I. I. Rapı: Zeits. f. Phys., 49, 507 (1928). 
(4) P. CALDIROLA: Rend. Ist. Lomb. Scienze, 77, 173 (1943-44). 
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Ponendo ancora: 


eH,n 

P,(r) = exp |- = i ie Od AGA 
si hanno le seguenti equazioni per le R,(r): 
a W 2H). ‚fi Cam dk, 
| 2 + MC r er + (or + pe) + MB +r Ar =0, 
+e EEO ate ys LHe yap, yp Ae _ 
7 | L Mo Fes i p ibe R,— (m;—1)R; PRE str ar 
i Ww a dk, 
j (— ne Mo à rk, = (5 Po, Je, ++ Mf, + ar — 0 
a W 24, DE ARE leH, , dk, 
= in EE ; rk, (5 + rl (N, Er h oa I +r de —10 


Si possono integrare queste equazioni col solito sviluppo: 
p q P 


Binz CHE 


s>0 


e, senza difficoltà, si trova che si hanno soluzioni accettabili solo se & verifi- 
cata la condizione: 


ME ; (= me + = 0 i (i at al 2N 
| 0 i (4 = Me — =) ee ip — FS 
| h | = a a ; | Lure =| 2 
a ip +2 0 (TT me + 1 
e be 


con N — 0,1, 2,2% 
E questa un’equazione biquadratica in W che risolta dä i seguenti quattro 
valori possibili dell’energia: 


W, = + V Meet + 12H? + 2N Hyehe + 2AH men , 
eo W, = — V Mie + 22H? + 2NH,ehe — 2AH mc , 
W, =— V Mic: + AH? + 2NHeñc + 2AH me*n , 


Wy = — V Me + 2H? + 2N Hehe — 2AH ym ecn , 
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Ove Si è posto: 
e? | NHeh 
201 — m°ei 9262 Ener 2 p2 a | 0 
Me mic! + pe Te pa mie + pc? , n=Vı+ st, 
0 


x 


e il cui significato fisico è immediato. 
Si noti che nell’approssimazione newtoniana (p — 0) si ha: 


ME V/mgor + & + (20 a + 1" Je 


Arc \/ m} ay % m + 


per cui la particella verrä ad essere dotata del momento magnetico seguente: 


eh m 
11 IN = oe 
20) Ge Arme ie m 


che puö essere scomposto in tre parti: 


eh eh m 
ate 2N +1 AL e 
(Tt) ( u) Arme ae Anme 2 m ’ 


dovute rispettivamente al moto di precessione della particella (carica) attorno 
alla direzione del campo, al momento magnetico intrinseco legato allo spin e 
alla carica della particella di Dirac e al termine supplementare di Pauli. Si 
noti che quest’ultimo nella teoria di Corben risulta ridotto, rispetto al valore 
che esso ha nella teoria di Dirac, nel rapporto m,/m: in particolare nell’appros- 
simazione in cui ci siamo posti, per la quale ehH,/Anme< M,c?, esso si annulla 


\ 


se la massa della particella à completamente legata alla carica. 


Particella in un campo centrale e campo magnetico costante. — Calcoleremo 
ora il momento magnetico di una particella che soddisfi all’equazione di Corben 
(con termine supplementare di Pauli) e che si muova in un campo di forze 
centrali (non necessariamente elettromagnetiche) e sotto l’azione del campo 
magnetico costante H, orientato come nel problema precedente. 

Il calcolo é analogo a quello eseguito da MARGENAU (5) per una particella 
di Dirac (senza momento magnetico supplementare) e da uno di noi (°) per 
una particella di Dirac-Pauli. 


(5) H. MARGENAU: Phys. Rev., 57, 383 (1940). 
(5) P. CALDIROLA: Phys. Rev., 69, 608 (1946); Nuovo Cimento, 4, 89 (1947). 
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L’equazione fondamentale del problema é la seguente: 


(12) Nae ne axp + Bme + 2 0 + Alu, 
€ be 
ove si è posto: 
e A 
(13) A=—~ (axA) + (SxH). 


L’equazione (9) si integra col metodo delle perturbazioni considerando A come 
termine perturbativo. In approssimazione zero (A = 0) si è ricondotti a risol- 
vere l’equazione (4). L’azione della perturbazione provoca una variazione del- 
l’energia Win, l, j) di uno stato stazionario data da: 


AW= e u >, urn, 1, 5) Au.(n, 1, j)dr = 


= — — | [ui — iy)u, — usi(a + Üy)us + site — iy)us — uzi(~t iy)u, |dT— 


— AH, ff Quiu, — U3 Uy — Ug Uy + usu,] At (dt = r? sin Ÿ db dy dr). 


Passando à coordinate polari ed eseguendo calcoli che non riportiamo per esteso, 
si perviene al seguente risultato: 


«eH, 16x(1—m)! (1 + m +1)! (1 + 1)(2m +1) 


va (21 + 1)(21 + 3) : 
3 1 
; 6 [ I DE ( j 5)» 
(14) Rose ne 
hie (mee a n) 
a q ities ! 22) Le ap. ASUS | 
AH 47 + m + 1)! (— m)! (2m + 1)(06 — 1) EP eae 


dove: 


x ne; OP Ou = [er dr, 
J i 


0 


sono grandezze legate fra loro dalla condizione di normalizzazione, che si pud 
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scrivere sotto la forma seguente: 


1 = (@ + y)4n(l + m +1)! (1—m)!(1 4 0), 


e dalla seguente relazione, deducibile direttamente dal sistema di equazioni 
radiali, 


| ePer2 dr = | (2622) G2721d7e 
0 0 
essendo: 
: Bein 
ge — 
e2 
V mer): + 7, 


E facile trovare per AW delle espressioni particolarmente semplici nei casi 
asintotici © >c e v —0. 
Nel primo caso si ha: 


«(t+ 1)(2m +1) ehe 
See vo) on tes os 
Aoc|/ (myer)? +7, 
(! + 1)(2m +1) MC? 
2 ES Fe fe 
RG DB a) RE 
| (mer): + 
e nel secondo: 
© 2 = NE 2 
AW = (l ss EN 2 a 2) Ich, x ehe + = at 1 15; Bite = mn 
(21 + 1) Piece un = 
x| (mac)? +, (moe?) + pe 


Corrispondentemente si hanno le seguenti espressioni per il momento magne- 
tico della particella: 


per ov ie: 
eee feat! ehe : MC? i) 
(15) jes (i ate 1 er. + (29+ 1) Be, i 2G +1) 
An |/ (moc2}? + |/ (mget)? + © 
b? b? 
per v<e: 
1 
jee 
1 2 _ eh 
(16) = (+50, ) ru 
j +1 4 e? 
7 | (moc?)? + be 
5 MoC 
1 — TE = —— À ; 
+ 25 +1 | e 
(moe?) + pe 
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quest’ultima, nell’approssimazione newtoniana (v = 0), si riduce a: 


1 2 


Si noti che, come nel caso precedente, nell’approssimazione in cui 4H,< moc?, 
se la particella ha massa interamente legata alla carica il termine in A si an- 
nulla; viceversa se la massa è tutta di origine non elettrica si ha w= A. 


Osservazione — E noto che il momento magnetico del protone puö essere scritto 
sotto la seguente forma: up = wy + Ap; essendo u, = eh/4nm,c il magnetone nucleare 
© Ay = 1,793u,; mentre quello del neutrone & dato da uy = Ay con A, = — 1,913w. 
D’altra parte se riteniamo che la massa del protone sia di origine parzialmente elet- 
trica, si dovrà avere per la (17): 


(18) Ap = À ——— 


con VU? = M° — 1836m,,. Si pud allora porre il problema di esaminare se & pos- 
sibile interpretare la differenza tra il momento magnetico supplementare del protone 
e il momento magnetico del neutrone (o meglio dell’antineutrone) come dovuto all’esi- 
stenza nel protone di una massa parziale legata alla carica. Supposto pertanto 
À = — Ay, Si puö determinare il valore di M, che, sostituito nella (18), porti al valore 
osservato di À. Si trova M, ~ 203 m,,. Puö essere interessante notare come il valore 


della massa protonica legata alla carica sia circa eguale al valore della massa del 
mesone u. 


Conelusioni. 


I calcoli, eseguiti nei paragrafi precedenti, mostrano come gli autovalori 
per l’atomo idrogenoide che si ricavano dall’equazione di Corben coincidano 
(previa normalizzazione della massa) con quelli deducibili dall’equazione di 
Dirac. Questo risultato é del tutto generale, in quanto si verifica immediata- 
mente che, per una particella non dotata di momento magnetico supplemen- 
tare, l’equazione del 2° ordine che si ottiene iterando l’equazione di Corben 
coincide con quella che si ottiene iterando l’equazione di Dirac. Esistono perd 
differenze nelle espressioni delle autofunzioni, le quali possono, in qualche 
caso, avere interesse fisico. 

Differenze anche negli autovalori delle due equazioni si hanno invece se 
alla particella a spin 1/2 si attribuisce un momento magnetico supplementare, 
introducendo nelle equazioni un termine di Pauli. In questo caso infatti le 
proprietä magnetiche della particella dipendono in modo notevole dalla natura 
della massa della particella stessa. 
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Uno dei due autori (P. G.) ringrazia il Consiglio Nazionale delle Ricerche 
per il conferimento di una Borsa di studio che gli ha reso possibile la perma- 
nenza presso l’Istituto di Scienze Fisiche dell’Universitä di Milano. 


SUMMARY 


The results obtained by a new ondulatory equation for a spin 1/2 particle, recently 
proposed by CORBEN in an unitary field theory, are compared with those araising from 
the ordinary Dirac equation. Particularly the problem of the fine structure for hy- 
drogen like atom levels is considered and the relation between the magnetic moment 
and the nature of the particle mass is examined. 


IL NUOVO CIMENTO Vor, EX. N.9 1 Settembre 1952 


Sull’operatore U(t) di Dyson- Feynman. 


S. FUBINI 


Istituto di Fisica del Politecnico - Torino 


(ricevuto il 17 Luglio 1952) 


Riassunto, — Si dimostra che l’operatore U(t) di Dyson-Feynman diago- 
nalizza la Hamiltoniana totale anche nel caso in cui esso non & definito 
mediante il metodo delle perturbazioni. L’operatore U(t) è perd, in gene- 
rale, solo unitario in un senso ristretto, il che ci permette di estendere 
le nostre considerazioni anche al caso in cui esista uno spettro discreto 
di autovalori della Hamiltoniana. 


1. — Introduzione. 


E stato dimostrato da FERRETTI, in una serie di suoi recenti lavori (1), che 
l’operatore U(t) di Dyson-Feynman genera una trasformazione canonica che 
diagonalizza la hamiltoniana totale del sistema. 

In questo lavoro ci proponiamo di estendere le considerazioni di FERRETTI 
al caso in cui non sia valido un calcolo di perturbazioni. 

La dimostrazione che daremo si basa sull’ipotesi adiabatica, la quale à 
indispensabile nella forma attuale della teoria dei campi, poiché senza di essa: 
non è possibile definire l’operatore U(t) di Dyson-Feynman (?). 

Infatti l’equazione (5): 


dU( 


t) s 
(1) ET iV(t)U(t) = 0. 


(1) B. FERRETTI: Nuovo Cimento, 7, 79, 783 (1950); 8, 109 (1951); F. J. Dyson: 
Phys. Rev., 75, 486, 1736 (1949); R. P. FEYNMANN: Rev. Mod. Phys., 20, 367 (1948). 

(?) B. A. Lippmann e J. SCHWINGER: Phys. Rev., 79, 469 (1950); H. SNYDER: Phys. 
Rev., 83, 1154 (1951). 

(5) Nel corso di questo lavoro useremo un sistema di misura in (ui sia 7 = FE 
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assieme alla condizione iniziale 
(2) Dos) AN 


€ equivalente all’equazione di Volterra: 
t 

(3) ete | V(t!) U(t!) au’. 
—00 


L’equazione a secondo membro ammette soluzioni solo se l’integrale (nel 
senso di Riemann) al secondo membro esiste; questo avviene solo se per 
t > — oo, l’integrando tende a zero più velocemente di 1/4. 

E percid opportuno, per definire l’operatore U(t), ricorrere al seguente pro- 
cedimento: si moltiplica V(t) per un fattore di convergenza A(«, t) tendente, 
per ¢ > — oo, a zero più in fretta dil/te tale che lim A(«, t)=1; lim dA(e, t)/dt —0 

a—o0 


a—> 0 
(si potrà, per esempio, assumere A(a, t) = exp [— at®]); si definisce quindi un 
operatore U(«, t) mediante l’equazione : 


dU(«, t) 


aie Aa, t) V(t) O(a, t) = 0, 


(4) 


e la condizione iniziale: 
(5) Ua, — ©) =1. 
Per valori di t finiti si fa tendere « a zero, e si pone: 


(6) (EKG) nen Wee OY 


a>0 


2. — Proprieta dell’operatore U(a, t). 


Le considerazioni che faremo si applicano a tutti i casi in cui l’energia di 
interazione V, nella rappresentazione di Schrödinger, non dipende dal tempo. 
Esse si applicano perciö alla teoria quantistica dei campi localizzati, in cui 
la Hamiltoniana di interazione dipende dal tempo solo attraverso gli operatori 
dei campi. 

L’interazione V(t) che compare nelle (1), (3), (4) € espressa nella « rappre- 
sentazione di interazione » e si ricava dalla V della rappresentazione di Schrö- 
dinger mediante la trasformazione canonica: 


(7) V(t) exp [cHt]V exp [— cH], 
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da cui segue l’equazione: 


dv 
(8) ae = ı[E V(t)]. 


Sia U(a, t) l’operatore coniugato hermitiano di U(x,t), esso obbedisce alle 
equazioni: 


(4') ee Ue A(x, t) (a, t) V(t) = 0 
(5') Ua, — co) =0, 
(6') lim U(a, t) = U(b). 


Dalle equazioni (4) e (4’) si ricava: 


d 
(9) a (Um NU (a, 1) =0, 


da cui, ricordando le (5), (5’), (6), (6'): 


(10) U(a, t)U(a,t) =1, 


(11) DUG We 


Osserviamo che non si pud dimostrare allo stesso modo la relazione 
U(t)U(t) =1 (4) la quale non é in generale vera. Questa osservazione ci sara. 
utile in seguito. 

Dalle equazioni (4) e (8) si ricava: 

COUNT) gg (La) 


| dV (t) 
de 37 Mlle cde lt 


V(t)U(a, t) + Aa, t) —. — Ula, t) + Aa, t)V( 


dU (a, t) 
u | 


ne 


- CE U(a, t)] = — À(x, t) A) U(a, t) — 1A(x, t)V(t)[E Te, t)]. 


(12) a: nt Ua, Dia Cat), 


(4) M. Neuman: Phys. Rev., 83, 671 (1951). 
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si ottiene: 

d 
(13) N erden ir, À out, t), 
(14) C(a, — ©) — 0. 


Per & — 0, tende a zero il secondo membro della (13); quindi dalle (13), 
(14) per ragioni di continuità Segue: 


(15) lim O(a, t) = 0. 


a— 0 


3. — L’operatore W. 

Consideriamo l’operatore: 
(16) We, t) = exp [--ıEt]U(e, t) exp [Et]. 
Derivando si ottiene: 


(17) ee = exp [—ıHt]0(«, t) exp [ıHt], 


da cui, ricordando la (15): 


(18) lim Wa, t) =W, 


a>0 
dove W (5) & un operatore indipendente dal tempo. 
Ci proponiamo ora di dimostrare che l’operatore W diagonalizza la hamil- 
toniana totale 7 
HWat)=(E+VWWeat)=(E + V) exp — ıEt]U(«, t) exp [ıEt] = 
= exp [— «ILE + V(t)]U(«, t) exp [ct] = exp [— ıEt]U(«, t) exp [ıEt]E + 


ae ea = exp [— Ht] 0e 2 exp [cHt] + exp [— ıEt][E U(«, t)] exp [Ht] = 
x, t) dt 
Wo Bie m exp [— cHt] { C(a, t) + [1 — Aa, t)]LE U(a, t) )] } exp [ıEt], 


Ale, t) 


() L’operatore W da noi definito con metodi dipendenti dal tempo coincide con 
l’operatore Y definito da C. MÖLLER con metodi stazionari. Dan. Mat. Fys. Medd., 
23, 1 (1945). 
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facendo tendere « a zero per valori di ¢ finiti si ottiene: 
(19) HW= WE, 

e prendendo i coniugati hermitiani di ambo i membri: 
(19°) WH= EW. 


Dalla (19) segue che l’operatore W trasforma una autofunzione di Æ in 
una autofunzione di H corrispondente allo stesso autovalore. Infatti se ®, 
è una autofunzione di # corrispondente all’autovalore Æ, dalla (19) si ottiene: 


(20) HW®, = E,W®,. 
Dalla (10) segue: 
(21) ws 


cioé il sistema delle W®, & ortonormale. 

Ci si pud domandare se il sistema delle W@, sia completo. La risposta dovra 
essere negativa poiche H puö avere degli autovalori (livelli energetici dei si- 
stemi legati) che E non ha. 

Possiamo renderei conto di questo mediante le seguenti considerazioni: 

Poniamo: 


(22) WW—E, 

per la (21) deve essere: 

(23) et = €; 

da cui segue che € ha per autovalori 0 ed 1. Dalle (19), (19’) segue inoltre che: 
(24) [en 0% 


Esiste percid un sistema completo di autofunzioni comuni ad € e ad H. 
Le autofunzioni W®, sono autofunzioni di € corrispondenti all’autovalore 1. 
Infatti: 


(25) EWD, = WWW®, = W®,. 
Le autofunzioni comuni di H e di € corrispondenti all’autovalore 0 cor- 


rispondono percid allo spettro discreto di autovalori di H, mentre le auto- 
funzioni corrispondenti all’autovalore 1, corrispondono allo spettro continuo. 
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Ringrazio il prof. G. WATAGHIN per il suo cordiale interessamento e per 
le utili e chiarificatrici discussioni avute sull’argomento. 


SUMMARY 


The Dyson-Feynman operator U(t) is shown to diagonalize the total Hamiltonian 
independently from any perturbation-theoretical treatment. The fact that, in general, 
U(t) is unitary only in a restricted sense is used to extend this study to the case in 
which the Hamiltonian admits also discrete eigenvalues. 
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LETTERE ALLA REDAZIONE 


(La responsabilita scientifica degli seritti inseriti in questa rubrica è completamente lasciata 
dalla Direzione del periodico ai singoli autori) 


Sulla capacita di un condensatore sferico. 


F. BERTOLINI 


Istituto Nazionale per le Applicazioni del Oalcolo - Roma 


(ricevuto il 15 Giugno 1952) 


La capacità ec, di un condensatore sferico, di cui l’armatura interna abbiæ 
1% 0 
raggio r e quella esterna raggio R, è data dalla nota formula 


2 > = +RIR—r), 


nel sistema C.G.S. elettrostatico di unita. In pratica, l’armatura esterna di un con- 
densatore sferico presenta un’apertura, onde permettere il caricamento dell’armatura 
interna, e quindi la (1) deve considerarsi affetta da un errore: vogliamo stabilire 
una limitazione (inferiore e superiore) di tale errore. Essa si consegue facilmente, 


utilizzando una nota proprieta di media delle funzioni armoniche in uno strato 
sferico (1). 


1. — Si abbia uno strato sferico di centro O, raggio interno r e raggio esterno R; 
sia w la superficie sferica di centro O e raggio 1, dw(Q) elemento d’area di questa 
nel punto Q; per individuare un punto dello spazio, ne daremo la distanza o da O e la 
proiezione P da O sulla superficie w; sia w(o, P) una funzione armonica nello strato 
sferico accennato. Secondo l’indicata proprieta di media, riesce 


(2) u(o) = = | wo IZ GWayG@2) = Gs 22 = per rS oS R. 


wo 


dove a e b sono costanti. 
Ne segue 


b 
SER 


? 


Sr 
u'(e) jme P)da(P) 


@) Cfr. M. PICONE: Appunti di Analisi Superiore (Napoli 1940), p. 96. 
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e quindi 
1 
AT. 


[0] 


(3) = | 0”u.(0, P)do(P) = b, per r<o<R, 


D’altra parte, dalla (2) seritta per @ — r e per 0 = R si ricava 


ee 
(4) D | [w(r, P) — u(R, P)]do(P) , 
L TC 


e quindi (dal confronto delle (3) e (4)) il teorema: 

Il’ flusso di una funzione armonica in uno strato sferico, attraverso una super- 
ficie S regolare semplice chiusa contenuta nello strato e contenente nel proprio interno 
il centro dello strato, vale 


| "du rR 
d 


(5) oo | [u(R, P) — u(r, P)]do(P) = 42 a [u(R) — u(r)]- 


dn, IE, 
Ss w 
La (4) e la (5) possono anche dedursi (assai meno rapidamente) dalle formole 
che danno una funzione armonica in uno strato sferico, della quale siano noti i 
valori sulla frontiera dello strato stesso (2). 


2. — Venendo ora allo scopo di questa nota, sia O° il condensatore sferico ideale, 
costituito di una sfera 7, di raggio r e centro O, e di uno, strato sferico concen- 
trico 7%, di raggio interno R (> r) e spessore arbitrario; sia 8, la superficie della 
sfera 7,, 8, la superficie interna dello strato T3. Sia © il condensatore sferico 
reale, costituito della sfera T,, e del settore di strato sferico T, ottenuto asportando 
da 73 tutti i punti interni ad un cono circolare di vertice in O e semiapertura 6; sia 
Z la traccia di tale cono su 8,, sia y la traccia di tale cono su o. 

Com’é ben noto, se portiamo 7, al potenziale 1, T, al potenziale 0, in tutti i 
punti dello spazio si genera un potenziale u(e, P), funzione armonica nell’insieme 
complementare di ©, e verificante le condizioni 


(6) ur, P)=1, uk, P)=0 su S,—2, limu(oe, P) = 0; 


0 >09 


d’altra parte, la capacità di © & data dalla formola 


ER NT : 
ee pean 
8; 
e quindi (per le (4), (5), (6)) 
(7 er N re te 
) WER ae uk, P)do(P)| ; 


y 
l’errore relativo da cui & affetta la (1), quindi, vale 


(8) “—- =— | u(R, P)do(P) = 5 [mcr y) sen pdp > 0, 


v 0 


(?) Cfr. M. PICONE: op. cit., Ds 371: 
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avendo indicato con ¢ la colatitudine del punto P in un sistema di coordinate polari 
di polo in O, avente per asse polare l’asse di simmetria di 2 diretto nel verso che 


AiO vata ale). 
Per ottenere delle limitazioni dell’errore relativo (c, —¢)/¢ bastera ora trovare 


delle limitazioni per u(R, 9). 


3. — Consideriamo la funzione (r/e) — u(e, P): essa vale zero su Sj, è positiva 
sulla frontiera di T,, infinitesima all’infinito, armonica nell’insieme complementare 
di ©, quindi & sempre positiva in tale insieme: ne segue, in particolare, 


(9') 0<u(R,p)<r/R per 0<p<ô. 
Consideriamo ora la funzione 


(r/o)? cos m — (r/o)(r/R) cos ö 
— = (0, P) : 
1 — (r/R) cos 6 
essa 6 non positiva sia sulla frontiera di T, che sulla superficie S,, infinitesima 
all’infinito, armonica nell’insieme complementare di 0, quindi 6 sempre non positiva 
in tale insieme; ne segue in particolare 


(r/R}?[cos @ —- cos à] = Bo) ey, 
- —_____— < w(R, er Ss N 
b= (r/R) cos Ô F2 P = 


(9) 


Dalle (8), (9’), (9/’) segue la limitazione cercata 


P)\2 gin4 ) Peer 
iR (02) BET eyesore 


(10) ne 
1 — (r/R) cos 6 Cp 


Da verifiche eseguite in qualche caso particolare, avendo determinato l’errore re- 
lativo sperimentalmente con il calcolatore analogico del prof. PERSICO, (¢ — c)/c, 
è risultato approssimativamente eguale alla media geometrica tra i due limiti 
indicati, e preeisamente (per r= 4, R = 18, n/4> à > x/15) una decina di volte 
maggiore di quello inferiore, una decina di volte minore di quello superiore (*) (*). 


(#) Data la simmetria assiale delle condizioni (6), anche la funzione w(o, P) godrà di tale 
sim metria assiale, ossia, nel sistema di coordinate polari scelto, sara indipendente dalla longitu- 
dine; se @ & la colatidudine di P, porremo (0, p) = wo, P). 

(4) Non ci @ noto che altri abbia trattato matematicamente il problema della determina- 
zione di (c, — c)/c,; tuttavia in una nota di E. B. Rosa ed N. E. Dorsey [4 new determination 
of the ratio of the electromagnetic to the electrostatic Unit of electricity, in Bullettin of the Bureau of 
Standards, vol. 3, nu. 3-4, pp. 481 e sgg. (1907)] questa grandezza viene determinata speri- 
mentalmente. 

(*) La Redazione del Nuovo Cimento ringrazia il Dr. CH. Moon Capo della Sezione Induttanza e 
Capacità della Divisione Elettricità del Bureau of Standards per aver fornito l’indicazione biblio- 
grafica riportata nella nota precedente (N.d. R.). 


IL NUOVO CIMENTO Vor. IX, N. 9 1 Settembre 1952 


Sulle soluzioni dell’equazione di Dirac-Corben. 


A. LOINGER 


Istituto di Fisica dell’ Universita - Pavia 


(ricevuto il 15 Luglio 1952) 


Per equazione di Dirac-Corben (1) intendiamo l’equazione 
(1) [Yp%, + haey, —iMec]y = 0, 


ove 


relativa a una particella di carica e, spin 1/2, massa m = (M? + a?e?h2/e2)U2. 

La (1) & stata applicata (?) allo studio di alcuni problemi di interesse fisico ed 
è stato fra l’altro osservato che l’equazione del secondo ordine che si ricava da (1) 
coincide coll’equazione del secondo ordine cui da luogo Vordinaria equazione di 
Dirac 


a (Yu = UNC) Wp = 0, 


In questa lettera vogliamo mostrare che se sono conoseiute le soluzioni y, del- 
l’equazione di Dirac (2), si puö facilmente risalire alle soluzioni della (1). 

A questo scopo passiamo a considerare l’equazione di Dirac-Corben del secondo 
ordine. Applicando a sinistra l’operatore y,x, + haey, + Me otterremo da (1) 
l’equazione 


he à 
(3) (air, Lo MP ae NR AE re) 07205 


. () H. ©. CorBEN: Nuovo Cimento, 9, 591 (1952). Le nostre notazioni sono quelle stesse di 
tale lavoro, al quale pertanto senz’altro rinviamo. 
(?) P. CALDIROLA e P. GULMANELLI: Nuovo Cimento, questo stesso fascicolo. 


856 x A. LOINGER 


\ 


che, come si & detto sopra, & identica all’equazione del secondo ordine di Dirac. 
La soluzione y dell’equazione del primo ordine di Dirac-Corben sara pertanto 


(4) y = (Vp, + haey, + iMe)®. 


Ma @ si pud ottenere da y,, soluzione (che supponiamo nota) di (2), assai facil- 
mente; infatti 


(5) Y = (Yu. + ime)®, 


e inoltre, considerando accanto a (2) anche la 


(27 Sr imC)Ph = 0, 
sarà 


(6) Pp = (Yup —ime)® . 


Da (5) e (6) si deduce subito che la soluzione ® pit generale dell’equazione del 
secondo ordine € una combinazione lineare delle y,, ,. 


E pertanto le soluzioni y,, 9, possono utilizzarsi per costruire, mediante la for- 
mula 


fete hae 
i ic M + m) u 


che si ottiene dalla (4) con qualche passaggio, la più generale soluzione dell’equa- 
zione del primo ordine di Dirac-Corben (1). 


IL NUOVO CIMENTO Vor. IX, N. 9 1 Settembre 1952 


Sull’anomalia nella distribuzione angolare intrinseca dei rami 
nelle stelle prodotte dalla radiazione cosmica. 


G. LOVERA 


Istituto di Fisica dell’ Universita - Modena 


(ricevuto il 28 Luglio 1952) 


In una precedente nota (1), & stata rilevata una anisotropia nella distribuzione 
angolare intrinseca dei rami di stelle prodotte dalla radiazione cosmica nelle emul- 
sioni nucleari, almeno per le stelle con piccoli numeri di rami. Precisamente, si 
6 osservato che le frequenze sperimentali degli angoli tra le proiezioni successive 
dei rami delle stelle sul piano dell’emulsione, si scostano, per le stelle a 3 ed a 4 rami, 
senza recoil e senza rami relativistici, dalle frequenze relative ad una distribuzione 
casuale, isotropa, e si lasciano invece inquadrare soddisfacentemente nella stati- 
stica degli eventi correlati, sviluppata nella nota stessa. 

Le formule (9) e (10) ivi riportate, sono casi particolari della 


1 


fra — EL di 


(1) PÉTER) Se 


1 
Fer] —tym-1dı 
F 
per 7 — 055; 1,0; 1,5; 2,0; 6 m= 2: 3. 

Lo studio statistico suddetto & stato esteso anche a stelle con n = 5; 6; 7 rami 
(non relativistici), osservate in lastre G5 Ilford esposte per alcune settimane a 
3500 m (di cui alla precedente nota), con i seguenti risultati. 

Per le stelle con n = 5 rami, prive di recoil visibile, la distribuzione presenta 
una netta discordanza dall’isotropia. 

Precisamente, per m = n — 1 = 4, dalla (1) si ottiene: 


Sr 
(2) r=05 PwW(I,4,1)=1— me (5 =, i), 


(2") r= P (1; 4,1) = (1— J, 


j 9 009 1 10 10 5 1 
(2") r= 1,5 P&:5(1; 4, 1) 1 Dee mn 1—12 + 7 13 ar D Pp — — r) , 


(2°) eh) PORN) ZI) se nr 


() G. LOVERA: Nuovo Cimento, 9, 134 (1952). 
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La formula per r = 1 corrisponde al verificarsi di una distribuzione casuale, 
di eventi non correlati (isotropia). Dalle (2')-(2/"), con le modalita indicate nella 
nota sopracitata, si calcola la probabilita pi (p?; pl) che l’angolo tra due rami suc- 
cessivi della proiezione della stella sia compreso tra py? e p! gradi sessagesimali (ta- 
bella I). 

Le frequenze sperimentali per le stelle a 5 rami (senza recoil) sono state deter- 
minate su 125 eventi (625 angoli): esse sono raccolte nella tabella I, e poste a con- 
fronto con le probabilità teoriche per r = 1,0; 1,25; 1,5 (i dati relativi ar = 1,25 
sono stati calcolati per interpolazione). 


TABBLLA I. — Frequenze sperimentali e probabilita teoriche per n = 5. 


PS (pis pa) 


Oe TPS Freq. sper. | 
fay eC) | ru 25 FPS) 
| | 
| O= 36° | 0,307 0,344 | 0,306 | 0,274 | 
| XG. || 0,268 | 0,246 | 0,275 | 0,299 | 
[22108 0,193 0,170 | 0,193 0,212 | 
| 108 = 144 0,110 0,110 | 0,115 0,119 | 
| - 144 = 180 0,082 0,067 | 0,064 0,061 
| 180 + 216 0,027 0,037 | 0,030 | 0,025 
|. BIG = Wis 0,013 0,018 | 0,012 | 0,008 | 
252 ~ 288 0,000 0,006; | 0,004 0,002 
288 + 324 0,000 0,001, 0,001 0,000 
324 = 360 0,000 0,000 | 0,000 0,000 


Appare evidente che le frequenze sperimentali si scostano nettamente da quelle 
previste per il caso dell’isotropia (r = 1), e sono, invece, in accordo soddisfacente 
con quelle calcolate per r = 1,25. 

Invece, le deviazioni sistematiche delle frequenze sperimentali dall’isotropia si 
attenuano, sino a diventare praticamente irrilevanti, per le stelle, senza recoil, con 
n = 6 rami (44 stelle, 264 angoli), e rispettivamente con n = 7 rami (34 stelle, 
238 angoli), come risulta dalla tabella II. Questo risultato appare sufficientemente 
netto, pur essendo minore l’entità del materiale statistico disponibile. Percid l’inda- 
gine non è stata estesa alle stelle con n > 7. 

Dai dati della presente e della precedente nota (!), si deduce che l’anomalia della 
distribuzione dei rami rispetto all’isotropia è particolarmente spiccata per le stelle 
senza recoil a piccolo numero di rami, e si attenua al crescere di questo sino a scom- 
parire praticamente per n — 7. Essa si estrinseca in una maggiore frequenza degli 
angoli appartenenti alle classi intermedie, a scapito degli angoli piccoli e grandi, 
rispetto alle frequenze inerenti alla statistica degli eventi casuali, indipendenti: si 
osserva cio® una simmetria nella distribuzione dei frammenti nucleari maggiore di 
quella comportata da una loro emissione puramente casuale. Siffatta anisotropia 
è peculiare delle stelle senza recoil: in altri termini, il fenomeno & presentato essen- 
zialmente dalle stelle provenienti dalla disintegrazione dei nuclei leggeri presenti 
nell’emulsione. 
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TABELLA Il. — Prequenze sperimentali e probabilità teoriche per n = 6, e n — 7 
| w—=6 nm = 7 
| | 2 | : | 
Pi + Ye Freq. | Pa? (915 93) Freq. 2 (P25 98) 
Sper. ss 1,0 es 1,25 | sper. = 1,0 
| | 
0 — 360 0,401 0,409 0,371 0,471 0,469 
36 — 72 0,286 0,263 | 0,300 0,265 0,269 
72 = 108 | 0,146 0,160 | 0,176 | 0,132 0,145 
108 — 144 | 0,098 0,090 0,088 | 0,111 0,071 
144 — 180 0,053 | 0,047 0,044 | 0,021 0,031 
180 = 216 0,008 0,021 0,015 | 0,000 0,011, 
216 — 252 | 0,004 0,008 0,005 | 0,000 0,003 
| 252988 | 0,000 0,002 | 0,001 | 0,000 0,000; 
| 288 = 324 | 0,000 0,000 | 0,000 | 0,000 0,000 
324 = 360 | 0,000 0,000 0,000 | 0,000 0,000 


Questi risultati inerenti alle stelle senza recoil visibile sono in accordo con le 
conclusioni di analisi analoghe, condotte con criteri statistici differenti, da DELLA 
CoRTE e GIOVANNOzzI (2) (limitatamente an = 3) e da Hopason (3) fino an = 6: 
i dati contenuti nella presente nota permettono di accertare la presenza dell’ano- 
malia almeno fino a n = 5, con evidenza ed efficacia descrittiva maggiori rispetto 
al eriterio del y? (adottato da Hop&sox). 

E, infine, opportuno rilevare che, per estrapolazione dei risultati ora esposti, 
anche per le stelle a due rami senza recoil e’è da attendere una distribuzione stati- 
stica angolare (degli angoli in proiezione) non coincidente con quella casuale. 

Dalla (1) per m = n — 1 — 1 si ottiene 


2 1 —7\12 
(3°) m 00 POS; 1,1) == acte | | ; 
Tt N 
(3°) Del BO (pail Ti) 1 Say 
2 DIN 2 
(3”) r= 1,5 HEE) (Uae) ant | = | + = (1 —71)/2702(1 — 21), 
T ‘ l 7 
(3””) T—2 P® (1; 1, 1) =3(1 — I? — 2(1 — 1}, 


e le corrispondenti distribuzioni di probabilitä per i dieci intervalli angolari sinora 
considerati sono contenute nella tabella III. 


; (*) M. DELLA Corre e M. GIOVANNOZZI: Nuovo Cimento, 8, 741 (1951). 
() P. E. Hopason: Phil. Mag., 44, 190 (1952). 


860 G. LOVERA 


TABELLA III. — Probabilità teoriche per n = 2. 


| | PP: 99) 
| FL + | ‚ | 
| | r=05 | r=10 | r=18 r = 2,0 
0 = 360 0,205 | 0,100 0,052 0,028 
Do 0,090 0,100 | 0,090 0,076 
724-108 0 0,074 0,100 | 0,110 0,112 
1082144 0 0,067 0,100 | 0,121 0,136 
144 = 180 | 0,064 0,100 | 0,127 0,148 
180 — 216 0,064 | 0,100 | 0,127 0,148 
216 + 252 | 0,067 | 0,100 0,121 0,136 
252 + 288 | 0,074 | 0,100 | 0,110 0,112 
288 + 324 0,090 0,100 | 0,090 0,076 
324 = 360 0,205 0,100 | 0,052 0,028 


Considerando i dati di quest’ultima tabella, appare evidente che la distribuzione 
statistica angolare potrebbe essere profondamente alterata se, invece di r = 1, ciod 
della distribuzione casuale, valesse, invece, ad esempio il caso r = 1,5 (come accade 
per n = 3, per le stelle prive di recoil). 
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Intensity measurement of Cerenkov light. 


P. Bassi, A. M. Brancur and C. MANDUCHI 


Istituto di Fisica dell’ Universita, Istituto Nazionale di Fisica Nucleare - Sezione di Padova 


(ricevuto il 12 Agosto 1952) 


According to the theory (1) the in- 
tensity of Cerenkov radiation is propor- 
tional to 1 — 1/n2ß2 so long as the de- 
pendence of n on the frequency is not 
considered. The dependence of the 
observed light on 8 may vary according 
to different experimental conditions: 
e.g. coherence might affect the results. 

In this letter we are publishing the 
first measurements taken with an ap- 
paratus the purpose of which is to 
study the intensity of Cerenkov light 
as a function of the energy of the inci- 
dent particle. 

The apparatus consists of four trays 
of G.M. counters I, II, III, IV and a 
Cerenkov light detector. Coincidences 
I, II, III, IV and anticoincidences HS TT 
III-IV can be registered. The detector 


consists of an EMI 5311 photomultiplier 
at the bottom of a plexiglass cilinder 
16 em high and 4 cm in diameter, coated 
with MgO. 

The coincidence or anticoincidence 
output signal goes to a syncroscope 
which triggers a C.R.O. (Phylips DB13/2 
the vertical plates of which are con- 
nected with the photomultiplier ampli- 
fier. The amplifier has a rise time equal 
to decay time (1,5 us) and a gain of 350: 
besides the usual advantages of these 
type amplifiers, a convenient delay of 
the peak value can be obtained. The 
pulses are photographed and their, height 
measured. 

The anode and heater supplies are 
stabilized. 

In order to select particles of diffe- 


TABLE I. 
| em Pb | cm Pb | cm Pb à ©. W 2 2 
| gjem? air | | 
| | 
| 22.5 0 15 40 — 100 141 22.46 0.70 8 
I IN 15 52 140 — 500 | 144 24.92 0582 10 
15 15 52 500- co | 250 | 32.18 0.70 10 


@) H. J. Som: Quantum Mechanics (New York, 1950). 
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rents B, various thicknesses of the lead 
absorber B and © were used. Absor- 
ber A has the purpose of increasing the 


plexiglass 


10 cm 


intensity in the interval selected by the 
counter system; it also absorbs secun- 
daries from the roof and the air. 
Results are plotted in fig. 2a and 
Table I where N is the number of photo- 
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Fig. 2. 


graph, W the observed intensity of ra- 
diation in arbitrary units, q and q re- 
spectively the standard error and the 
mean error of a single measurement, 
according to the Poisson distribution. 
In fig. 2a the average values of energy, 
in fig. 2b the theoretical curves of the 
radiation intensity W, of the emission 
angle 4 ans of the ionization energy 
loss J are reported, all of them in arbi- 
trary units. 

For high momenta our results do 
not agree with theory. 

Beyond some hundred grams per 
em? 9 and J are practically constant, 
so that it appears impossible to justify 
the observed increase both with geo- 
metrical effects and fluorescence (?). In 
the same range of energy the contribu- 
tion from knock-on electrons should not 


(2) P. Basst and A. M. BrANCHI: In corso 
di stampa sul Nuovo Cimento. 
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depend, according to theory, on energy. 
As a matter of fact knock-on electrons 
are few as it was possible to deduce 
comparing two series of measurements 
in one of which absorber A was close 
to counter I and another in which ab- 
sorber A was 250 cm higher. The ratio 
of the average pulse heights in the two 
different conditions is 9 = 1.00 + 0.03. 

Experiments to confirm this disagree- 


ment of our results with theory are in 
progress: such a disagreement might in- 
dicate a new method of measuring high 
energies. 


Thanks are due to Professor A. Ro- 
STAGNI of Padua University for many 
stimulating discussions and to Ing. E. 
Garti of Lab. CISE in Milan for useful 
advice on technical problems. 


IL NUOVO CIMENTO Vor, EXC Ned 1 Settembre 1952 


On the Oscillating Cylinder Viscosimeter. 


A. A. K. IBRAHIM and A. M. KABIEL 
Physics Department, Faculty of Science, Alexandria University - Alexandria (Egypt) 


(ricevuto il 14 Agosto 1952) 


When a non-elastic fluid of density o in g cm~* and coefficient of viscosity nin 
gems is placed between two co-axial cylinders, the outer of inner radius b 
(in cm) is set in small oscillations at various frequencies n (in s-1), while the inner 
of radius a (in cm) and moment of inertia J (in g em?) is suspended by a torsion 
wire of restoring constant t in g cm? s”?, then the amplitude @, of the oscillation of 
the inner cylinder is calculated from the following formula: 


(1) (3) = {(0? + D?) exp (6d/R)} | {(0® + D?) + (y2w? — D?) sin? 2Dd + 

0 
+ ([3D/2R] — yoO) sin 4Dd + 
+ (yo? + 20?) sinh? 20d + 


+ ([30/2R] + ywD) sinh 40d)} 


| 


where 9, is the amplitude inexorably imposed on the outer cylinder, 


© = {(9/8R?) + ([81/64R4] + Beyu2yıla 
D = { — (9/8R2) + ([81/64R4] + BeyuzYılz 
R = (a + b)/2 
d = (b — a)/2 


B = (non/n)"? 


and 


ie} 5 10 15 20 2a 


Fig. 1. y = Low? — t)/2aLnao? , 


being the angular velocity of the outer cylinder and Z the length of the inner 
cylinder from the surface the of liquid to the lower extremity. 

The theoretical details of equation (1) will be published later. 

From the figure we see that in the case of liquid paraffin B.P. of density 0,877 g/ml 
and viscosity 1,76 poises at 25 °C the experimental results of OLDROYD, STRAW- 
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BRIDGE and Toms (!) agree with our theoretical equation (1) for all values of 0,/8,. 
The agreement appears to justify the neglect of end effects as well as to indicate the 
limitation in our proposed theory (2). 

Also the figure shows a peak in 0/2; not exceeding unity, at a value of n almost 
exactly equal to the frequency n, of free oscillation of the wire and inner cylinder. 


It is a pleasure to thank Professor M. A. EL-SHERBINI for helpful suggestions. 


(*) J. G. OLDROYD, D. J. STRAWBRIDGE and B. A. Toms: Proc. Phys. Soc., 64 B, 55 (1951). 
(*) A. A. K. IBRAHIM and A. M. KABIEL (unpublished). 
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